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PERTURBATION THEORIES OF WEAK INTERMOLE CULAR FORGES * 

Phillip Robinson Certain 

Under the supervision of Professor Joseph 0. Hirschf elder 
ABSTRAGT 

This thesis is a contribution to the development of a workable 

exchange perturbation theory for intermolecular forces. It is divided 

into three parts. The first part develops a perturbation formalism 

for degenerate and almost degenerate energy states. The formalism is 

related to methods of Van Vleck, Kato, Bloch, Hir schf elder Kirtman, 

and Lowdin and can have a greater range of validity than the Rayleigh- 
u 

Schrodinger perturbation theory. The second part extends the formalism 
to exchange problems and leads to the Hirschfelder-Silbey perturbation 
theory. A method of solving the first order equation is developed 
which reduces the many electron equation to one and two electron 
equations. The third part applies four different perturbation 
formalisms for exchange forces to three model problems: the hydrogen 

molecule at internuclear separations R = 4, 6, 8 a ; a harmonic 
oscillator model of the hydrogen molecule ion; and che- delta- function 
model of the hydrogen molecule ion. 


•This work was supported by the National Aeronautics and Space 
- Administration Grant NGL 50-002-001. 
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FOREWORD 


ii 

V. 


This thesis, is a contribution to the development of a workable 
"exchange" perturbation theoiy for intermolecular forces, which starts 
from the isolated molecules as the unperturbed state, but which takes 
account of the full symmetry of the interacting system, including the 
Pauli exclusion principle. 

In Chapter One the mathematical difficulties associated with 
exchange perturbation theory are discussed, and it is shown that many 
different expansions of the wavefunction are possible. ‘A brief his- 
torical survey of the formalisms which have been proposed is given. 

In preparation for exchange perturbation theory, a perturbation 
formalism for degenerate and almost degenerate energy states is devel- 
oped in Chapter Two. The formalism is related to methods of Van Vleck 
Kato, Bloch, Hirschfelder , Kirtman, and Lowdin. It has the feature 
of giving the energy as the root of a secular equation. The matrix 
elements of the secular equation are assumed to be analytic in the 
perturbation parameter, but not necessarily the energy. Thus, the 
treatment can have a greater range of validity than the Rayleigh- 
Schrodinger perturbation theory. The connection of the formalism with 
Lowdin’s partitioning technique is investigated. It is shown that if 
the two methods are truncated at the same order, they give the same 
results for the energy and wavefunction, except for higher order terms 

In Chapter Three the extension of the formalism of Chapter Two 
to exchange problems is shown to be equivalent to the Hirschfelder- 
Silbey perturbation theory. It is also shown that the first order 
perturbation equation can be separated into a polarization equation. 



iii 


which introduces van der Waals correlations, and an exchange equation, 
which introduces ionic terms. In the many electron case, these equa- 
tions can be reduced to one and two electron equations. 

In addition to the Hirschfelder-Silbey (HS) procedure, exchange 
perturbation formalisms have been proposed by Eisenschitz and London, 
van der Avoird, and Hirschfelder (EL-HAV) ; Murrell and Shaw, and 
Musher and Amos (MS-MA); and others. A Sternheimer (RS) expansion is 
also possible for one and two electron problems. Assuming convergence, 
each perturbation expansion yields the exact energy and wavef unction if 
carried to infinite order. In any practical application, however, 
the energy is computed to low order only, and at this level of approxi- 
mation, the different approaches give different results. It then 
becomes of interest to ask the question: "Which formalism provides 

the best second order energy and the best expectation value of the 
hamiltonian using the wavef unction accurate through first order?" 

In Chapter Four an attmept is made to answer this question on 
the basis of model calculations. The energies of both the ground state 
( £ ) and the first excited state ( ZI u ) of the hydrogen molecule 
are computed at internuclear separations R - 4, 6, 8a Q , using the 
EL-HAV, MS-MA, HS, and RS formalisms. The results, discussed in detail 
in Chapter Four, do not show that any of the schemes is outstandingly 
better than the other three. The second order energy in the EL-HAV 
expansion, hox^ever, appears to approach one-half of the second order 
polarization (dispersion) energy at large values of R . The other 
formalisms give good values for the energy in the region of the van 
der Waals minimum for the Z^ u state, although at shorter distances 



the error increases. Also considered are a harmonic oscillator model 


and the double minimum de It a~£ unction model of the problem. 
These calculations verify the long-range behavior of the EL-HAV 
second order energy. 



ACKNOWLEDGEMENTS 


It is a pleasure to express appreciation to all those who have, 
directly and indirectly, contributed to the development of this thesis. 
Special thanks are due to 

Professor J. 0. Hirschfelder , who suggested the problem of 
exchange forces and who was a constant source of novel ideas; 

Professor S. T. Epstein, who suggested the delta-function model 
calculation and how to solve it and who was a patient and reliable 
critic; 

The taxpayers of the United States, who provided financial support 
through the National Science Foundation and the National Aeronautics 
and Space Administration; 

Drs. A. S. Dickinson, R. E. Johnson, X. J. P. O'Brien, and 
J. D. Power, who generously shared with me both their professional 
competence and personal charm; 

My wife and family, who made it all worthwhile; 

The secretarial staff of the Theoretical Chemistry Institute, 
who typed and reproduced this thesis and cheerfully performed many 


other chores. 



vi 


TABLE OF CONTENTS 

Abstract i 

Foreword ii 

Acknowledgements v 

I . INTRODUCTION 1 

II . PERTURBATION THEORY FOR DEGENERATE OR 

ALMOST DEGENERATE STATES 11 

2.1. Brief Review of Rayleigh-Schrodinger 11 

Perturbation Theory 

2.2 General Formulation 15 

2.3 Degenerate Perturbation Theory 26 

2.4 Almost Degenerate Perturbation Theory 35 

2.5 Example of Almost Degenerate Perturbation. Theory 43 

III'. EXCHANGE PERTURBATION THEORY 5-1 

3.1 The Interaction of Ground State Hydrogen- Atoms 51 

3.2 Extension to Many-electron Systems 67 

IV. MODEL EXCHANGE PERTURBATION THEORY CALCULATIONS 70 

4.1 Interaction of Two Hydrogen Atoms 71 

4.2 Harmonic Oscillator Model 95 

4.3 The Delta-Functi'on Model 99 

APPENDIX A. Notation for Chapters Two and Three 103* 

APPENDIX B. Algebraic Manipulations of Perturbation 

Equations 104 

APPENDIX C. Partitioning Technique Approach to 

Degenerate, Almost Degenerate and Exchange 
Perturbation Theory 114 



vii 


APPENDIX D. On the Modification of Exchange 

Perturbation Theory 121 

APPENDIX E. Direct Calculation of Second Order Energies 

for the Delta~Function Model 124 

APPENDIX F. ■ Other Publications of Phillip R. Certain 129 

REFERENCES 144 



I. INTRODUCTION 


The study of intermolecular forces is essential to the .under- 
standing of most of the physical .properties -of matter, such as equa- 
tions ‘of state, transport properties, scattering cross-sections, and 
-so on. For intermolecular forces’ 1 ' govern all-, of the attractive 'and 
repulsive Interactions between atoms and molecules, as opposed.’ to 
-chemical, forces, -which hold a singlfe molecule together and. prevent, it 
from separating, into its constituent atoms* Of course, there is no 

precise difference .between these two kinds of force, except in orders 

2 

of magnitude. According to Pauling, 

there .is a chemical bond between two atoms or groups of 
atoms in case that the forces acting between them are such' 
as to lead .to the formation of an aggregate .with sufficient 
stability to make it convenient -for the chemist to consider 
it as an independent molecular species., 

Intermolecular forces., on the other hand., -are sufficiently weak so 
that any aggregates that are formed, are transitory in. nature and 
easily destroyed by thermal motion. It is Best, however,, not to em- 
phasize the difference between intermolecular and chemical forces too 
much.. 

Quantum- mechanical perturbation theory , has . long'’ been used: to- 

3 

systemize the study of intermolecular forces. Although it is -con- 
venient to speak of - forces , a more basic theoretical quantity is- t-he 

electronic interaction energy, ( E ~ £ CoJ ). (The Born-Opperiheimer 
4 

separation of the nuclear and- electronic -motions is assumed'- through- 
out this, thesis-, and relativistic effects are neglected.) The. force. 


1 - 
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on a .particular molecule is the negative gradient of the interaction 
energy with respect to the nuclear coordinates of the molecule. 

The electronic energy E is obtained as an eigenvalue of the 
Schrodinger equation 


C H - £ ) =■ 0, 0.1) 

where the total hamiltonian H is a sum of t-f the electronic 
hamiltonian for the isolated molecules (i„e. separated .from each other 
by infinite -distances),, and , the -sum of all of the coulombic 

interactions between electrons .and nuclei belonging to different mole- 
cules. The energy is the electronic energy of the isolated 

it 

molecules and is an eigenvalue of the Schrodinger equation 

( - 6 <w ) = 0 } . (1 -2) 

where is a simple product of isolated molecule wavefunctions. 

The simplest application of perturbation theory to the -calculation 
of (E — € io ^) is the polarization expansion, which is a straight- 
forward application of Rayleigh-Schrodinger perturbation theory. The 
unperturbed problem is Eq. (1.2), and the theory leads to an -expansion 
of the interaction energy in powers of the perturbation = H ~ H # 

If the terms are further expanded in inverse powers of the separations 
between the molecules, each term may be given a classical or quasi- 
classical interpretation as arising -from interactions between either 
permanent or instantaneous multipole moments of the separated mole- 
cules. In this way are identified electrostatic forces, due to the 
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permanent moments of the molecules; Induction forces, due to the 
permanent moment of one molecule inducing a moment in the other; 
dispersion forces, due to the instantaneous moments of the molecules; 
and resonance, forces, which can arise if the level is degenerate, 

and which behave as electrostatic forces. Such a classification is of 
•unquestioned utility in understanding a wide variety of experimental 
results on a qualitative basis, and in providing a starting point for 
semi-empirical theories. 

The polarization expansion is applicable only when the molecules 

are well separated, however, and does not predict chemical binding 

7 

between molecules. As Heitler and London showed, binding is obtained 
with an approximate wavef unction which is an antisymmetrized product of 
isolated molecule functions, and which satisfies the Pauli exclusion 
principle; i.e. s which changes sign upon the exchange of any two electron 
labels. The added terms in the expression for ( £ - are 

called exchange forces; from the standpoint of perturbation theory, 
such terms result from a first order treatment. 

Mathematical Complications . Difficulties arise when perturbation 
theory is applied to the accurate, ab initio determination of inter - 
molecular forces in case that it is necessary to go beyond first order 
and also to take explicit account of the exclusion principle. At 
inter nuclear separations short enough so that the electronic clouds 
of the molecules overlap, the total wavefunction is required to 
have symmetry properties corresponding to the total hamiltonian H , 
which is symmetric with respect to the exchange of any pair of electrons 
between the two molecules. The most natural choice for the unperturbed 
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hamiltonian, , is the sum of the hamiltonians for the non- 

interacting molecules. However, this implies associating particular 
electrons 'With definite molecules. Hence the symmetry of H with 
respect to electron permutations is greater than that of The 

difficulties to which this gives rise are known as the exchange problem. 

Another basic difficulty is that the order of the perturbation 
terms is not uniquely defined. If A is the operator which project^ 
the component with the symmetry of the desired .total wavefunction, 
then A commutes with H . However, A does not commute -separately 
with' either R^or the perturbation but rather . 

La, h (05 ] » L h"' , a 1 d.3) 

In any conventional perturbation scheme, the left hand side of Eq. (1.3) 
is zeroth order and the right hand side is first order. This equation 
means that "order in " -is not a well defined concept, and that 

the apparent order of various terms in a perturbation expansion can be 
arbitrarily shifted 0 

Related to the nonuniqueness of order is the difficulty - of de- 
fining a symmetrized basis set for the expansion of 12? . A natural 
choice consists of the set } where the are the complete 

set of eigenfunctions of 

C h co) - e£> ) ^ 0) » o, d.4) 

where is the isolated molecule hamiltonian and ^ ^ is a simple 

A JL*°> 

product of isolated molecule wavef unctions,. The basis n^rL 1S 
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nonorthogonal, however, whereas nondegenerate eigenfunctions -of a 
hermitean hamiltonian are necessarily orthogonal. Thus it is impos- 
sible to define a single, hermitean unperturbed hamiltonian, of which 
every A ^ is an eigenfunction. In particular , since and A fail 

to commute, the * are not -eigenfunctions of H*?* 

Another lack of uniqueness is due to the fact that the complete 
set of symmetrized functions a<t are linear ly .dependent . To 
prove this, consider a function which has symmetry different from 
A , i.e. A£L- O . Since the are complete, has the 

unique expansion 

= E ^ d k 5 (1.5) 

where the C ^ are constants. By hypothesis, AtQ* vanishes, so thal 




(1.6) 


- O = £ A+? 

k ^ 

which is a statement of linear dependence of the set For a 

simple example of*Eq. (1.5) and (1.6), consider the interval — I ^ l 

and let A project the symmetric -component of any function : 

A-fcx) = . A complete set of unsymmetric functions 


on. this interval i-s 

if tQ>= x 3 


so that c k 


± r i 

L d** J x=1 




then J2 = <4>p + % + 3.4.^ + £ 


and A X. 3 = o - : A <f>. 




A4 Z ' C 


a 4 ; 


‘°>+irA4 to h 


Moreover, the linear dependence of the is non- trivial in 

case that it is possible to find a function Sh such that A tSh- O 
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and each in Eq„ (1.5) is nonvanishing. (In the previous example 

put SunTTK/iX, Then it follows that it is impossible to construct 

a linearly independent set by excluding a finite number of from 

the original set. For examp le^ can be eliminated from Eq. (1.6) 

by the expansion 

A*r - s <c«» * , 

or 

m* = r h"[^]. 

ki f 

Substituting 'Eq. (1.7) into Eq. (1.6) yields 

o = Z A4? - -r^t, 7 . 

k*\ 

Thus the set A *r 3 h ^ l 9 Is also linearly dependent. 

This means that no unique expansion of the total wavef unction of 
the form 


v - 2 n" <* 

is possible. 

The above considerations do not imply that it is impossible to 
develop a perturbation expansion for intermolecular forces that takes 
full account of symmetry., but rather that many different approaches 
are possible. 
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.storical Survey of Exchange -Perturbation Theories . The oldest 

perrurDation theory which takes the exclusion principle into account 

8 

was developed by Eisenschitz and London., The formalism was recently 

g 

put into a more modern notation by van der Avoir d. These authors 

assumed the expansion Eq„ (1„8). Of course, they recognized the lack 

of uniqueness of such an expansion, but nevertheless resolved the 

Schrodinger equation into an infinite set of perturbation equations 

which can be solved in a well defined, though arbitrary, way. Their 
1 

expression for the first order interaction energy agrees with the 
Heitler-London result,, The -second order energy is expressed in a sum 
over states form, .which when evaluated by the Unsold method^ gives 

the second order polarization energies, modified by the effects of 

\ 

exchange. 

Dalgarno and Lynn^ have introduced the Unsold approximation to 
the second order energy of a Brillouin -Wigner expansion and obtained 

the same result as Eisenschitz and London,. 

•12 

Van der Avoird has recently given an elegant wave operator 

formalism which gives the same expression for the first and second 

13 

order energy as the Eisenschitz-London expansion. Hirschfelder has 
also derived van der Avoird* s equations without the introduction of 
the wave operators. 

Musher and-Salem^ have also assumed the expansion of given 
by Eq. (1.8). These authors used a Eeenberg^ iteration technique to 
evaluate the coefficients, however, and obtained a different expression 
for the second order energy. This approach has the feature of not 
requiring H to be separated into +■ H* 11 
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Other examples In which the expansion Eq. (1.8) is used is the 

16 17 

work of Murrell, Randic and Williams and of Salem. The former 
authors assumed .that Eq. (1.8) consists of a finite number of terms 
which includes both covalent and ionic type functions. This allows 
questions of over -completeness to be avoided, and is equivalent to 
solving a finite dimensional secular equation by a perturbation expan- 
sion. Salem makes the assumption that ^ A 9 which 

does not hold for the functions defined by Eq. (1.4), but which is 
useful in assessing the significance of various terms in Eq. (1.8). 

A different type of expansion of was assumed by Murrell and 

Shaw, ^ who used a wave operator approach., and in an equivalent^ 

20 . 
treatment by Musher and Amos, who started from an infinite secular 

equation. These authors assumed that the zeroth order component of 

has proper symmetry, but that the remainder can be expanded in the 

i (O) 

unsymmetrized functions <p^ » Hence 

V = A4 ,0> + t C 

k 

The expansion in this set of functions is unique, although the higher 
order terms do not have definite symmetry properties, order by order. 

A different class of approaches to the exchange -problem involves 
different ways of defining and computing a "primitive function" whose 
projection onto the space of desired symmetry is the total wavefunction 
, in the same sense that the zeroth order function A<$ <0> is the 
projection of <£ fo) . Hirschfelder and Silbey^ propose that there 
is a physically significant primitive function whose symmetry 
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projections -correspond to all the wave functions for the family of 

states arising from a single electron configuration. Related ap- 

22 23 

proaches have been discussed by Herring, Musher and-Silbey, and 

... «. 24 

Kir t man. 

25 26 27 

'Jatisen^ (see also Byers Brown and Ritchie ) has explicitly 

constructed an operator which operates on a symmetrized function 

ASli to produce -a function in which specific electrons are assigned 

to particular atoms. With _A_ he is able to construct an unperturbed 

Vlabel free” hamiltonian whose eigenfunctions are *<C- 

28 tt 

Cor ina Ides i replaced the Schrodinger equation by a "modified 

Schrodinger equation" in which the wave function is represented by a 
vector,, each component of which represents a different assignment of 
electrons to molecules* The modified equation is then solved by a 
perturbation expansion and the true wave function is obtained as a linear 
combination of the vector components. 

Applications of Exchange Perturbation Theories . Numerical 
applications of some of the formalisms have been reported previously* 

g 

In their original paper Eisenschits and London considered the inter** 

action of ground state hydrogen atoms. More recently, the hydrogen 

29 

molecule has been considered by Alexander and Salem, who used the 

formalism of Musher and Salem, ^ and by Gorina Idesi.^ Jansen and 
30 

coworkers have applied his theory to the calculation of a wide 

31 

range of crystal properties* Van der Avoir d has treated adsorption 

32 

on metal surfaces by his method. Murrell and Shaw have computed the 

33 

interaction energy of two helium atoms; and Duijneveldt and Murrell 
have treated problems involving hydrogen bonding* 
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In all of these treatments , hox<?ever a approximations have been 
made to various terms in the perturbation expansions , thus obscuring 
the efficacy of the exact perturbation series « Also^ for most of the 
applications the exact answer to the problem is unknown^ so that it 
is difficult to judge the convergence properties of the expansions. 

Exceptions to this are the calculation of van der Avoir of the H2** 

34 

potent ialj the solution of a spin model by Johnson and Epstein^ and 

21 + 

the application of the Hirsch£elder~Silbey procedure to M 2 by 

35 

McQuarrie and Hirschfelder „ 
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II- PERTURBATION THEORY FOR DEGENERATE AND AIMOST DEGENERATE STATES 

In this Chapter, a formalism is developed for treating the effect 
of a perturbation upon a set of degenerate energy states -of an un- 
perturbed hamiltonian. The formalism,, which is related to methods of 
Van Vleck, 36 Kato, 37 Bloch, 38 Hirschfelder , 39 Kirtman,^ 0 and Lowdinf 5 
has the feature of giving the energy as the root of a secular equation. 
The matrix elements of the secular equation are assumed to be analytic 
in the perturbation parameter, but not necessarily the energy., Thus 
the treatment can have a -greater range of validity than the Rayleigh- 
Schro dinger perturbation theory. If the formalism of this .Chapter 
is expanded in powers of the perturbation parameter, however, the 

usual Rayleigh-Schrodinger theory is obtained. The connection of the 

65 

formalism with the partitioning technique is established in 
Appendix C. 

2. 1 Brief Review of Rayleigh-Schrodinger Perturbation Theory . 

Consider the solution of the Schrodinger equation 

( HU) - E^x) ) ^(X) - 0 } (2.1) 

where the hamiltonian is the sum of two terms. 

Hex') = H Co) + x H co (2.2) 

The H ( ° * is the hamiltonian for the unperturbed system, and the 

is a perturbation. The parameter X in some cases has physical 
significance (e.g. field strength), but otherwise is a formal ordering 
parameter with physical value of unity. 
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The unperturbed hamiltonian Is assumed to possess a complete set 
of eigensolutions. 


c n (o> - e<r ) =o. 


/ (®0 

The perturbed wave functions can be expanded- in the set ^ , 


•«, = 2 t ( ° } b 


^(X) 3 


where 


TV 


The fundamental assumption of Rayleigh -Schrodinger perturbation 
theory is that both and are analytic functions -of the 

parameter X ' 5 hence. 


E a Crt = £, 


Co-) 


Z x" e ^ 0 


(2.3) 


and 


where 


1 „ ''O «v V' tC vV) 

b n ul = l n -V 2* x b n ? 


CO 

= itr , 


V"t— I 


<t" - z r 4, 


(2.4) 


Co) Cm) 

^ -f 

Substituting the expansions (2.2), (2.3), and (2.4) into (2.1), and 

setting the coefficient of each power of A. to zero yields an 

— I (n) Cm - ) 

infinite set of inhomogeneous equations ror the and : 
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n-2 


( H lo> - + ( H‘“- Ep^" -0 = Z ^ 


(vj-m) 7 | (m> (2.5) 

h • 


For ft =. l , the upper limit of the sum in Eq. (2.5) is greater than 
the lower limit. The convention adopted in this thesis is that when 
this occurs,, the sum vanishes. 

The perturbation energies are determined by multiplying Eq. (2.5) 

from the. left by and integrating. Then by a series of algebraic 

r 

manipulations involving the perturbation Eqs. (2.5). it is possible 
41 

to show that 


p (2*0 


k-\ 


• < V"-° I «“>!+“> - £ f f.Z 


b=l 


( 2 . 6 ) 



Eqs. '(2.6) and (2.7) are completely general regardless of the degeneracy 
and the normalization of the wavef unction. These equations demonstrate 
the Wigner theorem that a knowledge of the wavefunction accurate 
through CPC >*) is sufficient to determine the energy through OC)?'*'*')', 
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Normalization conditions are required to complete the specif ica- 
tion of the perturbed functions. There are two types of normalization 
in common use. For full normalization^ the exact wavefunction is 

required to be normalized to unity. Substituting the expansion (2.4) 
into yields 

Z < V _kl * V°> = 0 > < 2 * 8 > 

k~o * 1 


For intermediate norma lization_, the requirement is that 

^ 4" ^ \ - ^ ^ tci) J l|> * { > OT * expanding in powers of 

% % ' % % 


i ° i. 

In practice the series (2.3) and (2.4) are computed up to some 
finite order "only. It is. convenient to define the partial sums 

Hun) - 21 X 4^ 

P V"i=£> 

(2.9) 

and 

E- (n) = 21 X ^ 

V h'O 

The radius of convergence of the expansions (2.3) and (2.4) is 5 

in general., difficult to ascertain. In many cases the ^ and 

i (m) 

D are finite for all n . even though the expansions do not 

ry 

converge for required values of X . Then it can be shown that the 


perturbation series is an asymptotic expansion of both the energy and 
the wavefunction. Since the error in truncating an asymptotic ex- 
pansion is of the order of magnitude of the last term retained , the 
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expansions (2.3) and (2.4) are expected to provide good approximations 
to the true energy and wave function of a state if the quantities 
are small. For n = g , the solution of Eq. (2.5) gives 


iO 

b n 




H 


to iV 5 > 


r (o^ 


C <0 > 

t 


( 2 . 10 ) 


so that a rough criterion for accuracy is 


| < 4 t M I H ( ° I ^ to) > I « I ^ I . (2- n > 

When the states ^ and -are degenerate or almost -degenerate,, 

the right hand side of the inequality (2.11) is - small, whereas the 

left side can be large,, so that the expansions (2.3) and (2.4) 

become unreliable. If the closely spaced levels are well separated 

from 'the remaining states, however , it is possible to modify the 

Rayleigh-Schrodinger -treatment to avoid the expansion of if 

-ja > and are degenerate or almost degenerate, while retaining the 

expansion for widely separated levels. This modification is developed 

in the next Section. An alternative approach, the 'partitioning 
65 

technique, -is discussed in Appendix C. 


2. 2 General Formulation . 

Consider the linear manifold which is spanned by £> eigen- 
functions % ^ , of the total Schrodinger equation (2.1). 

is defined to have the property that if state k in J has the 

(O) , / 

zeroth order energy 6^ , then all states K which have zeroth 
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^(o> 


order energies 6.^ that are degenerate or almost degenerate with 

, are in . The sequential labelling of the states by k is 
for convenience and does not imply that they are the £ lowest energy 
states of H . 


The 3^ constitute an orthonormal basis for manifold 


Any 


other basis in ^ may be defined by 


- L % C 




*Jl f 


(2.12) 


where the numbers are elements of a non “singular transformation. 

The basic idea in introducing the <$ L is to choose the coefficients 

K 

C Ak in such a way that; when (a) is expanded into powers cf A 
analogous to Eq„ (2„4), such terms as Eq. (2.10), with and both 
belonging to , do not occur. Thus the expansion of in powers 

of X , can be expected to have a greater range of validity than the 
expansion of *5, . 

R 

The basis functions <5^ satisfy the coupled equations 

0 0 

H = B Jl % e = ^ ^ i 

where 


(2.13) 


s' 




- £ e, c 


Jr! 


^Ak> 


and 


ft 

S-, - Z. & 


(2. 14) • 


- 1 

■y 


c 


j-k 


Given the basis functions *k' the ei gens elutions 
are recovered by solving the secular equation 
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I < <§ k l H - E. I ^ >1« 0 . t 2 - 15 ) 

Thus, Eq. (2.13) and (2.15) are equivalent to the uncoupled 
Schrodinger equations for Although <4^ and 

later he expanded in powers ^of X , it is not assumed that the roots 
of Eq. (2. 15) are analytic in X - This is* it is assumed that the 
matrix elements of H are analytic, but not necessarily the energy 
eigenvalues. An example Where such an assumption^ is valid occurs for 

the interaction of two ZA> or 2L-f> hydrogen atoms, as discussed by 
59 

Kim and Hirschfelder . See also the example in Section 2.5. - - 

At this point it is convenient to introduce a more compact nota- 
tion. Matrices will be denoted by "fat” symbols; e.g, ? 

4 = < ) , 

= (<$(0 = 

See Appendix A for a full explanation of the notation. In this nota* 
tion, Eq. (2.13) and (2.15) become 

W ^ & (2.13) 


and 


I <<*§ I H *. EL i.<$ >| - '0 


(2.15) 
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Alternative Derivation of Equation (2„13) , In order to clarify 
the meaning of Eq. (2 0 13) it is useful to give another derivation, 

based on Kirtman's treatment^ 0 ' of Van Vleck degenerate perturbation 
36 

theory. 

Consider a complete set of functions, which is split into two 
classes, and and which spans the Hilbert space of. the 

hamiltonian R (A} = H £0> + X . One class, to itfhich special 

attention is given, consists of the finite set of £ functions. 


„ t l(o ) 
Each 

hamiltonian 


4'” = C i“’, 4‘,“\ ... <#”) 

is assumed to be an eigenfunction of the unperturbed 

H l °\ 

( H (0> - 6 k ( ° 5 ) ^ 0) = 0 , \<k'« 


S 


If the -eigenvalue is degenerate, then it is assumed that all of 

4 ^o) 

, Fur ther , it 

| (o) 

is assumed that the are -orthonormal. 


<4 <t0 I «T 5 > = i, 

where k is the unit matrix. 

The remaining functions (in general, infinite in number) which 
complete the set are denoted by 

Y (0) = 


( , x “ , - 


) 
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Without loss of generality, it is assumed that 

I X <0) > = A 

and 

<r°M4 co, >= o. 


where <D is the null matrix. No other assumptions are made regarding 

■VlCO W Co) 

& ; in particular, the /Ll is not necessarily an eigenfunction 

of H to l 

4 to) 

is a basis of unperturbed functions for the ap- 
proximate calculation of 0 eigenstates ( ^ ) of the 

total hamiltonian H (», where 


E k U -<>)=• 

In general, to complete the ( 5 exactly, it is necessary 

to consider both ^ ^ and since the interaction elements 

^ < 4 >Co) I H l are non -vanishing. Then the E^ are 

roots of the infinite -dimensional secular equation 


\ trt | H-E l 4 “” > < 4 ‘°' I H- B \K Co) > 


<1t to) l H-El4 to1 > <% w iH-e.\xH 


- 0 . (2-. 16) 


Following Van Vleck^ and Kirtman, the solution of Eq. (2.16) 


is obtained by transforming the initial basis into the new set 
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4 cw) = 4 to ' + x 4 t0 + • • • + x* 4 (N \ 

and 

X ( m) = + x X co + - *• + x N X 1 * 0 . 

The functions 4 f fe ) 

are chosen so as to make the over lap and 
hamiltonian matrix elements -connecting ^ (N*) and vanish 

through terms 0C^) : 

I y-E i<$m)= O (*.*"). 

Then the leading contribution of the functions to the energies 

&k , \ ^ k ^ $ , is proportional to 

k4(«ih 0<rx 2W+a ) 

Hence the roots of the £> K secular equation., 

U<$(M) l H - E l <§CN))| = (2 * X ‘ 

are accurate through 0 ( A ^ ^ 

To obtain equations for <sf> ^ s it is convenient to 

introduce the projector onto the set 

cp - cp f = c? z 5 


(£) <$CN) = <5> (N). 
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The perturbation terms of $1(N) are given by 

which assures orthgonality between ^.Cn) and %(N) through 0Cy^)» 
The perturbation terms of <§M are determined by 


|<#tN) l Hl&ao>l= OOP*'). ( 2 . 18 ) 


The %M') 

4>(N) 

H€(n) 

Thus, Eq. 


is a .complete set of functions in the space orthogonal to 
If overlap of with all of vanishes, then 

must be expressible as a linear combination of the <3^(N) ♦ 
(2.18) is equivalent to 


WV-+ O (**+') (2.19) 

where the elements of (£ Cn) are numbers to be determined. In the 
limit N-^o o 9 Eq. (2,19) clearly becomes Eq, (2.13). 

Perturbation Expansion . Viewed as equations for the and 

€: ., Eqs.(2.13) and (2,15) are not well-defined, however, since 

all reference to the particular linear manifold^ has been lost. 

That is, Eq. (2.13) has solutions in any S-dimensional linear manifold 
defined by the exact eigenfunctions of H . Furthermore, within any 
particular linear manifold, there is an infinite number of solutions 
of Eq« (2.13) -corresponding to different choices for the coefficients 

C-jZfc 
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The lack of uniqueness in Eq. (2.13) is obviated by a perturbation 
oach to its solution* The following expansions are assumed: 


H = H to) + A H' 


<$ = Z V 


n 


OO y . -S 

£ = S ^ 

h = o 

The linear manifold^ is uniquely determined by specifying the 
zeroth order energies 


f c6 > _ r (2. 

Substituting the expansions (2*20) into (2.13) and setting the 
coefficient of each power of X equal to zero yields the infinite 
set of equations 


H4 W * 4 « £ V 


H Co, 4 <m ' 1 + H 40 ^ 0 ® Z4 t " _,V V w , n *i 


It is convenient to' define the sequence of partial sums^ 

X. .X _ 


4 CM) = s \ h 4 tn \ 

in=o > 


& Cm) = Z 4 




(2.23) 
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and the sequence of secular equations 

1<($(N) l H-E 1 CM)>| = 0 . (2.24) 

The perturbation equations allow Eq. (2.24) to be written 

0 = I <«i(w) l H-E- l \ “ 

(<§<*0 j (§CN)> (<£(N) - E) 

+ y?+'Z X n {<*| >£M) l H to l4 <w b 

v\~o 

n n+4»-h-i > 

-Z 7L <4 <10 i 

j>=o £= o 

*X 2M+, Z A h {<4 (W> ! H £ °l 4 (n; >^o 

h=o 

N-n-l f> n 

_ Z Z <4 fM+ t +, >i4 <N -f + ? ) >€ <N_ ? ) /, 

p° $=° 

This formula is derived in Appendix B. 

Thus j the total Eqs. (2.13) and (2.15) are replaced by the 
perturbation equations (2.22) and the secular equations (2.2^)^ The 
eigenvalues and eigenvectors of Eq. (2.24) for successively greater 
values of provide a sequence of eigensolutions which 5 assuming con- 
vergence., approaches the exact solutions 

Normalization conditions and the choice of the transformation C* 
must be given to fully specify ^j>^and (B The normalization 
of is completely determined by (JZ and Eq. (2.12)., assuming that 
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the original functions IB ^ are normalized- According to Eq- (2.14) 
the rows of C> are -eigenvectors of £ corresponding to the eigen- 
values . Thus, if none of the energies ^- s degenerate, the 

specification of (& uniquely fixes C- and, hence, the normalization 
of <§> . 

Relations to fix <& are given separately for degenerate and 
almost degenerate problems in the next two Sections. It can be antici- 
pated that (D .cannot be chosen in a non -arbitrary manner, since 
any set of S functions which satisfies Eq. (2.13) and (2.21) is 
sufficient to ^determine the solutions C %) • This is not to 
say that the choice of £ does not have important consequences,. A 
particular choice is C, = 4 , in which case the formalism reduces 

to a Rayleigh-Schrodinger development for each state C 0-^) 

separately. As the previous discussion has suggested, and as later 
discussions will make explicit, other choices of (L are advantageous 
in, treating .degenerate and almost degenerate problems. 

Nevertheless, in case that the Rayleigh-Schrodinger expansion of 
CEfe-, *- s possible, the present formalism, based on Eq. (2.13), 

(2.15), and any non -singular (C , is equivalent. More precisely, it 
can be shown that the roots of the secular equation (2.24) are the 
exact energies , plus terms Q(X LN * r ' tS ) , regardless of the 

choice of . Thus, if the eigenvalues and -vectors of Eq. (2.24) 
are expanded in powers of A for successive values of A/ , the 
unique asymptotic expansion of ( is obtained. 

The proof of the underlined statement is contained in the deriva- 
tion which led to Eq. (2.17). If H ^ (t^) — ) *h OC\^0 
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and is any basis for the space orthogonal fco <5(N);> then the 

interaction elements < <§> 00 | H l % CN)> are OC^0> irrespective 
of the values of * Hence if the dimension of the secular 

equation (2.24) is increased to include %(M) * original 
eigenvalues are affected by Since the inclusion of 

leads to exact energies of H 5 it follows that the roots of Eq. (2.24) 
are accurate through Q 

Eq. (2.24) can also be thought of as arising from the use of the 
variational method with the linear variational function. 

= 2 ^<-*0 \ 

where the are variational parameters. In this connection, the 

Hylleraas-Undheim theorem 5 ^ is relevant: if the roots of Eq. (2.24) 

are arranged in ascending order, they provide successive upper bounds 
to the corresponding exact energy eigenvalues of H of the same 
symmetry. Thus, in case that the states in are the lowest states 
of H corresponding to particular symmetries, the roots of Eq. (2.24) 
are upper bounds to the energies which they approximate. 

An alternative way to obtain energies which are accurate through 
0 i'K LHir ^ ) > b ^t which are not necessarily upper bounds to exact 
eigenvalues, is to solve the secular equation 

\ (£ (2N -M) — H \ — 0^ (2*26) 

* where is defined by Eq. (2*23). The roots of Eq. (2.26), are 

the exact energies > V* g x accurate through 0 C>> Zhi ^) 

regardless of the choice of (£. * This follows from Eq. (2.25) with 
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M replaced by jlN +* I • Then the sums in Eq. (2.25) are 0 
and hence may be neglected without affecting the roots except in 
Ocx aN+ *X Factoring <4>(2|J+I) I <4 (-2 N*-q> from Eq. (2.25) yields 
Eq. (2.26). 

The analog of the secular equation (2.26) in the usual Rayleigh - 

H 

Schrodinger theory is the partial sum of the perturbation energies * 
as defined by Eq. (2. 9). The secular equation (2.24) corresponds to 
the expectation value of H with the wavef unction accurate through 

0 C\ N ) • Even in the Rayleigh-Schrodinger theory^ it is difficult 
t 69 

to say which way of computing the energy is more accurate. In any 
case^ the roots of Eq. (2.24) differ from those of Eq. (2.26) by 

oa^ +z ). 

With these general results established^ it is convenient to con- 
sider separately the degenerate and almost degenerate cases in order 
to derive relations to fix the elements of (£ . 


2. 3 Degenerate Perturbation Theory . 

For a degenerate perturbation problem,, tne specrai set or states 
<sj> Co ^ is defined to be any linearly independent set which spans 
precisely the same space as the $ eigenfunctions of the S-fold 
degenerate level of the unperturbed hamiltonian 


(M to > - e<°')4^=o ? wfe<Sh 

It is convenient to choose the <sj> (c ^ to be orthonormal 

< 4 CC0 14 CO °) = iy 


(2.27) 


(2.28) 


but it is not necessary to assume any other special properties 
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Eq. (2.27) implies that in Eq. (2.22) ^ 'A • The knowl- 
edge o£ allows each term <£ Cn \ to be determined suc- 

cessively. The energies <£• ; are fixed by requiring each perturbation 
equation (2.22) to be mathematically consistent: for a general inhomo- 

geneous partial differential equation to possess a solution^ the 
inhomogeneity must be orthogonal to all solutions of the homogeneous 
equation. In the present examp le^ the homogeneous equation is (2.27)^ 
so that the consistency condition is 

<4 ( °> i 4<"-o> => £ <4 to) ' 4 Cn ~ fe) > <£ tte) 5 

te= i 

or 

n~i 

| H c,) l4‘ v '“ ) >-Z<4 fo, l4 cv '‘ B ' > ^ (2 ' 29) 

With 'the (& ^ thus defined_, the functions <2^^^ are obtained 
by solving the Eqs. (2.22). It should be noted^ however,, that the 
components <4 to, l4 ln) > are not fixed by any of these equations. 

This is analogous to the non-degenerate case where the indeterminancy 

is due to the arbitrariness of the normalization and phase of the 
52 

total wavef unction. The underlined theorem on page 24^ that the 

roots of Eq. (2.24) are invariant to (L through 3 implies 

that may be fixed in any consistent manner. A 

general approach is to leave the integrals as undeter- 

mined parameters in <Z^(n) to be fixed by minimizing the roots of the 
secular equation (2.24). This procedure is developed in more detail 
in Section 2.4 in connection with the almost degenerate case. 
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In the present Section it is convenient to fix ^<sj> | Cm?y 
bv -the full normalization condition (cf. Eq. (2.8)), 

£ <4» CB I «f> Un ' H-t0 \ = o 

2 < 4 Cte) 1 4 czn_ fe) > » - i <4‘ w 1 4 C "^ 

k=o 

This choice of normalization has the properties that 

i) {<5i<£> = & ; this follows simply by expanding 

in powers of A and using Eq. (2.28) and (2.30); 

ii) is a hermit ean matrix; this follows from i and 
Eq. (2.13), i.e„, 

(<SI H l^»> = <<4 13T>€ - <6 . (2. 3i) 

The first member of Eq. (2.31) is her mite an; hence the last is also. 
In Appendix B, it is shown that the perturbation terms are also 

hermitean. Intermediate normalization, A ^ leads to a 

nonhermitean (E ; 

iii) the knowledge of <2? accurate through O is 

sufficient to compute accurate through O + * This 

is the analog- of the Wigner theorem and follows from the formulas 
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gCZn+D = | H £, M4 Cn> > 

_t £ 21 { In+Ji+i-fe) ^<rwO>+ <^,(*-*>| 

k'i jM 

h zl . _ 

jp ^cE £te ^<4 4n ^ r " k) i4 CM) >-<4 c ^ ) l4 c ^ H ^>} 

2 Ma i=o 

w 

^tZ H {<<4‘ M ^ + '“ fe) I4 tn ‘^> - <4 CM ^i4 £vi ^'" b >>}e (fe > 


k-2. jl=6 


(2.32) 


where 


r *? 
KMt. 

L o 


4? 5 k even 


-* , k. odd; 

<£ <znl =i{<4 £ "-' ) iH w l4 £n >> + <4*)| (4 £i ’i 4^-o>j 


h k-1 


? S E{« <W <4 < ^ 4 W' , ' ht *>> + < 4 C»-fe-wo |«b £ n-«> <g c te>t 

b=i 0=.n \ f J 


fc=l £=0 

*-. m 


~\Yj Tj <<4 j, 4 Cw --« ) )l 

b=i l*o J 

HZ X{<4 (n ' jZ> l4 tr, ' te+J!) >-< < 4 ,tn “ fe ' M) 14 <vi ' s; >} s°° 


fc=i 


where 


r fe^z 

W ’ ti’ 


(2.33) 


%~ Z , k even 
fc odd. 

For V\ — 0 } I , -.recall the convention stated following Eq. (2.5). 
These formulas are derived in Appendix B by algebraic manipulations 

of EqSe (2.22),, (2„29) ; and (2*30) in a manner qompletely analogous 
52 

to nondegenerate theory. 
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It is at first sight surprising that Eq. (2.32) and (2.33) 
depend on a particular choice of normalization, in view of the proof 
on page 24 that is sufficient to compute the energies , 

accurate through O regardless of the normalization. There 

is no conflict in the two results, however, since, for example, the 
off-diagonal elements of do not contribute to the roots of 

Eq. (2,24) through OCX 11 ***). 

Solution of the Perturbation Equations . The n-th order perturba- 
tion equation (2.22) may be solved exactly in terms of the eigen- 
functions of the unperturbed hamiltonian H . This is facilitated 
by the introduction of the reduced resolvent. 



where the prime on the summation means that all the members in 

are to be excluded from the sum. Otherwise, the sum is over a complete 

set of eigenstates of . It may be readily shown that 

R <0) 4 fo) = O, 

(H w - ) R Cc) = I4' 0, X4‘ 0) J ~ / 

In terms of ^ the solutions of the perturbation equations which 
satisfy the normalization conditions (2.30) may be written 

= R <c0 - Xs < 4 C "'~ l0< £ (W ^ + 4 <0) (2.34) 

te=! 3 



where the matrix is determined by requiring that satisfy 

the normalization conditions Eq. (2.30). The explicit formula for 
&4 ^ is, in terms of lower order functions., 

= £ <40*) | ^CZn + i-jO^ 

3 

and (2,. 35 ) 


W <2h) = 4<4‘" ) !4 <n, > -g <4 m) 1 

That defined by Eq. (2.34) actually solves -the perturbation 

equation (2.22) and the normalization condition (2.30) is easily 
verified by substitution. 

Except for simple examples, the expressions given above are of 
formal interest only, since the summations over excited states cannot 
be evaluated. In general, however, it is possible to obtain varia- 
tional approximations to For example, the first-order 

functions may be determined by finding the stationary points of the 
functionals 


- <c 


l‘ 0) -e (o> l4fe 0 > + < 


b« 

r k 


io) x 

h y 



where 


LCD 


is a trial function corresponding to the exact function 
which belongs to . It is evident that if S *3^ 

W 'l/ /V 

vanishes for arbitrary variations i n ^^ ien 

satisfies Eq. (2.22) for VI s l « This is analogous to the Hylleraas 
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variational method in non-degenerate Rayleigh-Schrodinger perturbation 
theory. 58 ’ 52 

Calculation of the Energy . Given the solutions of the first H 
perturbation equations (2.22), energies accurate through O ( ) 
are obtained by solving either the secular equation (2.24) or the 
secular equation (2.26). 

As discussed previously, in case that the states in are t 1 
lowest energy states of H , the best approximations to the exact 
energies obtainable from Eq. (2.24) result from minimizing the 

roots with respect to the integrals It was also shown, 

however, that these roots differ by O (X****^) from the roots 
obtained using the full normalization condition Eq. (2.30). Then in 
particular cases, Eq . (2.25) simplifies considerably. For | and 

full normalization, the secular equation (2.24) becomes 

o = |'< <§(0 I h - e I £0>! = 


(It- A 2 (<£ (o, +- Acs <0 -E ) 

+ A 2 (k m A 3 & (3) 

+ i A 3 I <4 0) * 4 <0 >D 


(2.36) 


The root's of Eq . (2.36) are accurate through O ( A 3 }. 

Alternatively, the energy may be obtained as a root of Eq*. (2.26), 
which is easier to apply than Eq. (2.24) because the matrix elements 
are simpler and 5 appears only on the diagonal. The matrix 
may be diagonalized by a unitary transformation to give 



33 


eigenvalues E^ , accurate through & (X** 1 ***), and corresponding 
eigenfunctions (^0 “ CN)^” , accurate through 0(x N ). 

For jS|sj , i-t may be verified explicitly that the roots of 
Eq. (2.24) differ by from the roots of Eq. (2.26). 

Factoring (A + f<<4> a >l4°>>) from Eq. (2.36) yields 

0 - I <<§CJ) I H- E I <$(.)> I - 



« < ° , + A<£°>+.'aV* , + > 3 c f3) - E 

+ [> £,) 3 <4 u> i4 w >] + 0tiV 


or 


| g (s) - E x 3 jV 0 , <4‘°l4 t0 >] + 0-W)l - 0 . 

When is diagonal, the diagonal elements of the commutator 

vanish, so that the commutator does not contribute to the eigenvalues 
through O(^). Hence the roots of Eq. (2.36) agree with the eigen- 
values of (£(3) through 

Connections with other Degenerate Perturbation Formalisms . The 
first order perturbation equation (2.22) is equivalent to the first 

o/: 

order equation of Van Vleck degenerate perturbation theory and of 
37 38 

the treatments of Kato and Block. For example, in Kato’s method, 
one solves 


K4 £o) ! h, - e k J 4 Co) >! = o. 
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where 


H a = € <e) + X H (,) + X 2 H 10 R fo> H 0> +- 

and 


K a = ! - A* H <0 (R t0> f H i0 +■ ... 5 

which is seen to be equivalent to Eq. (2.36) through OC A a ) Thus 
the present formalism represents an extension of these, methods to 


arbitrary order. 

If the secular equation (2.24) or (2.26) is solved by expansion 
in powers of X s the usual Rayleigh-Schrodinger perturbation theory 


for the states in 


is obtained. One must then determine the correct 


zeroth order wavefunctions , which depend on the order in which the 
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degeneracy is lifted. Hirschfelder has given a very thorough dis- 
cussion of the complicated equations which result. The present 
formalism is much simpler 'because it does not contemplate the expan- 
sion of the secular equation, which may be solved by other procedures^ 
The connection of the present formalism with the partitioning 
technique for degenerate problems is considered' in Appendix C. It is 
shown that the secular equation (2.26) results from a factorization 
of the corresponding partitioned secular equation which leaves the 
roots invariant through For W>1 ^ Eq. (2.26) is 

simpler ta solve because the partitioned secular equation contains 
the energy E in a nonlinear fashion; 
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2 . 4 Almost-Degenerate Perturbation Theory . 

For an alraost-degenerate perturbation problem, the special set 
of states is defined to be $ states C 7 of t ^ ie 

unperturbed hamiltonian, 


( H 


(o) 


€ 




K“ - o. 


■which are almost degenerate in the sense that the ratios (see Eq..(2.10)) 


<4^1 h (i M4v) 

er - £& 


, k,k' 


(o') 


are large. For simplicity, it is assumed that each of the levels £j! 
is itself non-degenerate and that < 4 to > l af> <°>> = i , 

Again the general equations to be solved are Eqs. (2.22) and 
(2.24) or (2.26). These equations are similar to the degenerate 
case, except that now is diagonal but it is not a multiple of 

the unit matrix. Furthermore, the perturbation energies are not 
fixed by requirements of mathematical consistency. To see this, 
consider the first order equation (2.22) for the component : 

( h‘*> - e‘°’ )4>‘° + H u> = 2 $ #> * (2 ” 37) 


r 

0) 


The only general condition to fix the is that the inhomogeneity 


must be orthogonal to the eigenfunctions of 


i(<*> ... - <©> 

t with energy 


i to) 

By assumption, this consists of the single function <p b ; hence, 

.<«>* 


multiplying Eq. (2.37) by 


and integrating yields 
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That is, only the diagonal elements of €*** (and, in general, ) 
are determined. 

Thus the off-diagonal elements of <& are not fixed by the 
formalism. Each choice of these numbers results in a different ex- 
pansion of and . Although the roots of the secular equation 

(2.21) obtained by different choices differ only by "higher-order 
terms", it is of interest to examine various ways of fixing the off- 
diagonal elements of (S . 
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The DE-FOP- VIM Formalism . A general method of handling the 

indeterminancy in <5* is to solve for each <jj>^ as an implicit 

_ Cn) 

function of the off-diagonal elements fcy ^ . Then these quantities 

can be determined by minimizing the roots of the secular equation 
(2.24) with respect to variations in the * 

To gain insight into this procedure, consider the solution of 
Eq. (2.37) with j arbitrary. The general solution may 

be written 
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+ H 0> 4£ 0> 





<< S * 7 i 4>r> = ° , 

and is an arbitrary normalization constant. Then the 

S-dimensional variational basis set 

^k c,) = V + *4k 7 > 


which leads to the SxS secular equation (2.36) with replaced 

by <$ C0 , is clearly equivalent to the 2S -dimensional basis 




, Ufe<S\ 


This leads to the 2Sx2S secular equation 


0 


<4‘°M h-e I4 to) > <4 f0 MH-El4 cl) > 

<4 to \ h-e\4 £0) > <4°M h-e.\4 U) > 


' (2.39) 


S roots of Eq. (.2.39) are identical to the roots of Eq. (2.36) 

if in the latter equation the ^ and Ct^ are varied freely 

and independently to a stationary point. 
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The Eq. (2.39) appears also in Hirschfelder ' s DE-VOP-VIM 
39 

formalism for degenerate and almos t-degenerate perturbation problems 
It is easily seen that the extension to fsf-th order of this method of 
choosing & leads to an g secular equation, j£> 

roots of which give the desired energies accurate through 

The DE-FOP-VIM method is simplest to apply in case that the 
states fe in are the lowest energy states of H corresponding to 
particular symmetries. Then, by Hylleraas-Undheim theorem, the 0 
smallest roots, of the DE-FOP-VIM secular equation can be identified 
as the approximations to the corresponding exact energies. In case 
that there are lower energy states of H than the states in ^ , 

the interpretation of the roots of the DE-FOP-VIM secular equation 
can become ambiguous. Clearly, the DE-FOP-VIM procedure cannot be 
applied to the secular equation (2.26) since this equation cannot be 
derived from the variational method. 


The Kirtman Formalism. 
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Another method of choosing the off- 


diagonal elements of is by analogy with the degenerate case. The 
full normalization conditions (2.30) uniquely fix all elements of dr 
just as in the degenerate case. Furthermore, in the present case dr 
is hermitean and the formulas (2.32) and (2.33) for continue to 


hold if the wavefunction is required to satisfy Eq. (2.30). 
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This method of fixing <§■ was first discussed by Kirtman in his 
extension of Van Vleck .degenerate perturbation theory to the almost 
degenerate case. Kirtman considered in detail the calculation of the 
energy through third order, which in his formalism is given as a 
root of the secular equation 
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0 = 


*5* -E 

+ <4^n h u> i«p> v <4 Ui l 14 ®>> ^ 


(2.40)' 


^{<4°MH^l4 c,) >“t^W , M4° i >-i:<cs} cu H c, >>^ £ ' > ]' 


The first order equation (2,22) and Eqs. (2.32) and (2.33) can be used 
to simplify Eq. (2.40) to 


\ <£ ie> + *<£ 0> * A 2 £ ( * 5 + > 3 <# J - E. 1 « G a (2.41) 

which is Eq. (2.26) for NS « 5 « 

Hence the Kirtman formalism will denote the treatment of almost 
degenerate perturbation problems by the secular equation (2.26) with 
the full normalization conditions (2.30). 

The modified Kirtman formalism will be used. to label the treatment 
based on the secular equation (2.24) with the full normalization 
conditions (2.30). 

The advantage- of the Kirtman or modified Kirtman formalism is 
that the energy is obtained as a root of an SxS secular equation, 
rather than the CN-h)£ 2 X (N* 0*5 secular equation for the DE-FOP- 
VIM formalism. The roots of the DE-FOP-VIM secular equation are 
necessarily more accurate than the corresponding roots of the modified 
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Kirtman secular equation in case that the lowest energy states of ff 
are being treated. In any case the corresponding roots of the three 
different equations differ by . Thus, if the almost 

degenerate block of states is well separated from the remaining unper- 
turbed states, the energies obtained by the three methods differ by 
terms which are, by hypothesis, negligible. 

Transformation to an Exactly Degenerate Problem. The lack of 
uniqueness of the off-diagonal elements of can be avoided by 
defining a new split of H into an unperturbed hamiltonian and* a 
perturbation, such that the unperturbed limit is exactly degenerate. 

In some cases there is a natural choice for the new unperturbed 
hamiltonian and eigenfunctions. In general it is possible to define 

H*" = H <0) + !4 fo) > [c (0) - <g Co5 ] <4"l, 

and 

H t0 = + I4 w >r* <,9 -i w l<4«l, 

where 

_ gico ^ 

and e ^ is an average unperturbed energy ; e . g . ^ 

g“> = S“' Z 

fe-l 
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Then all of the states <$j> ^ are exactly degenerate with respect to 
-6^)4 <0) = ©, (2.42) 

so that the formalism of Section 2.3 can be applied. This method of 
handling almost-degeneracy is suggested in most textbook's . 

Let a bar denote perturbation terms in the expansion of ^ and 
<& , based on Eq. (2.42); 


— <^> ? •£» » « * 

and 

(£ = c£ ft> . + <$ co +- (B tz) + *** > 

Through first order, £ Co) does not contribute to (£ since, by 
Eq. (2.32), 

i<‘> =. <4 <0) I = A£ co + <B Ce/> — 

so that 

+- i c,) = +• u (l, t 

With the normalization (2.30), ^<4* lc ^ ^ — (0 # the fi rs t 

order equation ( 2 . 22 ) can be written 


(h ( °>- 6^)4 0) + AH‘°4 (e> = >4 to) C t,) 
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Comparison with Eq. (2.37) reveals that enters into 

and hence ( h , in a non- trivial way, so that the latter quantities 

are not simply related to <4> ^ and <£ . For this reason, 

development along these lines is not considered further here. 

Solution- of the Perturbation Equations . A formal solution to 
Eq. (2.22) may be given in terms of the resolvents 


R 


to) 

& 


r / 

• < a* (o) Co) 

A € h 


9 


7 


where again the prime means that no state belonging to <sj>^ is in- 
cluded in the sum. This set of resolvents has the properties 


R 


(o) 

k 


Co) -s 


( H'“- e" ) r" = K“K<n- l. 

The solutions of the perturbation equations which satisfy the normaliza- 
tion conditions (2.30) may be written 



+*“ = Rr h o> 4 >r , 


o e e> 


- rH h u K"'-£ 

m- i £=* jgbj 


^ 5 
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where M w is given by Eq„ (2.35). As in the degenerate case, 
variational principles may be given to allow the approximation of 
in practical calculations.^ 0 

The relationship among the DE-FOP- VIM, Kirtman, and modified 
Kirtman formalisms is best understood by means of an example. 

2 . 5 Example of Almost-degenera te Perturbation Theory . 

6 1 

The system of two coupled simple harmonic oscillators described 
by the hamiltonian 

= ~2 + if- ) + * + i ( i - a ) 2 ^+ 



provides a nice illustration of the formalism developed in the pre- 
ceding sections. The Schrodinger equation for this case is exactly 
solvable, and the eig'enenergies are 


E 




C 4) co ^ + Cm ) wjl- 


(2.43) 


where n ^ m - O ? 1 ? e . and 


60* 



\ +• ± y 4 a 2 * +- 


The theory will be applied to the almost degenerate hsl, fcvy=- <9 
and K-O, wi=M states. By squaring the sum and difference of 
and , and then taking the square root of the result, one can derive 
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Perturbation theory corresponds to expanding the radicals in 
Eq. (2.44). Assuming that A allows the inner radical to 

be expanded in' powers of X . For the upper sign, this assumption 
also permits the expansion of the outer radical. To expand the outer 
radical for the lower sign, however, requires a knowledge of the rela- 
tive magnitude of X , the perturbation, and A > the splitting 
between the frequencies. This is the characteristic feature of almost 
degenerate perturbation problems. At this point, the most general 
expansions possible are 



where the upper sign corresponds to \0; the lower, to 0\ 0 

The Rayleigh-Schrodinger treatment of this problem corresponds 

to expanding in powers of X » which clearly converges 

3 

only for X < A Nevertheless, through order A , the energies 

are given by 


E lo teg) 

oi 



(2.45) 
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For the application of the formalism developed in the last 
section, the unperturbed hamiltonian is taken to be 

= -i (& *^.)*^'*^* % +iO-4)Y a 

and the perturbation 

H t0 = • 

The eigensolutions of the unperturbed problem are 

= (H-A) -v- C»-A) } 

and 

4*n m ~ 

where is a simple harmonic oscillator eigenfunction. 

The theory t-” 1 ! be applied to the states H — l , m— o and 
H — O , im — f Hence, 
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and the solution of the first’ order equation (2.22) yields 

4 C0 * - (4^ , 4,‘r ). 


a 


CD 




{ 


2-i V-a 2 - 

O 


<k 




X 


4(4~A^) 

O 


o 


4-Cl~ A?) 


(& 


C3) _ 


V^tt-A’-D^ 2 - ° 


To this order j the energy of the Kirtman formalism is given by 
the roots of the secular equation (2.41) which yields 


.jo (k) 2 


_Xl 

4(»-A a 5 




^3 / _ r i2 

(2.46) 


Alternatively the modified Kirtman secular equation (2.36) 
reduces for this case to 


|(4+ A z <4 f 44 U5 >)(£‘ o7 +^ C0 -- E.) + ^ z < £ + A s € C3 4= O. 
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which has the roots 


E (0 CH K) = 2 

oi 


XV 


?'0 a A 4 )>•>■*• 


(2.47) 




a 3 ' 


ZClr-A 2 )'^ + ZC»-A*)^t€(»-A*)+>*J 


Since the 10 and 01 states are the lowest energy states with their 
symmetry, the roots (2.47) are upper bounds to the exact energies. 
Clearly, the roots (2.47) differ from the roots (2.46) by terms of 
order . 

For the application of DE-FOP- VIM to this example, the secular 
equation to be solved is Eq. (2.39) which' can be shown to be -equi- 
valent to 


(H - E $ | = <D 


where 


IH* 








< 2 ) 


2^ Co Vx<e £O }< 4 C0 14 Cu > 

+ X <4‘0 [.<£« + ^ U) l 
->, 2 -<£ c * 5 - t - 
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Substituting in the proper quantities -and transforming to an -ortho- 
normal set of functions yields 

0 = I $”' /z H - £i\~ 


Si 

2 + 6 ~ E 2(i-a‘) v * 


O 


au-A 2 -)’^ 




2“ A *“ £ 




o 

X 


•+ A~ E 
A 


C 

i(T--A a )' /s - ( i-a*)^ 4-a-E 


(2.48) 


££ this secular equation is now solved by the- matrix version' of the 
same formalism that was used to obtain Eq. (2.46) and-(2;47) it can 
be verified that the identical result is obtained, except: for 
"higher-order terms". Thus if the perturbation expansions (2.20) 
are rapidly converging, the 4x4 secular equation, (2.48) yields only 
a small improvement in the accuracy of the energy over the 2x2 
secular equations. 

The various perturbation treatments are -compared numerically 
in Table 2.1 for the splitting of the unperturbed frequencies - 

As expected, the Rayleigh -Sch rod inger result, Eq. (2.45), is very 
inaccurate for X> & \ Of hhe other methods, the Kirtman formalism 
provides the- best approximation, although the roots are not neces- 
sarily upper bounds to the exact energies. The DE-FOP- VIM results 
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must be more accurate than the modified Kirtman results since the 
states under consideration are the lowest energy states of the same 
symmetry. 



TABLE 2.1: Errors' In • Perturbation Energies of Coupled Harmonic Oscillators. 






^('approximate') - 

- E(exact') 





Rayleighr 


Modified 



X 

E(exact) 

Schrodinger 

Kirtman 

Kirtman 

DE-FOP-VIM 



Eq. 2.43 

Eq. 2.45 

J Eq. 2 .46 

Eq. 2.47 

•Eq. 2.48 


0.00000 

2.10000 

0 

0 

0 

0 


1.90000 

0 

0 

0 

0 


0.03980 

2.10158 

0.00002 

0.00000 

0.00000 

0.00000 

1.89762 

-0.00002 

0.00000 

0.00000 

0.00000 


0.07960 

2.10611 

0.00029 

0.00000 

0.00000 

0.00000 


1.89068 

-0.00028 

0.00000 

0.00000 

0.00000 


0.10000 

2.10941 

1.88552 

0.00069 

-0.00067 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 


0.11940 

2.11306 

0.00134 

0.00002 

0.00002 

0.00002 


1.87971 

-0.00131 

0.00002 ' 

0.00002 

0.00002 


0.15920 

2.12177 

0.00383 

0.00005 

0.00007 ' 

0.00006 


1.86532 

-0.00372 

0.00006 

0.00008 

0.00007 


0.19900 

2.13166 

0.00834 

0.00012 

0.00017 

0.00014 


1 . 84809 

-0.00809 

0.00014 

0.00019 

0.00017 

x>&. 

0.23880 

2.14222 

0.01538 

0.00025 

0.00034 

0.00028 

1.82844 

-0.01484 

0.00029 

0.00040 

0.00036 


0.27860 

2.15311 

0.02529 

0.00045 

0.00063 

0.00052 


1.80668 

-0.02428 

0.00056 

0.00076 

0.00069 


0.31840 

2.16405 

0.03835 

0.00077 

0.. 00107 

0.00087 


1.78300 

-0.03660 

0.00098 

0.00132 

0.00120 


0.35820 

2.17483 

0.05477 

0.00124 

0.00171 

0.00139 


1.75753 

-0.05193 

0.00161 

0.00217 

0.00199 


* A = 0.1 ; For each X 3 the upper entry gives results for n=1 ..m=0 state: che lower entry 3/ 
for the n s =0,m= ! l state. 




III. EXCHANGE PERTURBATION THEORIES 


The mathematical difficulties associated with the development 
of exchange perturbation theory have been discussed 'in Chapter One-. 

It was shown that the concept of order of perturbation terms is not 
rigorous and that many different expansions of the wavefunction are 
possible.' In Chapter Two the perturbation formalism for degenerate 
and almost degenerate problems was developed with the intention of 
using it as a guide for exchange perturbation theory in the present 
Chapter. The essential idea .is to construct, by perturbation methods , 
a finite basis set for the expansion of the wavefunction. The basis 
is defined to have small hamiltonian matrix elements coupling the 
basis with any function which is orthogonal to the basis. The energy 
and properly symmetrized wavefunction are obtained by diagonalizing 
the hamiltonian in the finite basis. The advantage of this procedure 
is- that the basis functions are not required to have the symmetry of 
the total wavefunction, but simply form a basis for a reducible rep- 
resentation of the symmetry group of the hamiltonian. This approach 
to exchange perturbation theory leads to the Hirschfelder-Silbey 
formalism. 

To avoid the added complications which arise when even the 
separated atom wavef unctions are unknown, the treatment developed in 
this Chapter is limited to the interaction of hydrogen atoms. 

3 . 1' The Interaction of Ground- State Hydrogen Atoms. 

Consider the interaction of two ground’ state hydrogen atoms, a 
and b, separated by a distance R. There are two molecular states 
which arise from this separated-atom state: ( $ E ^ ) 3 the *2 
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ground • state; and C M E ^ U ) , the ^21^* repulsive excited 

state. The exact spatial wavefunctions for these states satisfy 

( H - £ e ) = €? s 

and 

Aj- ® , i = 3 eir «•» 


The hamiltonian H is defined by 


O 


X 


^ 2 , 



and 


is a symmetry projector defined by 


= 


2, 




permute 


The upper sign is for symmetry; the lower, for U 
the electron labels, *1^ permutes the nuclei labels, CT reflects 
the electronic coordinates across any plane which includes the inter' 
nuclear axis, and is the projection operator for the M s 0 

eigenfunction of the total azimuthal orbital angular momentum 



jr M — of& 

n>o M 


A is seen to be a compound projector for the symmetry group of 
H The first factor in Aj projects onto the space of either 
singlet (upper sign) or triplet (lower sign) functions; the second. 
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onto either gerade or ungerade functions; the third, onto + functions 
and the last factor projects onto the space of ZI functions. 

Zeroth Order Primitive Functions . When 'the atoms are infinitely 
separated (r-*« 0 , an exact wavefunction for the system is 

= a 0 o) b 0 (2) > 

where is a Is orbital centered about nucleus a, and b 0 is a Is 

orbital centered about nucleus b. 'The subscript "1" indicates the 
arbitrary assignment of electron 1 to atom a and electron 2 to atom b 
Because of the physical indistinguishability of electrons, ° * 

is degenerate with the configuration "2", 

fz 0) = fiz = k(0 ( 2-3 

in the >sense that 

<4,“"’ i h 1 4,“’) = <<triHHr>. 


This is a different type of degeneracy from that discussed in 
Chapter Two in that ^ ^ are not degenerate eigenfunctions 

of the same unperturbed hamiltonian. In particular, 

( h, m - 6 ( °n4r = o, (3.D 

where 


H/ 0) * c„ 


6 <0) = 2e a ; o. 2 ) 
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and 



(3.3) 


and is the energy of a ground state hydrogen atom (= - % Hartree) , 

but 

( ^ 

where 


^ = ^<0 + iaCz) . (3-4) 

Eq. (3.2) and (3.4) give H, W ~ + ^ ^ J • 

4 (o ) f fo j 

J and are "primitive 

functions", i.e. s they are a basis for reducible representation of the 
symmetry group of R . Since <# 0) and <t>? dp not have pure 
symmetry, there is. no Rayleigh-Schrodinger expansion of 
which starts from ^ as the unperturbed states and pro- 

duces a wavefunction that has the symmetry of the exact functions, 
order by order . 

43 A A 

I-Ieitler-London Wavefunctions . The projectors allow 

the construction of zeroth order functions which do have the symmetry 
of the exact wavefunction. Putting 




Aj <£, to) a 0 o)b a <2.) i-kb(0 a a (z) 

<A 3 4riA 3 c>' / *" << o) 1 fh i ,fe 
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(3.5) 

A u 4| <c> _ b 0 U) - hpco OcU) . 

<A„4,“IA„4, W >'' 1 ' -^ < I < - Fi a I +, < “>"* 3 ' 

yields the Heitler-London approximation to the wavefunctions . In 
perturbation theory, the energy through first order is simply the 
expectation value of H computed with the zeroth order wavefunction. 
Thus 

i E to> + ’e £0 = < <°> | h I i ^ co) >. 

Eq. (3.1) s (3.3), and (3.5) yield 

= 4 £ ^ o) = £ 


and 

u- ip (0) = 


and 


l £ <o = <^/ e) Uh-h^)(i±^)|^ 0> > 

<4/ D) I i ± p n | 4«>> 


(3.6) 


v~ * 

i ± C 
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where 


V oe — // [0oC0b o (zH [r;^ + 

V ol “ IK “> k<0 a * u > 4 « [ -fc “ i - vi, + i] ar ( 

5 0o = /«.<!> ^0-> q/r, # (3.7) 


43 53 

The integrals in Eq. (3.7) are all well known. ’ 

It is not possible to proceed in a straightforward manner to 
higher order, however, since s are not eigenfunctions 

of or H a • 

• Exact Primitive Functions . On the other hand, given the exacc 
spatial wavefiinctions and 3 it is possible to construct 

two exact primitive functions ^ j and which have precisely the 


same transformation properties as <$> ' and 

<5 - T> 

j-o 1 1-5 


That is, 


(3.8) 


1 i'L 


where 


= v & , 


(3.9) 
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and 

<i = ( , 'ij , 

5? * (’ti , 14 -Q) ,. 

Eq. (3.8) and the requirement that and be norma l ized 

implies 


C 




As R increases without limit, & approaches 

0 o •=■ Q/fctam. ^>© I?** 4^*|3 ^and approaches +r 

Exact primitive functions have been discussed previously by 
Herring, Hirschf elder and Silbey,^ and Musher.^ Intuitively 

one expects to have ‘electron 1 localized about atom a and 

electron 2 about atom b. Herring, who call's ^ the* "home base" 
function, enforces this intuition by asserting a set of -auxiliary con- 
ditions that ^ must satisfy., such as the requirement that 
approach when 'Q. { — O or ^ * (f is not obvious 

that this condition can be satisfied.) By ^contr.ast, Hirschfelder and 
Silbey enforce their intuition only in zeroth order by setting up -a 
perturbation sequence for «5j>j which starts from . This is 

the approach followed here. 

The set of coupled equations .satisfied by is 


5 


(3.10) 
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where 


d « <£"' E £ 


an * 


o 


o 

a E_ 


It is convenient to define two energies, the Coulomb energy 

E ( Ccrw, I ) = "£ ( % + W E. ) ^ 

and the exchange energy 

E (exek) * T ( * E - U E ) . 


Then, 


E (ccml) 
EL (exck^ 


E (e&cW) 
E (<W) 


Given s , 
obtained by projection 


, the exact wavefunctions ^ are 


l£> 


Aj dl i 

<A^$, ! A 


The energies are given by 


= <*$1 H I >. 
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Perturbation Expansion . The formalism developed in- Chapter 
Two suggests itself as a means of solving Eq. (3.10). It must be 
extended somewhat, however, since the previous analysis assumed that 
all of the unperturbed functions are eigenfunctions of the same 
hamiltonian. This assumption is no longer valid in exchange per- 
turbation problems. 

The zeroth order equations corresponding to Eq. (3.10) are 
defined to be 

( nr - = ■*, 

and (3.11) 

( H - 6 = D. 

Corresponding to ^ ° * and are the perturbations 

H,‘" = H - N,*"’ 

and 

Hi 0 - H “ H‘ 0> , 

respectively. * The superscripts are used in a loose sense to denote 
the order of perturbation. As discussed in Chapter One, "order" is 
not a well defined concept, however, since 


- [>,. ,h,“] * ( 


r b 


Eq. (3.12) states that a "zeroth order” quantity is equal to a "first 
order" quantity, so that the apparent order of terms can be arbitrarily 


shifted. 
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Nevertheless, Eq. (3.10) can be expanded into different orders 
in a consistent, though arbitrary, manner by defining 


d = 6 <0) 4 + £ V £ 




n-B 


and 


(3.13) 


# = 4“’ + Z, i" <tf"\ 




where X is a formal ordering parameter whose physical value is zero. 
The hamiltonian can be resolved in two ways, either as A 


, £©> 


u> 


*0 w ^ ** s °8 
Only for X- I does H 8 — H, 


or as H a « +■ X H<g_ 

The convention adopted here is that when H acts on , it is 

replaced by H ^ ~ 0 or £> . 

Substituting the expansions (3.13) into Eq. (3.10) and equating 
the coefficient of each power of X to zero yields, in addition to 
Eq. (3.11), for Sand 2, 


0 C-s“’K“ + 


z 


Z +re£' 


(3.14) 


c H k - e )\ + H k 4k - 4 fr % 6 Ak 


n^-L 


Given the solution of the first 
the energy is given by 


Vh=l J0- 
perturbation equations (3.14) 


(Aj j^fro] H j Aid>, (m)> 

(A i 5>,oo \ A t -3L, cn)> 
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where 

<£,<«) = L 4, <n \ 

8 h~o 


The justification for this procedure -has been developed in Chapter 
Two. If X is any wavefunction which is orthogonal to and '5^) 

then the hamiltonian matrix elements ( Ai I H I Ai xy 

vanish through O ( . Hence the addition of X to the basis 

<i(N) affects the energy only by OCx 1N * & ) terms . 

Equivalence to Hirschfelder-Silbey Perturbation Theory . It is 
easily 'verified that Eq. (3.10) is identical to the equation solved 
in the Hirschfelder-Silbey exchange perturbation theory for the inter- 
action of ground' state hydrogen atoms. Nevertheless, the solution of 
first order Eq. (3.14) is discussed in detail below in order (1) to 
demonstrate that, in analogy with almost degenerate perturbation theory, 
the off-diagonal elements of & are not completely fixed by the 
formalism; and (2) to develop a method of solving Eq. (3.14) which 
reduces the corresponding equation for many-electron systems to one 


and two electron equations. 


Solution of the First Order Equation . Eq. (3.8) implies that 
■ 

Eq. (3.14), which may be written 


^so that it is sufficient to consider .only ^ ~ j) 


e m (e*d04l o> « 

(3.15) 

Multiplying Eq. (3.15) from the left by and integrating yields 
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V„ 


gf ° (Coa!) 4* S PO E (exck) 


O ) 


(3.16) 


where is defined by Eq. (3,7). Thus, in analogy with almost 

degenerate perturbation theory, the requirement of mathematical con- 
sistency gives only one equation to determine the two unknown first- 
order energies . 

21 

In the Hirschfelder-Silbey procedure, a further relation 
between (CoaO an< ^ is obtained by arbitrarily 

imposing the requirement that the energy through first order agree 
with the Heitler-London result, Eq. (3.6). This yields 

'4 - S* V' 


E (£°aD ~ 


/- $z 


4 


and 


(3.17) 


- (!) 


(txcM) = 


So© Vo< 


which clearly satisfy Eq. (3.16). 

Another way to fix the unknown constants is by a variational 
method analogous to that introduced in the almost degenerate perturba- 
tion theory (see page 36) . It was shown there that such a procedure 
can be expected to have a small effect upon the energy in a rapidly 
converging perturbation expansion. Thus, the question of applying 
the variational procedure to the present case cannot be answered with- 
out testing the convergence properties of the formalism with numerical 
examples. For the present, it is assumed that EE *** (Csui) and 
are fixed by Eq. (3.17). 
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With the constants in Eq. (3.15) fixed, the can be obtained. 

To this end, it is convenient to separate Eq. (3.15) inco a polarization 
equation 

5 <f> 0i +■ (///°- \J oe> ) (3.18) 

and an exchange equation 

( «/" - e’">) u <n = «•“> 

where 





The polarization equation is the first order equation in the un- 
symmetrical polarization expansion^ of ^^3 or . For the case 

of , it has been solved to high accuracy by variational methods 

by Hirschfelder and Lowdin^^ and by Kolos. ^ The primary effect of 



is to introduce correlations between the electrons which give 


rise to van der Waals-'interactions . The solution of Eq. (3.18) will 


not be considered further. 


The solution of the exchange equation will be considered in 
detail. In terms of 'the orbitals d 0 and i> 0 , Eq. (3,19) may be 
written 


/ $ a o) - e 0 &±C2.)-G 0 )co l * t} Oj& - 
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Thus, Eq. (3,20) is a coupled, two electron, two center equation. 
Introducing the complete set of eigenfunctions of and ^ , 

(k - e k )a k = o i - e A ) b k = o j 

the solution of Eq. (3,20) is 


2 ) ~ ^ ^ ^ ^ 

k?o JL$o e J. 


(3.21) 


*1. 


I £ 0 


"k 


4 « +21 a 0 (» 

.. jf i 

Jilt 


$ 


where zz &, — &> . 

k k 

The infinite sums in Eq. (3.21) may be evaluated by a technique 
45 

which has been used successfully for the calculation of atomic 
polarizabilities, van der Waals coefficients, etc. The first sun 

in Eq. (3.21) is analogous to dispersion energies in the polarization 
equation and the last two sums are analogous to induction energies. 

The sums for d>o and /£> O may be evaluated. using -the^identity , 







— 4 — If— 1 — + T-ld 

a'~~zUJ[A+lU j2,-T,U J 


(3.22) 


which can be shown to be an application of the residue theorem. 
Putting 0^=- 6^ j & £ substituting Eq. (3.22) into 
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Eq. (3.21) and interchanging the order of summation and integration 
yields 



and (3 . 24) 


( l £4 ) Q'O ^>0 ^*0 J> 

^ 4? = (£?)*, tj » ( k})^ r ■ The k* and J| 4 

3 5 

are the solutions of Eq. (3.24) for U = 0 In the application of 

the Hirschfelder-Silbey formalism to rig the first order exchange 

function tuf’® is simply ^ defined by Eq. (3.24). 
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Thus, it has been shown that can be obtained by solving the 

one electron , two center Eq. (3,24) and performing the quadrature, 

Eq. (3.23). 

Eq. (3.24) is still a coupled two-dimensional partial differential 
equation, and an exact solution in closed form does not appear possible. 

Variational approximations to the exact solution may be obtained by 

46 

finding the extreme points of the functional 

( 1 ^+^ 

. ^ . (3.25) 


j are trial functions. 

The insight. into the nature of is provided by the Unsold 

10 L Q* 

approximation to g. . The trial function is 




i > 0 - S 00 a„ 


z 4 


+ ZVL 


(3.26) 


where is a parameter determined from the extreme point of Eq. 
(3.25). A short calculation reveals 


i - 


L 


- I- 


i 

R 


+ 0(e *) 0 


(3.27) 


. Substitution of Eq. (3.27) into Eq. (3.23) and integration 


gives 
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eo,° } = i> 0 (oa 0 (z) -t- ^-[a 0 (o a 0 iz)-hi ki a){^J 


$0 


Z 


(3.28) 


^cCO lo^Cz) . 


The first and last terms in CO f ^ yield multiples of A{ ^ 
upon symmetrization, so that the primary effect of is to intro- 

duce ionic terms into the wave function. 

The approximation corresponding to Eq. (3.28) for the polarization 

.A 

function is 

4} 0} ~ '[* {*«/ 4z * y*l ibz ) + ^ Zarfb-JJd^CO 4 (3-29) 

where are the cartesian coordinates of the electrons. 

The 0^ and C l* are variational parameters which vary as JR ^ 

Eq. (3.29) results from evaluating by -.the Unsold method the lead term 

A . 

in the multipole expansion of . 


j (O 

With the approximations Eq. (3.28) and (3.29) for <f>. , 


wavefunction through first order is given by 


At (4>y -h £/'>) - £ / + ^ ^ 


the 

><*> 

(3.30) 


* Ai [*(*, V+ y» <j h *) « V*> 


+• 


Sea £ 0) Ce*<Ux) A / ^ , t \ 

— ttj Ai (a 0 toa 6 c*> + 4 co 4 / 
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Eq. (3.30) has precisely the form of the Hirschfelder-Linnett (HL) 
variational wavef unction. If ol , o^^and are fixed by the total 
variational principle, rather than a first order calculation such as 
Eq. (3.27), the wavefunction Eq. (3.30) yields the HL energy. With 
'fixed by Eq. (3.27), however, the coefficient of the ionic term 
•in the wavefunction is too small since £ (&XCk) varies as 

exp( — 3 R ). The HL calculations show that the correct coefficient 
( y in their notation) varies as roughly R (— R . ) 9 

Work is currently in progress to obtain $0^^ accurately. There 
is reason to believe that polarized ionic states 2 ^ occur in . 

In Chapter Four, a good variational approximation to the total first 
order wavefunction ^ ^ (denoted there by <7^ ^ ^ % ) is obtained, 
although the form of the basis set does not make it convenient to 
separate out ^ an ^ . The energies (labelled HS) reported 

in Chapter Four demonstrate that both good cpulomb and good exchange 
energies can be obtained' with the wavefunction A £ ( ^ ^ ^ )o 

3 . 2 Extension to Many-electron Systems . 

The method of solving the first order perturbation equation 
developed in the preceding paragraphs may be extended to the case of 
interactions involving many electrons. To avoid the complications 
which arise when the eigensolutions of the separated-atom system are 
unknown, consider the interaction of A/ groundj-state hydrogen atoms. 

A zeroth order primitive function is 


i 


(e) 


o . 0 (d K (*) • * • yi e cm) 
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This function is degenerate with N ! — / other functions 

» « * yi c 

where is one of the N{ permutations of the electron labels. 
The first order equation analogous to Eq. (3.14) is 


(H,“>-6 M )<k ln + n,"‘4,“> - Z ?4.‘" &,[ 


(O 


This may be solved by puttinj 



which yields the polarization equation 

+ ( ///"- (tf'lHf'Htftyf-i J/3-31) 

and the exchange equations 

(h, u> - e “M‘" = <4 a 'm;'‘>i a> . 0-32) 


In these equations 



and £ < * — A / & 0 . * 
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The solution to the polarization equation (3,31) has the form 


4 


o> 



±21 
a* (<i 


■+ 



±21 

aoox^it) 


where 


(H a l‘) €o ^~ ( y^2. ~ >2a .~~ V b! * '? ‘~^oo)^o C °k (i > 

which is precisely Eq. (3.18). Thus ±/° in the many electron case 
is obtained by solving a set of two electron equations. 

The exchange equation (3.32) is a coupled ^-electron equation if 
permutes the labels "of M> electrons. However, the identity for 

a O a 2.'” a k>°} _ 

(ah afXait 14) * ■ • (a% * ufc ) j 




* Co 

(l \fk , 

= ur) J du **' 


CO 

<ii V 


permits Che reduction of,Eq. (3.32) to & onej electron equations of 
the form of Eq. (3.24). This identity -'follows from Eq» (3.22). 

Thus, in the many ‘electron case the first order perturbation 
equation may be obtained by solving one and two electron equations. 

Summary . The extension to exchange perturbation problems of the 
formalism developed in Chapter Two leads to the Hirschfelder-Silbey 
perturbation theory if the undetermined elements of ^ are chosen to 
yield the Heitler-London first order energy. The first order wave- 
function may be solved as the sum of the first order polarization 
function, which introduces van der Waals correlations, and an exchange 
function, which introduces ionic terms. Xn the many electron case, the 
first order equation reduces to a set of one and two electron equarions. 
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IV. MODEL EXCHANGE PERTURBATION THEORY CALCULATIONS 

In this Chapter four different types of perturbation theories 
for exchange forces are applied to three model problems: 

1) the ground and first excited state qf the hydrogen molecule 
at internuclear separations R = 4, 6, 8 CLq \ 

2) a double minimum harmonic oscillator model of 

3) the double minimum delta-function model of 

Since much of this Chapter has already appeared in print, the 
published articles are reproduced here -with Addenda where appropriate. 
The delta-function calculations have not been published, previously. 
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2.1 Interaction of Two Hydrogen Atoms. 


Reprinted from The Journal of Chemical Physics, Vol. 49, No. 1. 24-34, 1 July 1968 
Printed in V. 5. A . * * * § 


Exchange and Coulomb Energy of H 2 Determined by Various Perturbation Methods* 

P. R. Certain! and J, 0. Hirschfelder 
Theoretical Chemistry Institute, University of Wisconsin, Madison , Wisconsin 

AND 

W. Kolos! and L. Wolniewicz§ 

Laboratory of Molecular Structure and Spectra , Department of Physics, University of Chicago, 

Chicago, Illinois 
(Received IS December 1967) 

Four different types of perturbation theories for the exchange forces between two atoms are applied to 
tbe ground and first excited state of the hydrogen molecule at internuclear separations R=4,6,8co. The 
energy through second order and the expectation value of the Hamiltonian using the wavefunction accurate 
through first order are calculated to compare the theories. The results for the Hirechfelder-Silbey procedure 
are satisfactory. The Murreli-Shaw or Musher-Amos results are equally good with the exception of the 
Hamiltonian expectation values for both states at i?=6 and Soo, which are bad. The Eisenschitz-London, 
van der Avoird, or Hirschfelder CHAV) results are good at small separations but at large separations they 
give a second-order energy which appears to be about one-half the correct dispersion energy. The Rayleigh- 
Schrodinger treatment using a Sternheimer type of zeroth-order Hamiltonian gave the best energy for the 
ground state but not very good energy for the excited state. At the separations considered, deviations from 
the *virial theorem are unimportant. 


I. INTRODUCTION 


In a previous paper, 1 the mathematical problems 
associated with the development of a satisfactory 
perturbation theory for the exchange forces between 
two atoms were discussed. Many formalisms have been 


*Tbis work was supported in part by the U.S. National 
Aeronautics and Space Administration Grant NsG-275-62, 
National Science Foundation Grant GP-7774, and by Advanced 
Research Projects Agency through the U.S. Army Research Office, 
Durham, N.C., under Contract No. DA-31-124, ARO-D-447, 
AR PA ORDER 368. 

t National Science Foundation Graduate Fellow. 

t Permanent address: Department of Theoretical Chemistry, 
University of Warsaw, ul. Pastcura 1, Warsaw 22, Poland, 

§ Permanent address: Department of Theoretical Physics, 
Nicholas Copernicus University, Torun, Poland. 

* J, O, Hirschfelder, Cbem. Pbys, Letters 1, 326, 363 (1967).. 


proposed during the past few years and most of them 
give different results for the second and higher order 
energies. From the mathematical standpoint, each of 
these treatments is equally satisfactory. Thus, we ask 
the question: “Which formalism provides the best 
second-order energy and the best expectation value of 
the Hamiltonian using the wavefunction accurate 
through the first order?” In the present paper, an 
attempt is made to answer this question on the basis 
of accurate numerical calculations. 

Wc consider both the ground state ( , S B + ) and the 
first excited slaLc ( 3 2 H 4 *) of the hydrogen molecule at 
the internuclear separations 72—4, 6, and 8«o- Ihc 
energy through the second order and the expectation 
value of the Hamiltonian corresponding to ^(1, X) = 
(where the constant X is either set equal to 
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one or else energy optimized) arc calculated using four 
different types of perturbation treatments (EL-HAV, 
HS, MS-MA, and RS). A rather large basis set was 
used together with vibrational principles to determine 
the first-order wavef unction The accuracy of the 

calculations was limited by the rounding errors which 
resulted from the single-precision computational pro- 
gram. 

The results, discussed in bee. IV, are very interesting, 
but they do not show that any of the four perturbation 
schemes is outstandingly better than the oilier three. 
The Rayleigh-Schrddmger (RS) is best for the ground 
state but is not very good for the first excited shite. 
The Hirsch felder-Silbey (IIS) results are all satis- 
factory. The Murrell -Shiuv (MS) or Mushcr-Amos 
(MA) energies are equally good, with the exception of 
the Hamiltonian expectation values for 12=6 and 8n 0 
where tire values aie exceptionally bad. At small 
separations the Eisenschttz-London, van der Avoird, 
or Hirschfelder HAV (EL-HA V) energies are com- 
paratively good, but at large separations the EL-HA V 
second-order energy becomes approximately one-half 
of the accurate dispersion energy. 

The energy of the triplet state, added to or subtracted 
from the energy of the singlet state, yields the Coulomb 
or exchange energy, respectively. A surprise to us is 
that the perturbation procedures give better values 
of the exchange energy than of the Coulomb energy. 
To test whether this result is a property of the per- 
turbation expansions or is due to our use of inaccurate 
approximations to the first-order wavefimetions, we 
have also computed total energies using the Rayleigh- 
Ritz variational principle and the same basis sets as' 
were used in the perturbation calculations. These com- 
putations show that the perturbation expansions of 
the Coulomb energy are slowly convergent. 

In terms of computational efficiency, none bf the 
perturbation schemes offers any advantage over a 
standard Rayleigh-Ritz variational calculation of the 
total energy, especially if it is necessary to use the 
trial wavefunction f to obtain a reliable estimate of 
the energy. The ffiajor portion of computing time is 
spent in calculating* matrix elements. Since the same 
elements appear in both methods, a more accurate 
energy can be obtained by varying the total wave- 
function than by varying the first-order wavefunction. 
The matrix elements for the EL-HAV, HS, and MS- 
MA second-order energy calculation are significantly 
simpler, however, since the integrands are not sym- 
metrized. Perturbation techniques are useful in diag- 
nosing the defects in approximate wavefimetions, in 
calculating accurate expectation values of properties 
other than energy, and in determining upper and lower 
bounds. Their practical value in determining the energy 
depends upon our ability to determine sufficiently 
accurate solutions to the perturbation equations with- 
out recourse to the use of large basis sets. 


II. PERTURBATION FORMALISMS USED 
m THIS PAPER 

The oldest perturbation theory which yields an ex- 
pansion of the exact wavefunction in terms of sym- 
metrised products of atomic orbitals is that of Eisen- 
schitz and London (EL) * More modern theories which 
have the same first-order wavefunction and second- 
order energy have been developed by van dcr Avoird 8 
and by Hirschfelder 1 (HAY). Wc also consider the 
recent perturbation formalisms of Hirschfelder and 
Silbey (HS), 4 Murrell and Shaw (MS), 6 and Mosher 
and Amos (MA), 5 as well as the usual Rayleigh- 
Schrddmgcr (RS) procedure using a Stcrnheimer-type 
zeroth-order Hamiltonian. 7 -* The MS and MA formal- 
isms arc equivalent 9 so that we are dealing with four 
independent perturbation procedures. 

The general structures of the various theories have 
been fully treated elsewhere.*-* Here we simply Intro- 
duce a notation that Is suitable for treating the lowest 
*2 + and S 2 U + statesof H 2 , The exact energy and wave- 
function are denoted by *TS and respectively, with 
i~g for the gerade state and i~u for the ungerade; 
state. The functions and Y are eigenfunctions of 
the projectors 

(la) 

U*>) 

where P 13 permutes the electronic coordinates and I 
inverts the wavefunction through the midpoint of the 
intemuclear axis. 

The unperturbed Hamiltonian is defined to be 

M*= - ( 2 ) 

with eigensolutions («*, $&)» The functions are thus 
products of hydrogen atom wavefimetions. The per- 
turbation V Is 


* v-JtH-nrt-nr*--*?*. ( 3 ) 

The exact energy and wavefunction are resolved 

into perturbation series: 

» ( 4 ) 

. ( 5 ) 


1 R, Eisenschitz and F. London, Z. Fhysik 60, 491 (1930); 
A. van dcr Avoird, Chem. Phvs. Letters 1, 24 (1967). 

* A. van dcr Avoird, j. Chem. Phys. 47, 3649 (1967); Chem. 
Phys. Letters 1, 411 (1967), 

*J. O. Hirschfelder and K. Silbey, J. Chem. Phys. 45, 21S$ 
(1966). 

* J. n. Mm roll and O.’Shaw* J. Chem. Plus, 46, 176S (1967). 

* j. I. Masher and A. T. Amm (Preprint) ; A. T. Amos and J. I» 
M usher, Chctn. Phys, Letters 1, 1 19 (1967). 

7 R. hi. Stcmhcimer, Phys. Rev, 96, 951 (1954). 

* J, 0. Hirschfelder, W. Byers Brown, and S. T. Epstein, Advan. 
Quant. Chem. 1, 255 (1964).. 

f R. E. Johnson and S, T. Epstein, University of Wisconsin 
Theoretical Chemistry” Institute Report WISTCX-265, 3 
November 1967. 
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Table I. Definition of the operator Ti appearing in Eq. (7). 


Formalism 

r. 

EL-HAV 

A t {V-*EV>) 

HS 

AAV-oE^+AAV-^EM) 

MS-MA 

V 


Since our present calculations arc limited to 'E through 
the second order and expectation values of the Hamil- 
tonian for V'o+V (i) j we do not need to determine the 
wavefunction beyond the first order. 

The zeroth- and first-order energies are given by the 
same expressions in all theories considered: 

eo“ — 1; •£«>— (0 O | VAtfo)/ (0 O 1 vl,0 o ). (6) 

It is in expressions for the second-order energy that 
we encounter differences. The most direct way to 
compare the EL-HAV, HS, and MS-MA results is 
to employ the spectral expansion of i E (2 \ which may 
be written 

‘£«>= (<h | A,i, 0 )-‘Z 1 **><& i im, 

kpH) CQ—ei k 

( 7 ) 

Here T t is an operator which depends on the formalism 
considered. The explicit expressions for 2\- are given 
in Table I. 

Because of the form of *£ (2) , we find it advantageous 
to define a function 

<X- E [<^(K- '£«>)*, | &>*/(*-«)>*. (8) 

h/0 

Note that *x depends upon the symmetry of the state 
considered through the projector A i} although f x itself 
has no simple symmetry properties. The second-order 
energy in each formalism may be obtained by a single 
* integration involving *x. That is, 

<*xir,0o)/<^M^>. (9) 

The function f x is related in a simple way to the first- 
order wavefunction in the various formalisms. 

For the Eisenschitz-London or HAV scheme, 1 *- 3 

V (l) (EL-HAV) = A ( (0o l A fx )/ (0o [ Ai4> 0 )) A <0 O , 

( 10 ) 

The second term on the right in (10) assures the or- 
thogonality of V (1) (EL-HAV) to A ,0o, which is the 
normalization condition of van der Avoird 3 and Hirsch- 
felder. 1 

In the Hirschfelder-Silbey expansion, 4 is ob- 
tained by projecting a function of the proper symmetry 
from the sum ff x+ u x : 

(HS) = A i( e x+“x) • (11) * 

Finally, the Murrell T Shaw 5 and Musher-Amos 6 
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first-order wavefunction is { x itself, 

V (1) (MS-MA)- 4 X . (12) 

Note that this formalism has the unique feature that 
the first-order wavefunction does not possess proper 
symmetry. 

The contribution of continuum states makes the 
sum (8) difficult to evaluate. We shall therefore resort 
to a variational principle for *x* It follows from (8) 
that *x satisfies the differential equation 

(tfo-€o) *X+M V- *'£<«) 0 O = 0. (13) 

Since we are considering the states which arise from 
the ground state of H 0 , the Eq. (13) is equivalent to 
the variational principle 5*/= 0, where 

< j- <‘X | (# 0 — €o) *x >+ (fx | A i( V - \E (1) ) 00 ) 

4- (0o | ( V— *'jE {1 >) A f‘x ), (14) 

and *x is a trial solution of (13). Note that 'J does 
not depend on the component of 0o in *x. We use this 
freedom to satisfy the normalization condition 

(*xl0o)~O. 

We choose the trial function to have the form 

S-E'Cjx/, 05) 

i 

XJ= exp[-f 22]($i+fc+ ( 16) 

where £ and tj are the usual elliptic coordinates (£= 

( r a -\-rb)/R and ??= (r a ~~r b )/R) ; p— 2rn/R, where rn 
denotes the interelectronic distance; r , s, f, s are in- 
tegers in the range 0-4; and n—Q, 1, 2. The linear 
coefficients ‘Cy are variational parameters, and we 
have taken N=3 0, 50, 70 for i?=8, 6, 4co, respectively. 

For one- and two-electron systems where the wave- 
function is a product of a space function times a spin 
function, it is possible to use standard Raylcigh- 
Schrodinger perturbation theory in which a Stern- 
heimer type of zeroth-order Hamiltonian H oi , corre- 
sponding to the zeroth-order wavefunction At, 0 O , is 
employed. Here 

H 0 iA{(f)o—€oAif{) 0 and [77o,-, Af\—0. (17) 

The perturbation is then 

(18) 

The Sternheimer Hamiltonian is the only Hamil- 
tonian which satisfies Eq. (17) and can be written as 
the sum i7oi=A-f-£7o,- where K is the kinetic energy 
operator and £ 7 0 » is a local potential energy function. 
Jansen 10 and Corinaldesi 11 have developed Rayleigh- 
Schrodinger-type treatments which might be applied 
to many-electron systems but their zeroth-order Hamil- 
tonians are non-Hermitian and involve nonlocal poten- 

10 L. Jansen, Phys. Rev. 162, 63 (1967). 
n E. Corinaldesi, Nuovo Cimento 25, 1190 (1962): 30, 105 
(1963); E. Corinaldesi and H. E. Lin, ibid. 28, 105 (1963),- 
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tials. The Sternheimer potential energy function is 
?7 o.=€o— (KA^/iA^o). If now we remember that 
Ho<f>Q= eo4>o, that Bq~K-]-Uq, and that K commutes 
with A { , it follows that 

H„^K+(A t UM/{A4 o). (19) 

Then using our previous notation, V—U— U 0 , and 
remembering that U commutes with A { , it follows that 

Vi= U- {A&M/{A rf o) = (. AiV<h)/{A#t>o ) . (20) 

Following the usual Rayleigh-Schrodinger pro- 
cedure, the equation forty 10 (RS) is 

(^o-e 0 )*^»)(RS)+(F i ~*£<»)i4^o-0, ‘ (21) ' 

where is given by Eq. (6). Because of Eq. (20), 
we may also write Eq. (21) in the form 

<(tfo-*o)-HD4 f , Uo]WA x ^m {1) ( RS) 

( 22 ) 

Equation (22) may be compared with Eq. (13). Again 
making use of Eq; (20), the Rayleigh-Schrodinger^ 
second-order energy is 

*‘£< 2) (RS) - i V 0) (RS) )/(4o I A#*). 

(23) 

Of course we cannot solve Eq. (21) exactly, but we 
can use the Hylleraas variational principle 12 - 8 to deter- 
mine an upper bound for the second-order, energy 
*£ (2) (RS) and an approximation to the first-order . 
wavefunction ty<° (RS). .By .virtue of Eq. .(20), the 
Hylleraas principle can be written 

*0»(RS) fa, I A ^ .) 

= (^<‘>(RS) I (ff w -*)<?o>(RS)) 

+ <^ 0) (RS) | 


of basis functions used was N'= 50 ,‘ : 35, and 20 for R= 
4, 6, and 8ao, respectively. 

Thus, for all four types of perturbation procedures 
which we discuss, , the first-order wavefunctions are 
expressed as linear' combinations of the xy» The sub- 
stitution of Eq. (15) into Eq. (14), or Eq. (25) into 
Eq. (24) , together with the variation of the *‘Cy leads 
in the usual way to a set of simultaneous inhomogeneous 
linear equations. The solution of these equations is 
hampered by the near linear dependencies in the set 
(16) when the expansions (15) or (25) contain a large 
number of terms. Since the calculations were done in 
single precision, round-off errors became significant 
before convergence to the true solutions of (13) and 
(21) was obtained. At large values of R, accurate values 
of the integrals were difficult to calculate. Neverthe- 
less, it is believed that the second-order energies ob- 
tained possess at least two significant figures (which is 
sufficiently accurate for our purposes). The calcula- 
tions were made on the University of Chicago- Com- 
puting Center IBM 7094 Computer. 

Having obtained the approximate first-order wave- 
functions ty (1) , it is useful to Construct a trial function 
for the total Hamiltonian Hq+V of the form 

<?(1 = (26a) 

where X is a variational ‘parameter. The expectatior 
value of the'Hamiltonian corresponding to ty, 

f J(l,X). 

- mi, X) I (1, X) >/<¥(!, X) 1-^(1, X) ), (26b) 

should' then give the energy accurate 'through the 
rd order of the perturbation and give an -tippet 
und to the true energy of the system. 

HI. PREVIOUS PERTURBATION CALCULATIONS 
FOR H 2 


+ (A { (V-*E<»)<h I $»(RS) ). (24) . 

By varying $< 1 >(RS) to make \0< 2 >(RS) a minimum 
we obtain the best approximation to ty (0 (RS) and the 
best approximation. \£< 2 >(RS). Since* H 0{ and K,- com- 
mute with A,, we can choose the -trial function to have 
the form 

(25) 

where the xy are defined by Eq. (16) and the Sy+sy are 
restricted to be even integers. Because out trial func- 
tion varies linearly with the variational parameters 
'Cy, it is easy to show that for the optimum values of 
l Cj the *i5 <2) (RS) of Eq. (24) is equal to the second- 
order energy of Eq. (23) if the exact ty«>(RS) is re- 
placed by ty (1) (RS). In our calculations the number 

. A ; Hylleraas, Z. Physik 48, 469 (1928) ;• ibid. 65, 209 
(1930) ; also, “The Variational Principle in Quantum Mechanics '* 
Kept. No.' 1, Inst. Tbeoret. Phys., University, of Oslo, 1961.- * - 


In previous papers by Liu, Lyon, an< _ _ 

Goodisman 14 ; and Matcha ’and Byeis 'Brown, 15 the 
wavefunction and energy for the ground state of the 
hydrogen molecule have been "calculated by pertur- 
bation procedures. However, in this previous work the 
zeroth-order wavefunction was not taken to be the 
symmetrized product of the. separated atom wave- 
functions and therefore the mathematical problems 
were quite different from those which we encounter. 
Liu, Lyon, and Byers Brown 13 and Goodisman 14 take 
^o=exp[— c(^i4-^)] and consider, the range from R— 
1.35ff 0 -L45ffo. Matcha .and Byers ' Brown 15 lake the 
zeroth-order wavefunction to be the product of the 
ground-stale wavefunction for II 2 + for each electron 
and consider the range from R= 0-2.2ff 0 . With a five- 

13 B. Liu, W. D. Lyon, and W. Bvers Brown, J. Chcm. Phys. 44, 
562 (1966). 

11 J. Goodisman, J. Chcm. Phys. 47, 1256 (1967). 

15 R. L, Matcha and W. Byers Brown, J. Chem, Phys. 48, 74 
(1968)-, ■ ■ - ' 
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Table II. Dissociation results of Matcha and Byers Brown 
(Ref. 15). 


n Errors in D—D» 


1 - +2.7509 eV ' 22189 cm"- 1 

2 -0.2683 ”2164 

3 +0.0293 * 236 ' 

4 +0.0041 33 

5 +0.0011 9 


term basis set (with optimum scaling) Liu, Lyon, and . 
Byers Brown 13 determined the energy through second 
order at the equilibrium separation which corresponded 
to an energy of dissociation equal to 4.98S eV=40233 
cm" 1 or 105% of the experimental value; whereas this 
same basis set gives a Rayleigh-Ritz variational energy 
corresponding to 95% of the experimental value. 
Goodisman 11 used a 10-term function to calculate the 
energy through the third order to obtain a dissociation 
energy equal to 4.617 eV=37 241 cm -1 or 97.2% of the 
experimental. Matcha and Byers Brown 15 used a 50- 
term basis set and calculated the energy at R=*1.4a 0 
through the fifth order. The accurate Kolos and Wolnie- 
wicz 16 variational calculation for the dissociation energy 
D= £( oo) - £= 4.7474 cV= 38 293 cm" 1 . If D n is equal ; 
to £{ oo) minus the Rayleigh-Schrodinger perturbation 
energy calculated through the «th order, then Matcha 
and Byers Brown 16 obtained the results shown in Table 
II. Thus, it appears that the Rayleigh-Schrodinger per- 
turbation sequence converges very rapidly under con- 
ditions where the zeroth-order wavefunction has the 
same symmetry as the perturbed wavefunction. 

IV. DISCUSSION OF RESULTS 

The results of our calculations are given in Tables 
III-XHI. Let us examine separately each' of the tables. 

The first two tables show how the second-order 
energies of the ground state and the first excited state 
vary with the size of the basis set. In the remaining- 
tables, the numbers quoted correspond to the largest 
values of N listed. 

Tables V and VI compare the accurately calculated 
variational energies, 19 i £, with f £®=€o+ , ’£ a) + i £ (2) 
and i E(l, X) as determined from the various pertur- 
bation procedures. We note the following. (We classify 
an energy as “good” if it is within r^lO% of the accurate 
energy.) 

Ground-state energy through second order, RS is very 
good, and gives the best energy at 4 o 0 . HS and MS-MA 
arc both good and give virtually the same values. These 
schemes give the best results at 6 and KL-1 f AV is 
fair at 4oo, but is bad at 8f/ 0 . Examination of Table 
III shows that this is due to °E i2) ( EL-HA V) being 
equal to approximately one-half of the second-order 
energy of the other schemes at 8oa- This is discussed 
in detail below. 

19 W. Kolos and L .Wclmewicz J. Cfaem; Pfays. 41, 3663 (1964) . 
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J Expectation ■value of the Hamiltonian for the ground 
slate. RS is good at 4 and 6co, but poor at 8«o. The 
energy is almost unchanged when X is varied since the 
optimum X is nearly unity. HS is quite good at all 
separations. The optimum X is nearly unity. MS-MA 
is almost as good as HS at 4«o, but is bad at 6 and 8ao. 
Again, X is nearly unity. EL-HAV with X=*l gives 
values between those of HS and MS-MA. Variations 
of X give considerable improvement in the energy, 
especially at 8a 0 where the optimum X is almost 2. This 
is related to the behavior of the second-order energy. 

Excited stale energy through second order . RS is not 
good at any separation. This is a result of the pertur- 
bation F« of Eq. (20) being very large due to the node 
of A u p 0 , n HS and MS-MA again give virtually the 
same results, which range from bad at 4 Oq to very good 
at 8c fi . EL-HAV is bad at all separations. We again 
note the factor of § in the second-order energy at 8co- 

Expectation value of the Hamiltonian for the excited 
stale . RS is poor at all separations and not improved by 
variations in X. HS is fair at all separations. The 
optimum X is nearly unity, MS-MA ranges from fair 
at 4a 0 to bad at large separations. Variations of X fail 
to improve the energy. EL-HAV exhibits behavior 
parallel to that in the ground state. 

Tables VII and VIII compare the accurate 19 * 18 
Coulomb energy, £coui= ■!(*£+“£), and the exchange 
energy, EUeb =£(*£-- U E), with the results of the 
various perturbation schemes. The exchange energy 
is very good in all of the methods, especially the 
Rayleigh-Schrodinger. Surprisingly, the Coulomb en- 
ergy (which is essentially dispersion energy) is much 
less accurate for all of the treatments. 

It is possible that the poor results for Coulomb energy 
are a consequence of the limited size of the basis set 
used to compute the second-order energies, rather 
than of slowly convergent perturbation expansions. 
Thus, we have computed the total energy of the ground 
and excited states using the Rayleigh-Ritz variational 
principle and the same basis sets as were used in the 
perturbation calculations. The results are given in 
Table IX and show that very accurate values of the 
Coulomb and exchange energies may be obtained if 
total energies are computed with the basis sets listed 
in Tables X and XI. Thus we believe that our com- 
puted second-order energies are similarly accurate. 
Note that using the same basis set as in the perturba- 
tion calculations, at R — 8f? 0 we calculated the Rayleigh- 
Ritz variational energy = 11.74 cm" 1 which is 

0.15 cm" 1 belter than the best previous energy ns de- 
termined by Kolos and Wolniewiez. 10 

Tables X and XI give the coefficients *'C/ for the 
wavef unctions 'x and ^^(RS), respectively. 


17 We wish to thank Professor S. T. Epstein for stressing this 
point to us. 

18 J. O. Hirschfeldcr and W, J. Meath, Advan. Chem. Phys. 12, 

66 (1967). • ** * 
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Table III. ~°E m (in atomic units) as a function of N, the number of terms in the basis set 


N 

R 

_£a 

10 

20 

30 

40 

50 

60 

70 

4 EL-HAV 
HS 

MS-MA 

RS 

0.0014217 

0.0018221 

0.0018903 

0.0026307 

0.0024212 

0.0028701 

'0.0029819 

0.0034559 

0.0027206 

0.0030362 

0.0031994 

0.0036736 

0.0028814 

0.0030805 

0.0032347 

0.0037517 

0.0029628 

0.0030831 

0.0032586 

0.0037975 

0.0030960 

0.0031350 

0.0033050 

0.0031377 

0.0031898 

0.0033533 

6 EL-HAV 
HS 

MS-MA 

RS 

0.0000784 

0.0001459 

0.0001459 

0.0001811 

0.0001459 

0.0002465 

0.0002466 

0.0002510 

0.0001577 

0.0002585 

0.0002589 

0.0002782 

0.0001700 
0.0002673 
0.000267 6 
0.0002784“ 

0.0001897 

0.0002797 

0.0002803 



8 EL-HAV 
HS 

MS-MA 

RS 

0.0000062 

0.0000121 

0.0000121 

0.0000272 

0.0000159 

0.0000312 

0.0000312 

0.0000334 

0.0000181 

0.0000350 

0.0000350 





* Obtained with 35 terms. 








Table TV. —“£<*> (in atomic units) as a function of iV, the number of terms in 

the basis set. 


N 

R 

10 

20 

30 

40 

50 

60 

70 

4 EL-HAV 
HS 

MS-MA 

RS 

0.0005320“ 

0.0009621 

0.0009553 

0.0010077 

0.0007929 

0.0012752 

0.0012682 

0.0011602 

0.0009977 

0.0013367 

0.0013161 

0.0012020 

0.0011131 
. 0.0013271 
0.0013036 
0.0012091 

0.0012329 

0.0013622 

0.0013113 

0.0012483 

0.0013658 

0.0014077 

0.0013467 

0.0013952 

0.0014511 

0.0013775 

6 EL-HAV 
HS 

MS-MA- 

RS 

0.0000650 

0.0001328 

0.0001328 

0.0001432 

0.0000S67 

0.0001877 

0.0001878 

0.0001684 

. 0.0000907 
0.0001920 
0.0001922 
0.0001753 

0.0000927 

0.0001905 

0.0001906 

0.0001785* 

0.0000973 

0.0001877 

0.0001878 



EL-HAV 

HS 

MS-MA 

RS 

0.0000057 

0.0000116 

0.0000116 

0.0000263 

0,0000148 

0.0000301 

0.0000300 

0.0000312 

0.0000161 

0.0000331 

0.0000330 






• Obtained with 35 terms. 


Table V. Comparison of perturbation energies with accurate variational energies, (in cm -1 ). %+ ground state, H 3 .» 


Exact 


R 


eo-p£ (,) -*E 

EL-HAV 

HS 

MS-MA 

RS 






~°E 


4 

3592.6 

1115.7 

427.0 

415.6 

379.7 

282.2 

.6 

178.9 

67.1 

25.46 

5.70 

5.58 

5.99 

8 

11.6 

7.8 

3.78 

0.06 

0.06 

0.37 





»E( 1, l)—°E 


4 

5592.6 - 

1115.7 

243.8 

301.8 

338.5 

168.4 

6 

178.9 

67.1 

13.93 

7.56 

27.13 

7.51 

8 

11.6 

7.8 

1.87 

0.14 

3.88 

1.36 





*E(l,X)-*£ 


4 

1592.6 

1115.7 

177.8 

278.1 

323.5 

151.1 




(1.36) 

(1.20) 

(1.16) 

(1.16) 

6 

178.9 

67.1 

9.41 

7.51 

27.07 

7.48 




(1.39) 

(0.97) 

(0.96) 

(0.98) 

8 

11.6 

7.9 

0.12 

0.14 

3. 88 

1.29 




(1.92) 

(0.99) 

• (0.97) 

(0.90) 


* Here S E is the accurate variationaliy calculated energy (Ref. 16); wavefunction through the first order taking X«1 [see Eq. (26)J; and 

1b the energy calculated through the second order of perturbation; °E{ 1, X) is the Bame as °E(l, I) except that X is optimized to make 

*£(1. 1) is the expectation values of the Hamiltonian calculated with the «E(1. X) a minimum. The values of X are given in parentheses. 
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Table VI. Comparison of perturbation energies with accurate variational energies, (in cm" 1 ) . excited state, Hg. a 




f0+UjE (D- U £ EI^-IIAV 


HS 


->1453.3 

-41.2 

4.3 


475.8 

44.8 

7.3 


169.6 

23.47 

3.77 


«S «>— «jE 


157.3 

3.62 

0.05 


173.4 

3.61 

0.06 


201. b 
5.66 
0.42 


4 

-1453.3 

475.8 

6 

-41.2 

44.8 

8 

4.3 

7.3 


4 

-1453.3 

475.8 

6 

-41.2 

44.8 

S 

4.3 

7.3 


* Here “£ is the accurate variationally calculated energy (Ref. 16); 
j s the energy calculated through the second order of perturbation; 
«E(|, 1) is the expectation value of the Hamiltonian calculated with the 

Table XII gives the expectation values of the kinetic 
energy K for the EL-HAV and HS schemes using the 
trial wavefunction *$(1, X) of Eq. (26a) with the 
optimum value of X. This permits the determination 
of a scaling constant s for the approximate wavefunction 
so that the virial theorem is satisfied. However, as can 
be seen from Table XIII, even at 4c<>, the values of s 
are so near to unity that the improvement in *25(1, X) 
is very small. At larger separations, the improvements 
is negligible. The algebra involved in this scaling process 
is discussed in the Appendix. 


“£(!,!)-"£ 


78.1 

131.3 

145.7 

219.7 

13.40 

4.35 

23.69 

* 6.87 

2.01 

0.31 

3.87 

1.01 


«E(1,X) 



39.6 

129.4 

143.6 

218.7 

(1.42) 

(1.0S) 

(1.09) 

(0.94) 

4.41 

4.33 

23.69 

.6.85 

(1.89) 

(0.98) 

(0.99) 

(0.97) 

0.29 

0.30 

3.86 

0.99 

(1.98) 

(0.96) 

(0.97) 

(0.94) 


wavefunction through the first order taking X = 1 [see Eq. (26)); and 
”22(1, X) 13 the same as “22(1, 1) except that X is optimized to make 
“22(1, X) a minimum. The values of X are given in parentheses. 


Finally, we consider the behavior of i E (Si (EL-HAV) 
noted above. From Tables III and IV, we see that at 


» -«£< 2 ) 

EL-HAV 3.97 cm- 1 3.53 cm" 1 

HS 7.68 7.26 

MS-MA 7.68 7.24 

RS 7.33 6.85. 


At i?=8co, the accurate value of the second-order 
dispersion energy 10 * 18 (not considering the exchange of 


Table VII. Comparison of perturbation Coulomb energies with accurate variational Coulomb energies (in cm -1 ).* 


Exact 


R 

«0— Ecoul 

< 0 -f£coui°>-£coui EL-HAV 

HS 

MS-MA 

RS 





SCou» (S>— 

•Seoul 


4 

1069.6 

796.1 

298.3 

286.4 

276.6 

242.0 

6 

68.8 

55.9 

24.47 

4.66 

4.60 

5.83 

8 

8.0 

7.6 

3.78 

0.06 

0.06 

0.40 





■Eooultt, 1) 

— -Ecoui 


4 

1069.6 

796.1 

160.9 

216.6 

242.1 

194.0 

6 

68.8 

55.9 

13.67 

5.96 

25.41 

7.19 

8 

8.0 

7.6 

1.94 

0.22 

3.88 

1.19 





Should, X) 

— Ecool 


4 

1069.6 

796.1 

108.7 

203.7 

233.5 

184.9 

6 

68.8 

55.9 

6.91 

5.92 

25.38 

7.17 

8 

8.0 

7.6 

0.20 

0.22 

3.87 

1.14 


* Here Ecoul*=J (*£+"£); the Coulombic energy through the second the first order is EcoulO. X). The values given are half the sum of the 
order of perturbation is £coui W ; and the Coulombic energy obtained from corresponding values in Tables V and VI. 
the expectation values of the Hamiltonian using the wavefunction through 
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Table VIII. Comparison of perturbation exchange energies with accurate variational exchange energies (in cn* 1 ) 
Exact 

R -E,„ b &, xch W-E, tch EL-HAV HS ' MS-MA RS 


4 

2423.0 

320.0 

6 

110.1 

11.2 

8 

3.6 

0.3 


4 

2423.0 

320.0 

6 

110.1 

11.2 

8 

3.6 

0.3 


4 

2423.0 

320.0 

6 

110.1 

11.2 

8 

3.6 

0.3 


“Here = “£); the exchange energy through the second 

order of perturbation is S #I0 h W ; and the exchange energy obtained from 
the expectation values of the Hamiltonian using the wavefunction through 

electrons between the atoms) is —8,0 cm -1 . Thus, the 
values of 'E {2) calculated by the HS, MS-MA, and RS 
are approximately equal to the dispersion energy 
£<*>( dis). However the (EL-HAV) is only about 
half as large. This raises the question as to whether 
the *'.E (2) (EL-HAV) can be expected to approach the 
correct £ (2) (dis) values in the limit of large separations. 

Rigorously, the relation between. t 'E (2) (EL-HAV) and 
E (2) (dis) for large values of R is not a simple one, as the 
following arguments demonstrate. We first substitute 
(1) into (7) to obtain 

( EL-HAV) ]P i? 

frO €0 —£k 

+4E L (27) 

fro Co~ <Jt 

where the symbol ~ indicates asymptotic equality. The 
first summation is ££ (2) (dis). Although each term in 
the second summation is exponentially decreasing 
in R , the series may not be. 17 For example, 

£er*R»/»l=l. 

n«0 

A. Arguments That i £< 2) (EL-HAV)~A< 2 >(dis) 

In their oiiginal paper, 1 Eiscnschitz and London 
used the Unsold approximation to evaluate , 7£ (2) (EL- 
HAV). Taking average energies Ae 1 and Ae* from the 
first (direct) and the second (exchange) series, respec- 
tively, in (27) and using the relation P= 1 and the 
completeness of the states <£*, we obtain 

*£ (2) (EL-HAV) HAe^-FAe*- 1 ) | F 2 ^). ( 28 ) 


128.7 - 

129.1 

103.1 

' 40.2 

1.00 

1.04 

0.99 

0.17 

0.01 

0.01 

0.00 

-0.03 


•&„h(U) 

—R»«h 


82.9 

85.3 

96.4 

-25 6 

0.27 

1.60 

1.72 

0.32 

-0.07 ‘ 

-0.09 

0.01 

0.18 


-Eoxshflj A) • 

“’-EnMh 


69.1 

74.4 

90.0 

-33.8 

2.50 

1.59 

1.69 

0.32 

-0.09 

-0.08 

0.01 

0.15 


the first order is £e X eh(l> X). The values given arc half of the diiTerctw 
between the corresponding values given in Tables V and VI. 


Eisenschitz and London chose (Ae^+Ae*- 1 ) so thal 
the coefficient of R~ G in (28) equals the coefficient ol 
R~* in 2s (2) (dis). In this way they obtained the result, 
^(Ae t “ l +A6*~ 1 ) = (0.925)“ 1 . This result appears plaus- 
ible since the average excitation energy must be greater 
than 0.75, the energy required to excite the two hydro- 
gen atoms from the H to the 2s or 2 p states. In this 
xvay it would appear that at sufficiently large values 
of R , so that only the Rr* term makes appreciable con- 
tributions to the second-order energy, the \E (2) (EL- 
HAV) becomes approximately equal to the £ f2) (dis). 

Van der Avoird 19 has calculated for H 2 + the values 
of *E (2) and \E(1, 1) using his method (which is 
equivalent to the EL-HAV scheme) for R< 7o 0 and 
obtains excellent agreement with the exact values. 
However, for at R= 7 the exchange energy is still 
large so that he obtains 

*£< 2 >= -252 cm- 1 , “E< 2 >= -79 cnr 1 . 

For this separation Dalgarno and Lynn 20 obtained the 


Table IX. Comparison of the Rayleigh-RItz variational energies 
(in cm -1 ) .“ 


R 

*Et-pE 


f^Coul 

— -Seoul 


4 

14.1 

9.6 

11.8 

2.2 

6 

1.3 

2.7 

2.0 

-0.7 

8 

-0,1 

0.0 

0.0 

0.0 


* Here °E, “E, Econl. arc the accurate variational energies (Ref. 
16); e E', U E', E^Coul. -E »xch are the variational energies obtained with 
the same basis sets as were used in the EL-HAV. HS. and MS-MA per- 
turbation calculations. 

18 A. van der Avoird, Chem. Phys. Letters 1, 429 (1967). 
(1957) ^ garno an< * ^ ynn > ? roc - Phys. Soc. (London)' 70, '223 
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Table X. The coefficients { Cj for the wavefanction *x- 


Sc «c 


n 

r 

s 

f 

s 

R= 4 

R-6 

R=8 

11 

R=6 

R= 8 

0 

O' 

0 

0 

*0 ■ 

—0.051793 

-0,013709 

-0.071657 

0.15228 

0.10664 

-0.063921 

1 

0 

0 

• 0 

0 

0.45177 

0.43402 

0.18830 

0.0069818 

0.12027 

0.15467 

0 

1 

0 

1 

0 

-0.11001 

-0.031741 . 

-0.43256 

0.099046 

-0.018047 

-0.037959 

1 

0 

1 

0 

1 

-0.35631 

-0.67385 


0.60628 

0.20964 


0 

0 

1 

1 

0 

-0.16695 

-0.23975 


0.26337 

0.11657 


0 

0 

0 

0 

1 

0.091648' 

0.11947 


-0.34809 

-0.21667 


0 

0 

0 

1 

0 

-0.11157 * 

-0.24252 

-0.032870 

-0.12545 

-0.068830 

-0.024074 

0 

1 

0 

'0 

0 

-0.084751 

-0.16727 


-0.11369 

-0 065128 


0 

0 

1 

0 

1 

1.0504 

1.1978 

-0.19732 

-0.95234 

-0.64999 

-0.17513 

0 

0 

0 

0. 

2 

0.02836* 

-0.04032 

0.057268 

0.31171 

0.082430 

0.040280 

0 

1 

1 

0 

0 

-0.33335 

-0.13367 

-0.093795 

-0.33512 

-0.25195 

-0.066610 

0 

0 

2 

0: 

1 

0.60210 

0.62991 

-0.18095 

-0.55453 

-0.22399 

-0.14659 

0 

0 

1 

0 

2 

-0.57767 

-0.60717 

0.23027 

0.55846 ■ 

0.17938 

0.17494 

0 

0 

0 

1 

1 

0.14224 

0.38298 

0.055254 

0.35575 

0.17420 

0.042859 

1 

0 

1 

0 

0 

0.59150 

0.63977 

0.18317 

-0.25092 

0.037367 

0.14738 

1 

0 

0 

0 

1 

-0.63361 

-0.62502 

-0.17661 

0.23356 

-0.020123 

-0.14436 

0 

0 

1 

0 

0 

-0.23905 

-0.28296 


0.33933 

0.31067 


0 

0 

2 

0 

0 

-0.0S539: 

-0.08783: 


0.29472 

0.15259 


0 

0 

2 

1 

0 

-0.14S06 

-0.10394 

0.028112 

0.30545 

0.12946 

0.17588 

0 

1 

0 

0 

2 

-0.18076 

-0.039101 

0.014502 

0.29755 

0.14846 

0.0072865 

0 

0 

.1 

1 ■ 

1 

-0.14124 

0.03938- 


0.075088 

0.10264 


0 

0 

3 

0 

0 

0.052181 

0.14873 

0.0059048 

0.25660 

0.010843 

0.00040308 

0 

0 

0 

0 

3 

-0.082322 

-0.11987 


-0.26058 

-0.017253 


0 

1 

2 

0 

0 

0.10710 

- 0.03649< 


-0.24275 

-0.11048 


0 

0 

0 

1 

2 

-0.014495 

-0.11473 


-0.25695 

-0.069458 


I 

0 

2 

0 

0 

0.24485 : 

0.39467 

0.061095 

-0.19842 

-0.027291 

0.060231 

1 

1 

0 

0 

1 

0.015407- 

-0. 030071 

0.018071 

0.038845 

0.035779 

0.010151 

I 

0 

1 

1 

0 

0.012253 

-0.06406: 

0.0092641 

-0.020150 

-0.037781 

0.0072177 

1 

0 

0 

0 

2 

0.28491 

0.34072 . 

0.019744 

-0.19044 

-0.016595 

0.033411 

i 

1 

0 

1 

0 

0.014S9S * 

0.026726 


0.011230- 

-0.014249 


0 

0 

0 

3 

0 

0.015919 

0.023200 


0.0010312 

-0.0000181 


0 

0 

0 

2 

1 

-0.059994 

-0.039295 


-0.00020925, 

0.0020093 


0 

1 

0 

0 

1 

0.23606 

0. 19028 


-0.26170 

-0.16341 


0 

1 

0 

1 

1 

0.12773 i 

-0.0014177 


-0.14794 

0.017456 


0 

0 

1 

2 

0 

0.011025 

0.035926 


-0.036119 

-0.0093768 


0 

2 

0 

1 

0 

-0.029322 

-0.00052429 


-0.010982 

-0.0037828 


0 

2 

0 

0 

1 

-0.035706 

-0.012967. 


0.029166 

0.020917 


0 

1 

1 

0 

1 

-0.094705 

-0.14723 

0.068493 

0.077427 

0.16682 

0.043811 

0 

1 

1 

1 

1 

0.14525 , 

0.077035 


-0.14725 . 

-0.034027 


0 

0 

2 

0 

2 

-0.18694 

-0.043066 

0.21009 

0.35176 

-0.026340 

0.14416 

0 

3 

0 

0 

0 

0.011425 

0.0053026 


0.0010916 

-0.0030360 


0 

1 

2 

0 

1 

0.0026651 

-0.079804 


0,011013 

0.054345 


1 

0 

0 

2 

0 

-0.019882 

-0.022990 


-0.0059394 

-0.00068379 


1 

1 

0 

0 

2 

0.011149 

-0.010324 


-0.038701 

-0.041927 


1 

0 

1 

0 

2 

0.11888 

0.36121 

-0.12376 

-0.35296 

-0.11629 

-0.075821 

I 

2 

0 

1 

0 

0.0069946 

-0.0037171 


0.0013877 

0,0023460 


1 

0 

2 

1 

0 

0.0080525 

-0.013398 


-0.038144 

-0.036215 


1 

0 

2 

0 

1 

-0.11705 

-0.32618 

0.0866S7 

0.3S172 

0.10579 

0.054464 

2 

0 

1 

0 

0 

-0.052368 

.0.049678 


0.0067828 

-0.0054071 


2 

0 

1 

0 

1 

-0.0023228 

-0.076542 


-0.047024 

0.0023889 


0 

0 

0 

2 

2 

0.074855 



-0.039671 



0 

0 

0 

1 

3 

-0.070934 



0.19798 



0 

0 

0 

0 

4 

0.083193 



0.023867 



0 

0 

1 

1 

2 

-0.020887 



-0.011937 



0 

1 

1 

1 

0 

-0.14127 



0.14244 



0 

0 

3 

0 

3 

0.11344 



-0.12287 



0 

2 

2 

0 

0 

0.031675 



-0.044026 



0 

1 

3 

0 

0 

0.092418 



-0.19611 



0 

0 

4- 

0 

0 

D. 083328 . ■ 



0.023846 



1 

0 

0 

2 

1 

0.0048235 



0.0049736 



1 

0 

0 

1 

2 

-0.064391 



0.023862 



i 

0 

0 

0 

3 

-0.039594 


-0.012946 

-0.00032689 


-0.013326 

1 

0 

1 

i 

1 

0.025114 



0.040709 



I 

1 

1 

0 

1 

0.0047199 



0,041727 



1 

0 

2 

0 

2 

0.010608 


-0.11516 

-0.28168 


-0.066636 

1 

2 

1 

0 

0 

0.0022431 



-0.0062230 



1 

1 

2 

0 

0 

-0.040334 


-0.00085946 

0.025800 


-0.00093048 

1 

0 

3 

o , 

0 

0.038687 


0.013830 

0.000079602 


0.015427 

2 

0 

0 

0 

1 

0.046280 



-0.0099149 



2 

1 

0 

0 

0. 

•0.021116 



-0.0076620 



0 

2 

0 

0 

0 



-0.021504 



-0.016327 

1 

1 

0 

0 

0 



- 0.006054o 



0.0077624 

1 

0 

0 

1 

o- 



0.0072332 



0.0057015 
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EXCHANGE AND COULOMB ENERGY OF H 
Table XL The coefficients *’C/ for the wavcfunction V (1) (RS). 
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go «c 



r 

i 

r 

3 

R=4 

R=6 


4 

R=6 

2?=8 

0 

0 

0 

0 

0 

1.0795 

0.47203 

0.34950 

0.46901 

0.11475 

0.23760 

0 

0 

0 

1 

0 

-2.4267 

-1.2169 

—0.89202 

-0.59911 

-0:59469 

-0.64864 

1 

0 

0 

0 

0 

-0.025480 

0.50379 

0.28219 

-0.77282 

0.39560 

0.23453 

0 

0 

0 

2 

0 

0.50808 

0.21295 

0.22413 

-0,50208 

0.063150 

0.16355 

0 

1 

0 

1 

0 

0.75202 

0.24268 

0.12079 

-0.18459 

0.14579 

0.067388 

0 

0 

1 

0 

1 

0.20340 

0.71649 

-0.17025 

-0.59061 

0.18786 

-0. 19483 

0 

0 

0 

0 

2 

-0.87846 

-0.84960 

-0.56734 

-0.70433 

-0.44423 

-0.48116 

2 

0 

0 

0 

0 

-0.29709 

-0.13494 

-0.083826 

-0,40648 

-0.096717 

-0.057912 

0 

0 

0 

1 

2 

0.86379 

0.76862 

0.42837 

0.46724 

0.69628 

0.36722 

0 

1 

0 

0 

2 

1.1605 

0.66784 

0.39023 

0.59328 

0.40433 

0.34132 

0 

0 

1 

1 

1 

0.41920 

-0.99775 

0.36429 

0.65109 

-0.30486 

0.34082 

1 

0 

1 

0 

1 

0.012357 

0.50529 

0.020496 

-0.057532 

0.27165 

0.065340 

2 

0 

0 

1 

0 

0.014939 

0.005643U 

0.027311 

0.18825 

0.047829 

0.019018 

2 

0 

1 

0 

1 

-0.46478 

-0.22085 

-0.023530 

-0.096753 

-0.075474 

-0.015719 

1 

0 

0 

1 

2 

0.012780 

-0.0085408 

-0.12639 

-0.21729 

0.11233 

-0.11726 

0 

1 

0 

3 

0 

0.047405 

-0.022361 

-0.030178 

-0.038742 

-0.0079197 

-0.021086 

0 

2 

0 

0 

2 

-0.45030 

-0.068766 

-0.14566 

-0.39949 

-0.099838 

-0.16418 

1 

1 

0 

0 

2 

0.32757 

0.071650 

0,031499 

0.10156 

0.13111 

0,056315 

1 

0 

2 

0 

2 

-0.60541 

—0.32849 

0.027198 

0,00034635 

0.14522 

0.034654 

0 

0 

1 

2 

1 

0.24711 

0.43031 

-0,060443 

0.026398 

0.11720 

-0,091280 

1 

0 

0 

1 

0 

1.3183 



2.4656 



0 

0 

0 

3 

0 

-0.029792 



0,17192 



0 

1 

0 

2 

0 

-0.32184 



0.41337 



1 

0 

0 

2 

0 ■ 

-0.25421 



-0.62435 



1 

0 

0 

0 

2 

-0.40164 

-0.14295 


0.25896 

-0.41729 


3 

1 

0 

1 

0 

-0.37739 

-0.028536 


-0.76601 * 

-0.073163 


0 

0 

0 

2 

2 

-0.076728 

0.0052492 


0,15498 

-0.090497 


0 

0 

0 

0 

4 

0.017466 



0.12185 



0 

1 

0 

1 

2 

-0.70930 

-0.51504 


-0.43101 

— 0.3218S 


0 

2 

0 

2 

0 

0.0086381 

-0.040965 


-0.065316 

-0.022744 


0 

0 

2 

0 

2 

0.78541 

0.78985 


0:094660 

-0.057522 


0 

1 

1 

1 

1 

-0.39479 



-0.062413 



1 

1 

0 

2 

0 

0.042381 

0.068055 


0,19837 

0.039404 


1 

0 

1 

1 

1 

0.59609 



0.042489 



2 

0 

0 

2 

0 

0.020309 

-0.030833 


0.021284 

-0.026549 


2 

1 

0 

1 

0 

0.064036 



0.0091210 



0 

1 

0 

0 

4 

-0,032020 



—0.094058 



0 

0 

1 

1 

3 

-0.20691 

0.064291 


0.50592 

0.31121 


0 

0 

1 

3 

1 

-0.087062 



-0,039437 



0 

0 

2 

3 

0 

0.040929 



0.069411 



0 

0 

2 

1 

2 

0.061907 



-0.20895 



2 

1 

0 

2 

0 

—0.014838 



-0.012007 



2 

1 

0 

0 

2 

-0.028268 



0.0044078 



2 

0 

1 

1 

1 

0.10815 



0.045803 



0 

2 

0 

2 

2 

0.057438 



0.037841 



0 

0 

2 

2 

2 

0.26402 



0,085560 



0 

1 

1 

1 

3 

-0.094076 



-0.32568 



0 

1 

1 

3 

1 

-0.015374 



”0.00086511 



0 

1 

2 

3 

0 

0.0065240 



0.013058 



2 

2 

0 

0 

2 

-0.0015669 



-0.0084985 



0 

0 

0 

4 

0 


-0.028163 



-0.0092040 


0 

0 

1 

0 

3 


-0.21286 



—0.53014 


1 

0 

0 

3 

0 


0.055989 



0.025084 


1 

0 

3 

0 

3 


0.014101 



0.047087 


2 

0 

2 

0 

2 


0.000080012 



-0.018930 


1 

1 

0 

3 

0 


-0.0037390 



0.0010645 



Tabls XII. Expectation value of the kinetic energy corre- 
sponding to ty(l, A), Here X) is given by Eq. (26a) and 

A-yq, X) Vwi, x)l V(1. M ). 


Table XIIL Scaling parameter s for ‘\&(1» X) and improve- 
ment in <2(1, X) resulting from the satisfaction of the virial 
theorem, A! «=* -lrro- 


°KiX t \) u X(i,x) 


R EL-HAV HS EL-HAV HS 


1.0330 1.036S 
1.0036 1.0028 
1.0002 1.0001 


. f £(i,x) 

** *-£{lj X) nectar 


% + EL-HAV 0.993 7 cm- 1 

HS 0.987 25 cm -1 

**,+ EL-HAV 
HS 


4 0.93973 0.95241 
6 0,99540 0.99612 
8 ‘ 0.99970 0.99970 


1.001 

1.000 


0 

a 
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CERTAIN, HIRSCHFELDER, 
second-order polarization (or dispersion) energy 
22®(pol) = —222.61 cm- 1 . 

Actually, we would not expect a meaningful compari- 
son between the *22® and the 22® (pol) for separations 
smaller than the van der Waals* minimum for the 1<t u 
state which occurs 21 at 2?=12.55c 0 . (For H 2 , the van 
der Waals minimum for the 3 2« + state occurs at R= 
7.85co. Thus for H 2 + much larger separations are re- 
quired for comparable behavior.) Thus, we hope that 
van der Avoird will extend his calculations to much' 
greater separations. 

Our calculations indicate that *22® (EL-HA V) ap- 
proaches £22® (dis) , however. We point out that we 
have evaluated the second-order energy with a trun- 
cated basis set and we cannot state with certainty that 
we have converged to the true value of *22® (EL-HAV) . 
It may be necessary to use a set of terms in the basis 
set which has the behavior of continuum wavefunctions. 
Thus we can only say that it appears that the Eisen- 
schitz-London or HAV second-order energy is asymp- 
totically equal to one-half of the second-order disper- 
sion energy. 

B. Arguments that i E (2) (EL-HAV)^lE (2) (dis) 

In a companion paper 22 a model for exchange forces 
is constructed which consists of a particle moving in 
the one-dimensional potential 

U(x) = (8JP)“*(*- J?) 2 0H-jR) 2 . 

The potential U(x) roughly corresponds to the II 2 + 
potential, with the nuclear Coulomb attraction re- 
placed by a harmonic oscillator attraction. 

In a manner completely analogous to the molecular 
case, the application of the various perturbation 
schemes to the model leads to an expansion for the 
energy consisting of a power series in Rr l plus terms 
exponentially decreasing in R 2 , In this case the sum- 
mations in (27) may be carried out exactly and it is 
found that 

22® (dis) « - (11/32 R?) - (21/512 2?) +• • • , 


51 J. M. Peek, Sandia Cotp. Research Rept. SC-RR-65-77 
(1965). 

15 P. R. Certain, J. Chem. Phys. 49 - , 35 (1968), following article. 


KOLOS, AND WOLNIEWICZ 34 

and 

*22® (EL-HAV) (11/64 R 2 ) - (141/10247?) +• * • 

Thus the exchange series in (27) is of higher order in 
Rr l than the direct series for this model, so that 
*22® (EL-HAV)~|22® (dis) . 

APPENDIX: THE VIRIAL THEOREM 

If a scaling parameter s is introduced into the trial 
function *$(1, X) and optimized, the resulting energy 
*0(1, X) satisfies the virial theorem. 23 ’ 24 The equation 
which determines s is 

2^(5) +sV(q) +q{dE/dq) - 0, (29) 

where q—sR, K(q) and V (jq) are the expectation 
values of the kinetic and potential energy, respectively, 
computed with for R—q, and 

E=fK(q)+sV{q). (30) 

For simplicity we have deleted the superscript i desig- 
nating the symmetry of E. 

The last term in (29) may be estimated by 

q{dS/dq) = [(£-co)/(E-e 0 )]2?(^/d22) , (31) 

where E and R(dE/dR) are accurate variational 

values. 16 Since s is nearly unity for internuclear dis- 
tances considered, we can solve Eq. (30) by expanding 
(29)-(31) in power series in (1— s) through quadratic 
terms, holding q fixed at ?o( = 4, 6 or 8c 0 ) . This yields 
the scaled energy at R—qo/s. To obtain the scaled 
energy at R=qo, we expand this result in a power 
series in [2?— (<?o/s)] through linear terms, assuming 
that s remains constant over the necessary interval. 

We used this procedure to estimate s and the improve- 
ment in *22(1, X) (X optimized) at the three intcr- 
nuclear distances considered in this paper. At 22—6, 
8cr 0 we obtained s= 1 for all schemes; i.e., the unsealed 
energies satisfy the virial theorem. The results for R= 
4 g 0 are given in Table XIII. We conclude that the 
energies *22( 1, X) „ are not significantly improved by 
scaling. 


23 J. O. Hi rschf elder and J. F. Kincaid, Phys. Rev. 52, 658 
(1937). 

* P.-O. Lowdin, J. Mol. Spcctry. 3, 46 (1959). 
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Addendum to "Exchange and Coulomb Energy of Determined by Various 

Perturbation Methods". 

III. Previous Perturbation Calculations fpr Hg 

Other recent perturbation treatments of are the single center 

expansion calculations of Dvoracek apd Horak,*^ who obtained a dls- 
sociatiptv energy of 4.207 eV ;> and. the calculation of Kirtman and 
pecius,~^ which employed the same approach as Matcha and Byers Brown 
with similar results. 


IV. Discussion of Results . 

22 

As Her.ring has emphasized', the Heitler-London wave function 

/ \ ^ <j > 0 for the hydrogen molecule predicts the impossible 'result 

that the 3 Z*£ energy crosses the * 2 ? energy at large '(k ~ 50a- ) 

<? Q - 

internuclear separations^ The physical , reason -for this is easy to 
understand in terms of electron- .correlations : the electrons are 

prevented by the' "exclusion principle from getting too close to each 
other in the Heitler-London wavef unction 'for fhe triplet- stated but 
the singlet wavefunction gives a small, but finite, probability for 
the electrons to coalesce. This results in an error in the Heitler- 
London energy "for the singlet state whfch is greater than the error 
for the triplet state. At large distances, where the error is greater 

tb^n the splitting between the levels, a crossing occurs. As Table VIII 

‘ . 29 

shows, and as has .been found by Alexander and Sa.lem, accurate 


numerical values of E(exch ) can be obtained by a second order calcula- 
tion, ever} though it is not likely that the perturbation theory (in 
low order)* gives the correct analytical R-dependence of E(exch ) . 
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A. Arguments that Z E ^ (EL-HAv) ~ B tX> (eLis). 

31 

Van der Avoird has extended his calculations to R - 20 a . 

o 

At R = 15 a Q , he also observes the* second order energy being, too small 
by a factpr of two, but suggests that itr could be due to a lack of 
continuum functions .in his basis set. 

B . Arguments that 1 B C&1~ "V ^ 

-J** 

In Section 4.3, the double minimum delta-function, model for H_ 
is solved, exactly for the various schemes. In this case it. is found 
that 

E^Cdis) = - e.~** +. 0(e' 3/z ) 

and 

*E°*t£irmO = +- 

Thus, for large distance, where €- ^ can be neglected, 

l E- m tEL- mv) ~ ± E*> CdisX 

Thus, although it has not been proved < that the EL-HAV second order 
energy approaches one-half the second order dispersion energy at large* 
distances, this discrepancy has appeared in* every application thus 
far reported. The possibility of modifying the El-HAV procedure to 
remove the discrepancy is considered in Appendix D. 
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Summary of Exchange Perturbation 
Calculations for H 2 at Large Separations* 

JOSEPH O. HIRSCHFELDER AND PHILLIP' R. CERTAIN-]- 
The University of Wisconsin Theoretical Chemistry Institute, Madison , Wisconsin 
and 

The University of Florida Department of Chemistry, Gainesville, Florida 


Abstract 

This paper is a summary and appraisal of the recent work of the authors together with 
Kolos and Wolnicwicz which applied four types of perturbation procedures to the cal- 
culation of the energy of the gtound state and first excited state of the hydiogen molecule at 
large separations. All of the perturbation techniques gave good values for the exchange 
energy but less accurate values' for the Coulombic energy. The reasons why the second- 
order Eiscnschitz-London (or van dcr Avoird or Hirschfcldcr hav) energy approaches 
one-hall of the correct limit at large separations are discussed. There seems to be no unique- 
ness to thp orders in an exchange perturbation problem. 


'There is considerable interest in the development of perturbation procedures 
for exchange forces where the molecular wave function and the basis set ftave differ- 
ent symmetries. Different formalisms give different results for the second- and 
higher-order energies. The mathematical problems have been discussed elsewhere 
[1]. A basic difficulty is that the order of the perturbation terms is not uniquely 
defined. Thus, if A is the operator which projects the component with the 
symmetry of the desired molecular wave function, then A commutes with the molec- 
ular Hamiltonian H. However, A does not commute with either the' zeroth- 
ordcr Hamiltonian H 0 or the perturbation V = H — H Q . Thus, wc have the 
seeming-paradox, 

(1) AfI 0 -H 0 A= VA-AV 

Zeroth order First order 

In any.conventional perturbation scheme, one would suppose that the leftrhand 
side of-this equation is zeroth order and the right-hand side is first order. By the use 
of (1), the apparent orders of various terms in a perturbation expansion can be 

* This work was supported in part by National Aeronautics and Space Administration Grant 
NsG-275-62 to. the University of Wisconsin. 

t National Science Foundation Graduate Fellow. 
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arbitrarily shifted. How subtle are these shifts in the order may be seen from the 
following derivation. 

Let Q 0 be the zeroth-order wave function and let the set of functions Q 0 , Q x , 
* * * span the Hilbert space of H. The exact wave function can then be written 

(2) ¥-«k + 2>iGi 

MO 

The Schrodinger equation becomes 


(3) ' (H - E)Q 0 + 2 - B) Q, « 0 

Mo 

Multiplying (3) by a function %* and integrating, 

( 4 ) e — — i y 

fa I <?0> mo fa I Qo> 

But multiplying (3) by Qf and integrating gives 

(si * ^ MEeiIM _v 

A <Q/I H — E (Q,j 

. , JIMO 

Substituting (5) into (4),' and iterating, 

(6) E = + E& + EM q- . . • 


^•co) ^ fai ^ fal ^ i@o) 

(7) <x\Qo> ’ <x\Qo) 

MO fa I QoXQjl ElQi) 

The convergence of the sequence depends upon the non-diagonal elements 
(Qs I H — E JQ*) being small compared with the diagonal elements U ~ E 
i Qi)' Ir X =* Qo> (6) is the familiar Brillouin-Wigner series. For the present 
examples, let Q 0 = Af 0 and Q j ~f s for j ^ 0, where 

C 8 ) (#o - *o)/o - (#o - e,)/ # « 0 

The lowest orders of the el-iiav (Eiscnschilz-London [2], van dcr Avoird 
[3], or hav[1]) perturbation procedures correspond to setting# = Q 0 = Af 0 , Then 

EM ss € £»> </ol ^ i4/o) 

(9) </• I 4f%) 

B!U„ - -! rnSL rgg! gyw- B. '/•>■+ pip, 

MO </o 4/q>fo — *o) * 
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Here 0( V s ) are terms which appear to be “third order” or higher in the pertur- 
bation. In any case, the 0(V Z ) terms do not appear in the el-hav second-order 
energy. 

The lowest orders of the ms-ma (Murrell-Shaw [4] or Musher-Amos [5]) 
procedures correspond to setting % =f 0 . In this case, f? (0) and J S <1) are the same 
as before, but 


( 10 ) 


i*0 


(fo\ V\f t )<f s \A(V-m l/o) 

</o I 4fo)( e j ~ 6 o) 


+ 0'(V*) 


Here the 0'(V 3 ) apparent “third order” and higher terms do not appear in the 
ms-ma second-order energy. 

Although the derivation of the -Eel-hav an ^ the -Ems-ma are ver Y similar, 
the resulting formulae are quite different. The difference between these two 
“second-order” energies can be traced to “third-” and higher-order terms in the 
perturbation sequences. 

There have been many other types' of perturbation schemes proposed for 
calculating exchange forces. One of the most curious is the hs (Hirschfelder- 
Silbey [1,6]) scheme which proposes that there exists a physically significant 
“primitive function” whose symmetry projections correspond to the wave functions 
for the family of quantum states arising from a single-electron configuration. The 
•function f 0 is the zeroth order of the primitive function. 

Intuitively, the ideal perturbation treatment would seem to be the straight 
Rayleigh-Schrodinger (rs) procedure. For one-or two-electron problems where the 
wave function separates into a space function times a spin function, Sternheimer [7a] 
showed how to form a satisfactory zeroth-order Hamiltonian. If K is the kinetic 
energy operator and Af 0 is the zeroth-order wave function, then the “Sternheimer 
Hamiltonian” is 


(H) 


J3J-JC + 


gjg ] 
Afo -I 


Clearly, HqA/ 0 = q>4/Jj . The Sternheimer Hamiltonian together with Af^ forms 
the basis for setting up a Rayleigh-Schrodinger perturbation procedure. The 
Sternheimer Hamiltonian is the only possible zeroth-order Hamiltonian corre- 
sponding to the zeroth-order wave function Af 0 which has a local potential. 
Unfortunately the Sternheimer Hamiltonian cannot be generalized conveniently 
to many-electron systems [7b]. 

Since there is no way of telling a priori from the formalism which of the various 
perturbation schemes gives the best second-order energy or gives the best expecta- 
tion value of the Hamiltonian corresponding to the wave function accurate through 
the “first” order, we must resort to mathematical experimentation. The hydrogen 
molecule would seem to be the best “guinea pig”. Recently, Certain, Hirschfelder, 
Kolos, and Wolniewicz [8] have made calculations for both the' *2+ ground state 
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and the 3 £+ first excited state using a basis set of Coolidge and James type func- 
tions (70, 50, and 30 terms for R — 4, 6, and 8a 0 , respectively, for the el-hav, 
hs, and ms-ma procedures; and 50, 35, and 20 terms for R « 4, 6, and 8a 0 , 
respectively, for the rs). The first-order wave functions were determined by 
varying the coefficients of these basis functions so as to satisfy a variational prin j 
ciple. The zeroth-order wave functions are the usual Heitler-London expressions. 


Table I. Accurate variationally determined energies of the ground state (g) and the first 
excited state (u) of the hydrogen molecule. i? Coul and E cxch are given by (12) and (13). 


*/« 0 



e o — -®coui 

■^exch 

4 

3592.6 cm- 1 

1453.3 cm” 1 

1069.6 cm- 1 

2523.0 cm" 1 

6 

178.9 

41.2 

68.8 

110.1 

8 

11.6° 

—4.3 

8.0 

3.6 


° In the new variational calculations given in [8] the improved value of 11.74 cm” 1 was 
obtained. 


Table I gives the accurate binding energy of the ground state *2* of the hydrogen 
molecule, <r 0 — S E } and the accurate energy of interaction of two hydrogen atoms 
in the first excited 3 £+ state of the hydrogen molecule, n E — e 0 . These values 
were previously determined by Kolos and Wolniewicz [9] using the Rayleigh-Ritz 
variational principle together with a very large basis set. The corresponding 
accurate Goulombic and exchange energies are respectively [10], • 

(12) - ^coui = H(*o " *E) - ( u £ “ *)] 

(13) -£exch - *[(«o - g £) + ( u 2? - e 0 )] 

It should be remembered that e 0 — g E = 38293 cm -1 at the equilibrium separation 
R = 1.4 a 0 . Also, the minimum in the a E — € 0 occurs at R — 6.85a 0 . 


Table II, Contribution of the first-order energies to the interaction, 
Coulombic and exchange energies. 


R/a 0 

■£<>>/(«„-'£) 

“£<«/(“£ - «o) 

1 

$£ 

o' 

1 

f 

•®exch / ^exch 

4 

0.69 

1.33 

0.26 

0.87 

6 

0.62 

2.00 

0.10 

0.00 

8 

0.33 

0.70 

0.05 

0.02 


Table II shows the fraction of the interaction energy which is given by the 
first-order contribution. Here, in accordance with (9), e 0 -j- %E W and e 0 -f- u ii {1) 
are the expectation values of the Hamiltonian for the molecule using the Heitler- 
London wave functions for the *£+ ground state (g) and the 3 £* first excited 
state (u), respectively. The first-order Goulombic and exchange energies, 
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and E™ h are given by (12) and (13) with the first-order energies replacing the 
total energies of the system. This definition of Coulombic and exchange energies 
is slightly different from the one which chemists have used in connection with 
Heitler-London wave functions. Writing 

(14) e 0 + 1 £ w = (J±L)j(l±S) 

the chemists have regarded J and L as approximations to the Coulombic and 
exchange energies. With our definition, the corresponding Coulombic and 
exchange energies are ( J — SL)j(l — S z ) and (L — — 5 2 ), respectively. 

The two definitions differ since the overlap integral S is not zero. There is a 
considerable body of evidence from the theories of many-electron systems, the 
solid state, and ferromagnetism which shows that our definition is the more 
meaningful [1 1]. At R = 4, 6, and 8 <j 0 the value ofiS is 0.0357, 0.0022, and 0.0001, 
respectively. Because the value ofS is small, our conclusions would not be changed 
if we had changed our definition of E cxch and fs Coul . 

From Table II it appears that the interaction energy as determined by the 
first-order energy is, for the ground state, too small and, for the first excited state, 
too large. The first-order energies give only a small fraction of the Coulombic 
energy and indeed — £<*>/( e 0 — £ Coill ) should approach zero at large separations 
where the Coulombic energy becomes equivalent to the second-order dispersion 
energy. It is indeed surprising that E exch gives 90 % of the exchange energy with- 
out any appreciable variation with the separation. 

Tables Ilia and IV give the results of the Certain, Hirschfelder, Kolos, and 
Wolniewicz calculations corresponding to various types of perturbation procedures. 
Let us examine separately each of the subsections of Table III. 

(a) Ground-state energy through second order . 
rs is very good, best at R = 4 a 0 . 

hs and ms— ma are both good and give virtually the same values, best at R — 8c 0 . 
el-hav is good at R — 4 a 0 , and becomes bad at R — 8a 0 . 

(b) Excited-state energy through second order . 
rs is not good at any separation. 

hs and ms-ma again give virtually the same values, slightly better than rs at 
R sss 4 and 6a 0 but still not good, and very good at R = 8a 0 . 

el-hav is the same as hs and ms-ma at R = 4, but is very bad at R — 6 and 
8a 0 . 

(c, e, and Table IV) Expectation value of the Hamiltonian for the ground stale wave 
function through the first order. 

rs is good at R = 4 and 6a 0 but poor at 8a 0 . The energy is almost unchanged 
when A is varied since the optimum A is nearly unity. 



89 


130 JOSEPH O. HIRSCHFELDER AND PHILLIP R. CERTAIN 

Table III. The energy through the second-order, V 2) ; the expectation value of the 
Hamiltonian corresponding to the wave function through the first order, *2?(1, 1); the 
l£(l, 1) improved by optimizing the normalization of the first-order wave function, ^isl, A); 
and the corresponding Coulomb and exchange energies. 


* 1*0 

EL-HAV 

IIS 

MS-MA 

RS 

EL-HAV HS 

MS-MA 

RS 

(a) 

(«0 -'«<“)/(«» - 

*E) 


(b) 

(V 2 > - e 

o)/CE ~ <o> 


4 

0.88 

0.88 

0.89 

0.92 

1.12 

Ml 

1.12 

1.14 

6 

0.86 

0.97 

0.97 

0.97 

1.57 

1.09 

1.09 

1.14 

8 

0.67 

0.99 

0.99 

0.97 

0.12 

0.99 

0.99 

0.90 

(c) 


> !))/(*„ -*£) 


(d) 

mu i) 

” W * - 

' c o) 

4 

0.93 

0.92 

0.91 

0.95 

1.05 

1.09 

1.10 

1.15 

6 

0.92 

0.96 

0.85 

0.96 

1.33 

I.1I 

1.57 

1.17 

8 

0.84 

0.99 

0.67 

0.88 

0.53 

0.93 

0.09 

0.77 

(c) 

(«*-**a 

, *))/(« 

o ~ **> 


<0 

mu $ 

- * 0 )ICE - 

- e o ) 

4 

0.95 

0.92 

0.91 

0.96 

1.03 

1,09 

1.10 

1.15 

6 

0.95 

0.96 

0.85 

0.96 

Ml 

1.11 

1.57 

1.17 

8 

0.99 

0.99 

0.67 ‘ 

' 0.89 

0.93 

0.93 

0.10 

0.77 

(g) 

( £ o — £ coui)/( e o “ 

‘ -^Coul) 


(h) 

£ e\ch/^cxch 


4 

0.72 

0.73 

0.74 

0.77 

0.95 

0.95 

0.96 

0.98 

6 

0.64 

0.93 

0.93 

0.92 

0.99 

0.99 

0.99 

1.00 

8 

0.53 

0.99 

0.99 

0.95 

1.00 

1.00 

1.00 

1.00 

0) 


1) Co ul)/ ( c 0 ^coul) 

CD 

*®(^» Oexch/^exch 


4 

0.85 

0.80 

0.77 

0,82 

0.97 

0.96 

0.96 

1.01 

6 

0.80 

0.91 

0.63 

0.90 

1.00 

0.99 

0.98 

1.00 

8 

0.76 

0.97 

0.51 

0.85 • 

1.02 

1.03 

1.00 

0.95 

GO 

<*-*(1* 

^)coul )/( c 0 ^coul) 

(1) 

E(U ^)exeh/^cxch 


4 

0.90 

0.81 

0.78 

0.83 

0.97 

0.97 

0.96 

1.02 

6 

0.90 

0.91 

0.63 

0.90 

0.98 

0.99 

0.98 

1.00 

8 

0.97 

0.97 

0.53 

0.86 

1.03 

1.02 

1.00 

0.96 


Table IV. Optimum value of the parameter A which optimizes the energy *£(I, A). 


* 1*0 


Ground state 



First excited state 


EL-HAV 

HS 

MS-MA 

RS 

EL-HAV 

HS 

MS-MA 

RS 

4 

1.36 

1.20 

1.16 

1.16 

1.42 

1.08 

1.09 

0.94 

6 

1.39 

0.97 

0.96 

0.98 

1.89 

0.98 

0.99 

0.97 

8 

1.92 • 

0.99 

0.97 

0.90 

1.98 

0.96 

0.97 

0.94 
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HS is not as good as rs when R ~ 4a 0 , equal to 6% and much better at 8 a Q . 
Altogether hs is quite good at all separations. The optimum A is nearly unity $o 
that the energy is not improved by the variation of X. 

ms—ma is almost as good as hs for R = 4a 0 but is very bad at 6 and 8a 0 . Again, 
X is nearly unity and the energy is not improved by variations of X. 

el— hav with X — 1 gives values lying approximately half way between hs 
and ms— ma; it is good at JR. « 4e 0 , fair at 6a 0 , and poor at 8 a 0 . Variations of X 
make a considerable improvement in the energy, especially at R « 8^ 0 where the 
optimum value of X is almost two. The X — 2 at large separations corresponds to 
the second-order energy being one-half of the dispersion energy. For the opti- 
mised values of A, the el-hav energies are very good at all separations. 

(d, f, and Table IV) Expectation of the Hamiltonian for the excited state wave function 
through the frsi order . 

rs is poor at all separations and not improved by variations in X since the 
optimum values of X are nearly unity. 

hs is fair at all separations and not improved by variations in X since the opti- 
mum values of X are nearly unity. 

ms-ma is fair at R = 4« 0 but very bad at larger separations. These energies 
are not improved by variations in X since the optimum values of X are nearly 
unity. 

el-hav with A = 1 gives a good energy at R s* 4% but bad energy (between 
hs and ms-ma) at larger separations. Again, the el-hav energies are considerably 
improved by variations of A. At R « 6 and 8s 0 the optimum value of A is very 
nearly equal to two, corresponding to the second-order energy becoming at large 
separations too small by a factor of one half. With the optimum values of A, the 
EL-HAV energies are good. 

(g, i, and k) Goulombic energy, 

rs is poor for R = 4<? 0 and somewhat better for 6 and &z 0 . The results for 
the energy through the second order are somewhat better than is obtained with 
the expectation values for the Hamiltonian at 6 and 8a 0 . 

hs Coulombic energies are comparable to rs at R *s 4 and 6 , but become 
good to very good at , Again the energy through the second-order gives better 
results than the expectation value of the Hamiltonian at 6 and 8% . 

ms-ma Coulombic energy through the second order is virtually the same as for 
the hs; however, the expectation values of the Hamiltonian give much worse 
results, especially at large separations. 

el-hav gives very poor Goulombic energy through the second order. With 
A = 1, the expectation values of the Hamiltonian give poor results. However, 
with the optimum value of A, the Coulombic energies obtained from the expecta- 
tion values of the Hamiltonian are good at all separations. 
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(h, j, and 1) Exchange energy. 

All four perturbation procedures give very good values of the exchange energy 
at all separations. The values obtained using the energy through the second order 
are virtually perfect at R — 6 and 8a 0 . 


Conclusions 

We reach four general conclusions. 

[1] . Of all the schemes which we considered , the Hirschfelder-Silbey procedure gives the 
most uniformly good results. Nevertheless, we.are not enthusiastic about this method 
because it is difficult to understand the physical significance of the “primitive 
function” whose projections give the exact wave function for both the singlet 
(g) and the triplet (u) states. Furthermore, the variational principle for the first- 
order wave function minimizes 

</ol A I /„> «£<•■» + Ul AJ/J 

rather than *E i2) and u is ,2) separately. 

[2] . All of the perturbation procedures give good values for the exchange energy but not 
as accurate values for the Coulombic energy. Table II shows that most of the exchange 
energy is given by the first order but very little of the Coulombic energy is given 
by the first-order terms. One would expect that the Coulombic energy should 
agree accurately with the second-order dispersion energy at R — 8 a 0 and the 
agreement should still be fairly good at R — 6 and 4c 0 . The second-order disper- 
sion energy is 

(15) - I K/ol V I f,)\ V(e 0 - «=,) 

i* 0 

In devising a new type of perturbation scheme for exchange problems, it would 
be desirable to require that the £ Coul ec iual to the E aisp . One might suppose 
that the error in is Coul might be due to the basis set which was used rather than 
in the perturbation methods themselves. However, this notion is dispelled by 
Table V where it is shown that using the Rayleigh-Ritz variational procedure with 


Table V. Energies (designated by tilde) determined by Rayleigh-Ritz variational pro- 
cedure using the name basis Art ns was used for the kihuay, ns, one! MS-MA perturbations, 
compared to accutatc energies of Table 1. 


RK 

1 1 

I J 

g 3 

1 1 
o o 

I 1 

4 

1.00 

1.01 

0.99 

1.00 

6 

0.99 

1.06 • 

0.97 

1.01 

8 

1,01 

1-.00 

1.00 

1.01 
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the same basis sets gives extremely good values ofboth the Coulombk and exchange 
energies. Thus, we would get about the same order of accuracy for the energies of 
the singlet and triplet states if we followed the suggestion of Robinson [12] and 
took the exchange energy to be isg *k and the Coulombk energy to be * 
This would have the advantage of greatly simplifying the calculations. Of course* 
the precision would not be improved. 

[3]. At targe separations^ the second-order EE-hav energy seems to approach, me half 
the Cmlmnbk {or di$perswn)ewrgy* This is seen most -dearly ia Table VI. All of the 


Table YL Comparisons of the second-order perturbation energies with 
the accurate Coulombic energy. 


R/% 


-•mu, , 

*“* w 




*“ -dcoal) 


EL-HAV 

MS 

MS-MA 

ns 

EL-HAV 

m 

WHMA 

RS 

4 

0.64 

0,65 

0.69 

0.78 

0.29 

0.30 

0.28 

0.26 

6 

0.60 

0.89 

0.89 

0.89 

0.31 

0.60 

0.80 

0.57 

8 

0.50 

0.96 

0.96 

0 sz 

0.44- 

0,91 

0.91 

0.86 


other perturbation methods gives values of and — n £' {2 * which approach 

{% — 2? Csni ) at large separations. Closely related is the fact (as shown in Table 
IV) that the optimum value of X for the EWIAV approaches two at large separa- 
tions whereas the optimum value of X for the other perturbation schemes is dose to 
unity; Let us consider more closely the behavior of £ (2) (el-kav) as given by (9). 

For the singlet (g) and the triplet (u) states we have the projection operators 
A ff as (1/4) (I + Z)(l 4* /*$} and A v *» (1/4)(1 — l)(l — P 12 ) where permutes 
the electronic coordinates and /.inverts the wave function through the midpoint 
of the internuclear axis. Thej^ «* ® *he product of a hydrogen 

orbital for electron W I M about center < V ? times another hydrogen orbital for 
electron M 2 W about center **&**. We will let be the function f with electrons 
*T” and **2*’ interchanged. 

At large separations* there are simplifications in (9) which can be made since 
the second-order energy varies as and we can neglect terms which decrease 
exponentially with the separation. However* we must be careful since an infinite 
sum of terms* each of which decreases exponentially with. R, may vary as a power 
of 1/R, Before writing down the asymptotic expression for the second-order 
energy we should note that the overlap integrals, {f 3 j/^), the first-order energies, 
8 R a) and and the second-order exchange energy ail decrease exponentially 
with it Furthermore, Table IILh assures its that the ex.~hav exchange energy 
through Ihe second order behaves properly at large separations. Thus, at large 
separations, to the accuracy which we require, ff i^ 2> *» ix J2 {2) «* * Using (9) 
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to obtain the equation for E$ u] and expressing its limit as R becomes large, 

(16) £(*>(e^-hav) = p&>„ - 12 J -& 1 - * 

i^O € i — e 0 

Each term in the summation decreases exponentially with R. However, without a 
detailed investigation we cannot reach any firm conclusions about the behavior 
of the sum. It certainly seems unlikely that at large separations £ (2 > (el-hav) 
would have the correct limiting form equal to Z^fJ p . Our calculations for H 2 
indicate that it does not. 

Van dcr Avoird [I3a] has calculated for H* the values of 5 fs (2 > and l £ (1, 1) 
using his method (which is equivalent to the el-hav scheme). For R < 7 a 0 he 
obtains excellent agreement with the exact values. For these separations the ex- 
change energy is larger than the Goulombic energy, however, and wc would not 
expect a meaningful comparison between l E i2) and Ej^ for separations smaller 
than the van der Waals minimum for the excited state, which occurs at 12.55o 0 , 
At our suggestion, van der Avoird [13b] has extended his calculations to distances 
up to 20a 0 and finds that for large R , the second term in (16) vanishes with respect 
to the first. He also finds that highly excited states make a substantial contribution 
to the sum in (16) for large R, suggesting that accurate values are difficult to obtain 
with a truncated basis set. 

For an analogous double minimum one-dimensional potential, Certain [14] 
was able to evaluate the corresponding expressions rigorously and in this case he 
found that the sum in (16) was smaller than 7?^ by a factor of the order of R~ 2 . 

[4], If the perturbation schemes that we have considered are to ham practical value in 
determining energies, new methods must be developed for solving perturbation equations 
without recourse to the use of large basis sets [15]. Table III shows quite clearly that 
third- and higher-order contributions to the energy must be included to obtain 
accurate results. The effort involved in such a calculation would certainly be 
greater than a straightforward Rayleigh-Ritz variational calculation if it were 
necessary to employ large basis sets and variational principles to compute each 
order of energy. 

In the various schemes for perturbations involving electron exchange, for 
properties other than energy, on'e still has * ! Dalgarno-type” interchange theorems. 
However, the interchanged perturbation equations are themselves extremely 
difficult to solve. 

Further research is required to develop new techniques for solving perturbation 
equations or to establish a procedure for exchange problems which is more satis- 
factory than the ones which we have considered. 
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One-Dimensional Model for Exchange Forces* 
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Four different types of perturbation theories for the exchange forces between two atoms are applied to the 
lowest gerade and ungcrade state of a simple one-dimensional model. It is shown that the second-order 
energy in one of the expansions (the Eisensclutz-London, van der Avoird, or Hirschfclder HAV expansion) 
exhibits undesirable behavior in the limit corresponding to large internuclcar distances. 


I. INTRODUCTION 

The formal 1-7 and computational 8-12 aspects of the 
development of a satisfactory perturbation theory for 
exchange forces between "two atoms have received a 
great deal of attention recently. In an application 11 of 
three different types of perturbation treatments (EL- 
HA V, HS, and MS-MA) to the interaction of two 
hydrogen atoms, it was discovered that the EL-HAV 
(Eisenschitz-London, van der Avoird, or Hirschfelder 
HAV) expression for the second-order energy apparently 
fails to agree with the second-order dispersion energy 
(no electron exchange) at large internuclear separations, 
where exchange forces are negligible. No-rigorous state- 
ment regarding this discrepancy could be made, how- " 
ever, since the second-order energy was not evaluated 
exactly. 

In the present paper, we apply the following formal- 
isms to lowest gerade and ungerade states of a simple 
one-dimensional model: 

RS: An unsymmetrical Rayleigh-Schrodinger ex- 
pansion. The second-order energy is analogous to the 
second-order dispersion energy of the H-H interaction. 

EL-HAV : Eisenschitz-London, 1 van der Avoird, 5 or 
Hirschfelder HAV 6 expansion. These expansions give 
the same expression for the second-order energy. 

HS: Hirschfelder-Silbey expansion. 2 

MS-MA: Murrcll-Shaw 3 or Musher-Amos* expan- 
sion. These two formalisms arc equivalent. 7 
This model calculation has the advantage that the 

* This -work received financial support from the National 
Aeronautics and Space Administration Grant NsG-275-62. 
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EL-HAV second-order energy can be evaluated exactly. 
The discrepancy with the analogous “dispersion energy” 
is verified for this model. The results obtained here are 
also relevant to the question 6 : “Which formalism pro- 
vides the best second-order energy and the best expecta- 
tion value of the Hamiltonian using the wavcfunction 
accurate through the first order?” As in the application 11 
to the hydrogen-hydrogen interaction, we do not find 
any of the four perturbation schemes outstandingly 
better than the other three. 

H. A SIMPLE MODEL FOR EXCHANGE FORCES 

We consider the one-dimensional system described 
by the Schrodinger equation 

-.w( X )+u(x)t( X ) =m*), CD 

where 

U(x ) - (1/8R 2 ) (*-R) 2 (s-l-R) 2 . (2) 

The potential U(x) has minima at x=±R and a maxi- 
mum at x—0 of height R 2 /S. As R increases without 
limit, U(x) evolves into two simple-harmonic-oscillator 
(SHO) potentials centered at =bi?. Thus, at J?=«> the 
lowest state is doubly degenerate with E=\. As R 
decreases, this state is split into a gerade and an un- 
gerade state under the influence of two physical ef- 
fects. The first, a long-range effect due to the lowering 
of the barrier height, increases the space available to a 
particle in the lowest stale and decreases its energy. 
The geradc-state energy is affected more since sym- 
metry forces the ungerade-state wavefunction to have a 
small amplitude beneath the barrier. The second effect, 
a short-range effect due to the decrease in separation 
between the wells, restricts the particle to a smaller 
region and increases its energy. 

These considerations indicate that the energy of the 
lowest state decreases initially as R decreases from 
infinity. This stale splits into a g and a u state, with the 
g slate decreasing in energy more rapidly than the 
u state. Finally, at short range the energy of both slates 
increases to infinity. This is the same qualitative be- 
havior as is found for the lowest states of H 2 + and H 2 , 
and is the reason for the choice of the model. 

In the discussion which follows, will denote the 
exact energy of the lowest gerade (i=g) and ungerade 
(i=u) state. Although no closed form solution exists 
for (1), the method of Secrcst, Cushion, and Hirsch- 
35 
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felder 18 can be used to obtain f E to the required ac- 
curacy. 

III. RAYLEIGH-SCHRODINGER EXPANSION 

For R large the minimum in U(x) at x——R is a 
small perturbation of an SHO potential at x-R . To 
show this we change variables to q-R—x. Then (1) 
becomes 

(ffo+F-JSWd-Q, (3) 

where 

(4) 

V=-(i/2R)<?+(l/8R>)q*. (5) 

The eigensolutions of Ho are (4>l, where 

(V5- (6) 

and the H k ( q ) are Hermite polynomials. . 

The natural perturbation parameter is RrK This 
choice does not correspond to the symmetrical expan- 
sions, however, so instead we expand in powers of the 
total perturbation V. The usual formulas of Rayleigh- 
Schrodinger perturbation theory 14 yield 

£<»(RS) — 3/32P 2 , 

(RS) - - (11/32 R?) - (21/512P 4 ) , (7) 

and 

^ (0) (RS) 

mw-oox 1 ^ \ 1/3 \/2 , 6 1 ' 2 \ 

* (RS)_ J?\T* + 12^7 ^\32'* + 64 *)' 

This approach does not recognize the inversion 
symmetry (#->--#) of the total Hamiltonian and pro- 
duces a wavefunction which has neither gerade nor 
ungerade symmetry. A simple way to obtain such a 
function to act on ^ 0) (RS) -{-^(RS) with operators 
which project onto gerade and ungerade function 
spaces. These operators are 

A.=i( 1+2), 

A .=1(1-2), (8) 

where I is the inversion (x— >— x) operator. In this way 
we obtain a trial wavefunction 

^ (RS) = A (RS) +hA tfW (RS) , (9) 

where X,- is a variational parameter, and use it to com- 
pute the expectation value of the total Hamiltonian- 

•E(RS) - (tf(RS) I n<f(RS) )/<<?(RS) I V-(RS) >. 
( 10 ) 

15 D. Secrcst, K. Cashion, and J. O. Hirschfeldcr, J. Chem. 
Phys. 37, 830 (1962). 

11 J- O. Hirschfelder, W. Byers Brown, and S. T. Epstein, . 
Advan. Quantum Chem. 1, 265 (1964).. 


It has recently been shown 16 that the resolvent oper- 
ator for a Hamiltonian which is the sum of the SHO 
Hamiltonian and a polynomial perturbation \P(x) of 
degree exceeding two is not expressible in a convergent 
power series in X. Thus the energy is not analytic at 
X — 0. This is true even if the total potential -^-{-XPfc) 
is bounded from below, since in this case the potential 
\x 2 ~\P(x) is not bounded. The RS scheme is of this 
type for the model considered here, and will therefore 
not converge. Nevertheless, the wavefunction (9) is a 
valid trial function and (10) yields an upper bound to 
the true energy. 

IV. SYMMETRIZED EXPANSIONS 

The general structures of the EL-HAV, HS, and 
MS-MA formalisms have been fully treated else- 
where. l ~ 7 The zeroth- and first-order energies and the 
‘ zeroth-order wavefunction are the same in. all three 
formalisms: 

= (3/32JP)=Ff[(l+R ! )<- ,! 7(l±«- Ji! )3 1 (11) 

where the upper (lower) sign corresponds to i~g{u) } 
and 

o. ( 12 ) 

It is in expressions for the second-order energy and 
first-order wavefunction that we encounter differences. 
The most direct way to compare the results is with the 
spectral expansions of *E® and ty (1) : 

<£(2)= (0 O j (7 >o | <f> k ) {4>k I •S'*0o>/(€o— <fc)> 

Mo 

(13) 

ty (I) — X/ C (0* I &*■>/ Ua—€k) ]P,<£jfcri-a*j4i$o. (14) 

Mo 

Here P f is a symmetry projector, Si and T{ are opera- 
tors which differ in the various formalisms, and a { is a 
constant which fixes the normalization of Explicit 

expressions for Pf, Si, Ti, a { are given in Table I. 

The sums for \E (2) and ty (1> are easily evaluated 
numerically. Below we shall consider in detail the ana- 
lytical evaluation of »*£<*> (EL-HAV) for large P. 
Having obtained ty (1) we construct the trial wavefunc- 
tion analogous to (9), 

^=V' (0) +XfV m , (15) 

and compute the expectation value of II, 

{ B- J ZT$)/ | $), (16) 

with Xf optimized. 


16 W. M. Frank, J. Math. Phys. 8, 1121 (1967). . 
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A desirable feature to require of any symmetrized 
expansion is that it produces an expression for the 
energy which agrees, order by order, with that given - 
by the RS scheme for R large enough so that terms 
0(<r fi2 ) may be neglected. The results in Table II ’ 
show that the EL-HAV expansion fails to meet this . 
requirement. 

For R large we may write 

<£®(EI^HAV)~*£LM^ 

fc-1 €Q— 

(17) 

Jt— X 

where the tilde indicates asymptotic equality. The first 
summation is |£ (2) (RS). To evaluate the second, we 
use the completeness of the states 4>l to write 


=££ <*>1 (v-‘e»>)<i>,) 

O m-0 

<0m | (is) 

x — i €o— e* 


The advantage of this procedure is that (<£ 0 | ( V— *E (1) ) ) 

vanishes for /> 4. The overlap integrals (</>i | Ifo)-are 
given by the formula 




.l+X—2a"| 

/!/<:! J 


1/2 

e~ xl2 , (19) 


where min(Z, k ) is the smaller of l, k , and x=2R 2 . 
Using this expression and to— €k— — k , the summation • 
over k in (18) may be evaluated exactly, yielding 


V (», I /») <»!/».) „ 

fc-1 co“ci Vitfil/ 


•X Ei(tf) ~ 

L * 

1 * ” y£ -. 1 

a(l—a) \x? / a-\ ot(m— a) b^J’ 


( 20 ) 


where Ei(s) is the exponential integral, 16 y =0.57721* • • 
is Euler's constant, and M is the larger of l, m. The 
Ei(») has the asymptotic expansion 16 

*’ <r*Ei(s)~I2 ^ (21) 


'so that 

V' (<i>l j Ifo) {<f>K 1 


w *-* W/ 

The nonvanishing matrix elements of V in (18) are 
3 


I 


x i 3 ^ 1/2 » 

32 R 2 610 8R 5<1+ 16R 2 5,2 

6 1/2 

(23) 


where 8i m is the Kronecker delta. This leads to 


g. 1 l 

x- 1 



. (24) 


Thus, (EL-HAV) has the asymptotic expansion 


‘E m (Elr-HAV) (11/64E 2 ) - ( 141/1024.R*) + • • • , 

(25) 

so that for R sufficiently large so that terms OCR" 4 ) 
may be neglected, *£ (2) (EL-HAV) is one-half of 
£«(R S). 

Finally, it is noted that none of the perturbation 
schemes is strikingly better than the other three. In 
fact these results indicate that perturbation theory 
may not be a very efficient way of computing the total 
energy of systems of this type. 
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4.3 The Delta-Function Model . 

In this Section the different formalisms discussed in Section 4.1 

are solved to infinite order in the perturbation for the double minimum 
49 

delta- function model of the hydrogen molecule ion. The exact 

u 

Schrodinger equation and the polarization expansion have been con- 
49 

sidered previously. This model is of interest because the exact 
energy can be obtained in all the methods in analytic form. Further, 
all of the excited eigenfunctions of the “separated atom" problem are 
continuum functions, so. that to .the extent that the model is analogous 
to molecular problems, the role of the continuum in the evaluation 
of second order energies may be clarified. 

The model consists of the one-dimensional > Schrodinger equation 


( H -*E )*!$ = 0 5 (4.1) 

where 

H = - i + U (x) , 

and 

IT60 = - S(x- f) - $(x v %.) ? 

where §00 is Dirac’s delta function. The model has two bound-state 
solutions (“E , 1 5 ?) , where - the superscripts indicate symmetry 
with- respect to inversion 

As increases without limit, and' each tend to -% . 

This is analogous to the behavior of the lowest states of . The 


model differs from hkj* 


however, in that there are no algebraically 
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decreasing terms in ^fE(R) and, hence, the ungerade state has no 
van der Waals minimum. 

The following definitions are, introduced: 


H‘*’ - -tfc - S(xt-%) 

6 <e * = - t : ^ = 


- ix+ %. 


* - S(x- */*) . 

' J 

At - i [ i ± (x-+-x)] 0 


Further, the Sternheimer hamiltonian^ 7 for this model is 

1 H (V •= ~M (0) * “ -» <£* *■ -jr^=r U(x} 

To solve the various perturbation formalisms to infinite order, 
it is convenient to introduce an interpolation equation for each 
formalism. These equations .are - listed in Table. 4.1. Each contains 
a parameter A in such a way that for A"^ , the interpolation 
equation reduces to the unperturbed problem; and for X“ j , the 
equation is equivalent to the exact Eq. (4.i). The solution of the 
interpolation equation is expanded in powers of X to obtain the 
perturbation energies and wavefunctions . 
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Exact Interpolation Energies . The solutions ^ to the interpola- 
tion equations in the regions \ x| are given in Table 4.2 in 

terms of unknown constants. These constants are fixed by the con- 
tinuity ^ ( t % -S- £ ) " <1? ( * R/z. - £) 5 J , normalization, 

and cusp conditions at .the singularities. The last are derived by 
integrating the interpolating equations over the interval 

£ £ % -fr 6 ^ £> , and are listed in Table 4.3. 

In all cases, there are precisely enough equations to fix the unknown 
constants, including E , uniquely. 

The result of the above procedure is a transcendental equation 
for £ , as given in Table -4.4. 

It is noted that the Sternheimer and Eisenschitz-London- 
van der Avoird procedures yield identical interpolating energies. 

This is apparently a peculiarity of the delta-function model, and 
not a general result. 

Given the expression for the exact IE , the radius of convergence 

of an expansion in powers of °X can in principle be determined. 

49 50 

Although this has been done for the polarization procedure, 5 the 
complexity of the equations discourages the extension to other proce- 
dures. At any rate, from a practical point of view, interest centers 
on how well each of the IE , accurate through second-order, mimics 
the exact energy. 

Second-order Energies . Iteration of the transcendental equation 
for exact energy yields the expression 


*E « - 


Z 


+■ 


e - R * e~ 2R 

1 *■( i+po e~ R 


te-* R - 


o 
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The first term is the zeroth order energy. The second is the Heitler- 
London first order energy. The third term will be taken to be the 
correct asymptotic behavior of second order interpolating .energies . 

Expansion of the transcendental equations for the interpolating 
energies in powers of X leads to the results in Table 4.5 for. the 
asymptotic form of the second order energies. The algebra involved 
in these -calculations was checked by solving the firsts order equations 
directly, as discussed in Appendix E. 

Discussion . It is noted that the second order energy in the 
Eisenschltz-London-van der Avoird method differs from that in the 
polarization expansion by a factor of two, as has been found for the 
models previously discussed. For the delta- function model, however, 
this difference is necessary * to correct the wrong behavior of the 
second’ order polarization energy. The correct asymptotic behavior 
is also obtained by the Hirschfelder-Silbey and- Sternheimer procedures, 
but not by the Murrell-Shaw-Musher-Amos procedure. 



TABLE 4,1: Interpolation Equations for the Schrodinger Equation, (4.1) 


Equation 



5 

Polarization 


Eisenschitz-London, 9 van der 


or Hi rschf elder HAV 


Mur r till- Shaw or 


Musher-Amos 


Hir s ch f e lder~ S i lbey 


C H te> + H io - ± e ) ""S’ = o 
( H lo> + A H en - £(*) ) = O 


Avoird 15 l H to) - <? lo> *At(ZH e,> - *B ft) + £ M> 


* * 


■ 


( H <0> +; 3 - 0-X) l ’At , H°>]<f> (o> 

I At 14 i0> > = <4 ec> l A t !4 Co> > 


(tt (o> + A H. a> - *Fft>A,. - ~EMA_)±0-)=o 






< I | $ (X>)> 
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TABLE 4.4: Transcendental Equations for Exact Interpolating Energies, i Ert)= • 


Procedure 

Equation 

Exact 

s -V±R. 

>4 = lie 

Polarization 

- 2 K R ! 

- «) = 1 e 

EL-HAV 

0 t 1 ± a" * ±R ) 

V± ' It e- R 

MS-MA 

— R r - 2. R y, 

(» t -A)(»,-0 . * ^rs [7T^- e ,-~4 _ 1 

HS 

f 

-^(i+Wti+e-VRjir V- - i(i+iKi-£- 5r - R )l _ , , 

I + e.-^ R JL I - e - *- 1 * J" ? ° ; 

- O 1 i~e~ R ’T^~~0 o ~ e - ^ 1 

lve'^ R J \ve- R - L i - e -y " R J 

RS 

(it i ± e ***) 

| t gr-fc 


P20I 







TABLE 4.5: Asymptotic R-Dependence of Second Order Energies. 


Procedure 


Exact* 

, 

- 1 
2. e 

+* 

Ofe" 3 *) 

Polarization 

- e" 2R 

-t. 

Ofe -3 *) 

EL-HAV 

- J- e -2R 

+ 


MS-MA 

x e~ 3R 

+- 


HS 

- 

-v- 

0(er 3fi ) 

RS 

^ -!• e- 2 - 2 - 

>- A c. 

-f- 

O(er^) 
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APPENDIX A. Notation for Chapters Two and Three. 

The compact matrix notation used in Chapters Two and Three is 
called by Lowdin (see, for example, ref. 62), the “fat notation",, 

i)' Matrices are denoted by fat* symbols; e^g. in Eq. (2.13), 



ii) The sum and product of two matrices is defined in the usual 
‘way. The adjoint of a matrix and the multiplication of a matrix by a 

scalar is also defined in the ord-fnarv wav. 

iii) The- matrix' of an operator (£) in the space ‘is denoted 

by <$ I (£> I <§.> , where (<<£/© !<$>) w = 

iv) If is a set of nonorthogonal basis functions, the pro- 
jector onto the linear- manifold spanned by ^ 'is 




104 


APPENDIX B. Algebraic Manipulations of Perturbation Equations. 

B . 1 . Calculation of ^ CN) [ H E. I ^ * 

According to the perturbation equations (2.22) 

U <0) + AH (0 )4 fc) s +4^€ c % 

A H t0 <4 te) - £ U t0 4 (!) + 

* 

to 

(h“>.xh" > ))"4‘ w + >'h“ , 4‘ m> = x " + ‘ h‘ ,, 4"' :i 

+ 

A~o 

Adding the’se equations gives 

W$(n> - +£x k z 

k-° A-O 

Exchanging the order of summation over and Jl gives 

H4(u) = a m+1 H t0 4 m) + Z x k <£ (fe) 

JL^o 

= <^CN)& CM) 4 > ,i+l H t0 4 W 

-ZjV'2 >*<*“" 

i=/ A =*-J+i 
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Collecting terms according to order in ^ , 

wioo -$CM)CCM) + / +, i>V w) -Z X fe Z4 (6 ^H 

1 O 4=6 


Hence , 


<«$(n) I H-E. = 


= </z£CA) ( <SfN)— E. ) 


4 -> m+, { i> X h < 4 c *'MB to | 4 fN) > 

1 W=fe 


M N _i nMi-i 

s Z V* l <4‘*l4 tfc+JM) 

k'6 J? 5 ^ 
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Again collecting terms according to order in X yields 


(4M I ® l<Soo^(<£oy-E.) 


+} M+ ‘Zr f<4 tn) Ih w I 4 w > 


n ^ + f-* 

-2 2 . < 4 l ^i 4 trt+ ‘"‘H > ><s 0 "^ 

a n s. /-* ^ 


Z^f <4 CHj iH t0 l4 tN> >s* o 

h"o 

/s 2 <4 c r*»' w |4“-tnJ>« w i? 

f j- 


B .2 Perturbation Energy Formulas . 

In this Section, the Eq. (2.32) and (2.33), resulting from the 
full normalization condition Eq. (2.30) are proved. Also it is shown 
that (& ^ is a hermitean matrix. 

First, however, it is convenient to establish some intermediate 
results. By taking the adjoint of Eq. (2.30) and adding it to the 
original equation, one obtains 

£ < 4^1 4 c "- fe) > = o, 

k~o 


y \>6 


(B.l) 
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The following } relation will be needed- later : 


< 4 a '° | H 6 >i 4 6 "' 30 > - <4 UJ | H ( ° 1 4 1 > 

(B.2) 


■2 < 4 ^ 1 4^- £ '^>(g' fe) * 2 

fe~0 


%This is proved by ; multiplying, the order perturbation equation 

(2.22) from the left by <2^^ ^ and integrating, 

W ,n -* ) \ h'^ I 4° V > 1 + <4 tv1 ' £) i'H t0 14 <*-'>> 

k (B - 3) 

L <4 cv "- 4) 1 4 «-*»>> <£ c « 

k-o 


1 
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and multiplying the (n -j£ )-th order equation from the left 


by 


and integrating «, 


+ <4 a > | h Li) j 4 O'4-o^ 

(B.4) 

= z <4 U) l4 l ”-^ h) >(£ Ch) 

o 


Taking the adjoint of Eq. (B.3) and subtracting it from Eq. (B.4) 
yields Eq. (B.2). 

From Eq. (2.29) 




Application of Eq. (B.2) ^ times yields 

* < 4^1 4 i 


-X ^-<4 d) l4 ct ' >fe) >£ (t) - ^ <4 Crt i4 c ' t ' b> ><£ ( ! ) 

£-i k-o 
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Interchanging the order of summation gives 

% * * k 

- < 4 c i>i H to 1 4 c t'S-')> + 2 Z <4 c?) f 4 

A- & 


The final term vanishes by the normalization condition; Eq*. *(B,zL).. 
Hence j 


% <g-k 

6 tf) = <4^iH < °)4^r > >^^^ + 2<4 W) i4^ !l ' Ji) > 




J3=o 


*1' h (■4 <1> I4 ;+ ^' i ‘ , >« < ’ b . 


(B.5) 


The desired proofs are given by the following applications of 
Eq. (B.5): 

Put p=n, q=n- 1 : 

<L (tii - <4 U " ) I H ( 'H £0) > ~ Z <4 m, 14 Cv ''^>« Co> 

JL~* 


*’£*«* 2 4 U M4 . 



uo 


The coefficient of * is simplified by using:Eq. (2.30): 

a M ' - < 4 (n '° | h 10 1 4 <e5 > ■+ <4 te> |4 w > ^ 

§ <£<*> + < 4 C *~‘° 14 C0) > 

t*o- ‘ ? 

or 


Put p=2n ? q=n-l : 

& l "» - <'<4 C "~° I H°* I 4 t ” ) > 
+ 2 ! <s Cl ° 2 * <4 tA) l 


4^ 


-2 

k-O 


Z<4 w !4 tzh ^' k) > & (k) . 

J3=o 


(B.6) 



Ill 


gCzrt) _ H c.) I^evo^ 

" Z Z <4 C ‘'' fe+J2 ) ; | 4 Cw '^ ? 

k=l Jl-t 


rt-l H-k-l 

+ Z Z o I (E Cti) <4 C£) |4 C2v, - k - fi >> 

- <<4 W) ( 4 C2 ^- •«•--« yg Cv) j, 

" i 1 ^ Co> <4 Ch> 1 4 c>,) > - \ <4 ^ co) 5 


In the one -dimensional case, when .all of the .quantities in this 
equation are scalars,, the Eq. (2.6) results. In general, the 
equation is simplified by using <£ ~ Cz ' v ^ 4 ^Czvoi" ^ • 



112 


= i f <<4 lw ' l) l h < °\ 4 Cv ' > > + <4 w ia a M4 (H "'>] 



(B.7) 

+a£ X d tfe) {<4 W) i c 4 fZh '^>-<4 c ^- A >|c4a)^l 

h*l 1-6 1 


n-t VHa-i 

-■j x X {<4 w t4' 2,, ' ll ‘ i) > _ <4 c ^- fe -^)(4 «)>{ £ Ck ] 

k-\ $'o 1 


From Eq. (2.30) 


h~U-l 

X <4‘»l4 tew - 


r 


fe-Z 


,(-Ju( - )y 








m 

I st ° 




1 4 W " £) > 

fecdd < 


(B.8) 
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Substituting Eq'./ (B. 8) into Eq. (B.7) yields Eq. (2.33)-. 
Put p-2n-l, q=n : 

By steps similar to those used to derive Eq. (2.33), one 
Eq. (2.32). 


obtains 
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APPENDIX C. Partitioning Technique Approach no- Degenerate Almost 
Degenerate, and Exchange Perturbation Theory, 

It is the purpose of this Appendix to clarify the connection 
between the formalism of Chapter Two and the partitioning technique. 
Only the degenerate case will be treated in detail, 

36 

The basic notion of partitioning goes back to Van Vleck, 

Lennard-’Jones , ^ Brillouin,^ and others in the early days of 
n 65 

quantum mechanics . Lowdin and his coworkers have developed the 
partitioning techniques into a very elegant apd powerful approach to 
perturbation theory. 

It is convenient to .first quote some results of Lowdin, ^ Let 
and E. be an eigenfunction and eigenvalue of the hamiltonian H , 


(H - E) 5 = 


O 


(C.l) 


The Hilbert space of M is divided iptq two subspaces. The fir si: is 


spanned by che set of functions and is 

characterized by the projection operator ( 5 > ■ The second subspace 
is spanned by an infinite set of functions and is characterized by 
the projector T = 1- C8. These two operators satisfy 


‘j p^= p = p + • o 1 ? - Pc? — o. 


In the partitioning technique, the Eq. (C.l) is replaced by 

0 C H - B ) (9 » o f 


(C.2) 
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where 

S = h + w tce; h. 4 

and 

V(b) * PE*£* P(e-h.)P.T’ p, 

with c£ t^O but otherwise arbitrary. Eq. (C.2) is equivalent to the 
secular equation 

!<4 <0 ’ I H - E i| 4 (0) >ho. (c4) 

-This is- not an ordinary secular equation pince H is a function of 
& • The usual approach is to solve Eq. (C.3) by an iterative 

process. For the degenerate case, writing 

H =• 14 (0> + > H 

where 

allows V( E) to be written 

Vie) = ?[o<?0+ P(e'°EH lcl) )? - ?-(AH“ > +e' fl -E)p]] P, 

Now defining 

+ P P , 
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and using the "fundamental identity" 

(A-B) = -+ A"' 13 ( A - B)“ 1 

yields 

r CE) = R 6 + R 0 ()iH rt + e tc) -E)T(E), i (c - 4) 

This relation allows ^(E) to be approximated snqcessively . 

It is possible, as Lowdin^ in the nondegenerate case and as 
Choi^ and Goscinski and Lowdin^ in the degenerate case have shown, 
to obtain the Rayliegh-Schrodlnger expansion of the energy and the 
wavefunction by iterating Eq. (C.4) and expanding the roots of the 
secular equation (G.3) in. powers of X 

Of course^ the expansion (C.4) can have a greater radius of con- 
vergence than the expansion in powers of X of the roots E of the 
secular equation (C.3). Furthermore, given T P(E-) the secular equation 
can be solved by techniques which avoid power series expansions. 

Intermediate Normalization . In applications of the partitioning 
technique, the natural choice of normalization is "intermediate” 
normalization , 

<4 00 l < 4) > - = i 5 (c.5) 

rather than the full normalization chosen in Chapter Two. Thus it is 
necessary to derive the expressions for, and which are 

applicable for the normalization Eq. (C.5). 
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To avoid confusion with the results using, full norma liz at ion. 3 
subscripts rather than superscripts will be used to denote order in 
intermediate normalization . It is also convenient to absorb X into 
the definitions of the perturbation terms and to define 


V * 1 H° j 

V' * V- (b ~ £ e ) } 

K-\ ’ k " 

l r> k ^ I 

Expanding Eq. (C.5) in powers of X yields 


<4c I <4 h > * O , n » I . 


(C.6) 


Multiplying Eq. (2.22) from the left by 



and integrating gives 


- Wo 1 vI 4 h -, >. 


Eq. (C.6) allows to be omitted from Eq. (2.34) so than the 

perturbation wavefunctions are 


V' 


- 2 R* 4, 


I 
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Iterating this equation s it follows that 


* (RoV^r 1 ^ "Z E (RpV^R 


/C.8) 


4=0 b=i 

Expansion of . Iteration of Eq. (C.4)' (2N-1) times yields 


Tie) = . IE ) 4- (R 0 VO^Tte) (C.9) 


where 


%„.,<*> • &> 1 


(C.10) 


With the neglect of the second term in Eq. (C.9) the secular equation 
(C.3) is 


IHI + V + V T aw „ ( CE) V -E \ 4 o>\ - (c - n) 

It can now be shown that the roots of Eq. (C.ll) are identical, 
through terms O ) , to the roots of the secular equation 

(2.26) . 

This theorem is proved by using Eq. (C.8) and (C.10) ana tne 
fact that R 0 <4^ = © , to obtain 


AN~I 


, n-t n-£ 

Vf 10 H - £ + Z 2 Z 

i 1 


6rj=-o 
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Substitution of this result into Eq. (C.ll.) and the .use of Eq. (C.7) 
gives 

o = I I H 0 +v +■ v v - E 14> >1 




z <i 




E 


+ I^I £'S| ^ +l ^w+I- 4 -fe.^ 


Mfl n-} B-tf 


*£ 5 2 i-H i & So) 1, 1 


earranging the order of summation in the last two terms gives 


2 - E 

h-o 


0= 




>Z € 

te= 1 


fe 


aw-*: a.ta-n-2. 


2< Z Z Z<HwtRo(v' 


rv=i 


t 


i ss o 11=6 




k=i 
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It is now clear that, in the summations over fe, the upper limit 
may he replaced by if terms are neglected. 

Hence, 


1 


n 


0 = 


+ Zl 





|>-o 


£ 


|<t,jR.(v'R^I4 t . i+( . T > 


I 


su>s ■m 

Z <£ n - E ^ 0 (x ztl+<l ) 

h-o 

Assuming that the first determinant is nonvanishing yields 

u 


(G. 12) 


Hence, the energies determined by Eq, (C.ll) are identical with the 
energies determined by Eq. (C.12) through O ) . Further- 

more, it was shown in Chapter Two that the change from intermediate 
to full normalization also affects the roots of Eq. (2.26) by 
O C> a -" +a ') terms. Hence, it follows that the roots of secular 
equations (C.ll) and (2.26) differ by OU^) terms . 
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Received 21 October 1968 


A transformation of the unperturbed Hamiltonian in the Eisenschitz-London, Van der Avoird, or Hirsch- 
felder HAV exchange perturbation theories is suggested, which modifies the second order energy to produce 
the correct long-range behavior. 


1. INTRODUCTION 


T = R 0 +R 0 V'T , 


T =V-{E-E 0 ), 


In recent applications [1,2,9] of the Eisenschitz- 
London [3], Van der Avoird [4], or Hirschfelder 
HAV [5] perturbation theories for the exchange 
forces between atoms, it has become apparent 
that the second order energy in these schemes 
fails at large interatomic separations R to equal 
the second order dispersion energy (no exchange), 
even though "exchange effects" are negligible at 
these separations. In the previous calculations 
[1,2], the correct long-range limit was obtained 
by a variation-perturbation approach. In this note 
we suggest an alternative method for obtaining 
the correct limit, by introducing a transforma- 
tion of the unperturbed Hamiltonian of a type ori- 
ginally suggested by Feenberg [6]. 


2. PERTURBATION THEORIES 

In this section we review/ briefly the usual 
Rayleigh-Schrbdinger (unsymmetrized) theory 
and the Eisenschitz-London, Van der Avoird or 
Hirschfelder HAV theory for the long-range in- 
teractions between atoms. The problem may be 
stated succinctly by making use of the wave op- 
erator formalism [4,7]. Conventionally long- 
range interactions are computed from 

E =E 0 + <0 O |F+ VT V|0o> , (1) 

where 


Ro = P(B o - Ho)' 1 ? , H 0 <b 0 = Eotbo , * 2) 


^ = 1 -IVfyjl > 


where H Q is a sum of isolated atom Hamiltonians. 
0 O is product of isolated atom wavefunctions, and 
V contains the interatomic interactions. Substitu- 
tion of eqs. (2) into (1) and iteration yields 

E = E 0 + ei + €2 + ^3 + • • • . 

where 

e 1 = <$o I V\ <b Q ) , 

€ 2 = (<&ol VR q V\ <p 0 ) , 

€3 = (<t>o\VR 0 {y-ei)R o V\<b 0 ) . 

When exchange effects are not negligible, the 
above treatment is not sufficient because <b 0 does 
not have the symmetry of H = Ho + V. Let A be 
the projection operator onto the subspace of de- 
sired symmetry. Then the interaction energy 
may be computed from 


E - E o + 


(<b 0 \VA + VT V\ (b Q ) 

(d>o\A\<f> Q ) 


(3) 


where 

T = PR 0 P + PR 0 UT (4) 

and the former definition of P is replaced by 


* 7 his work was supported by Iho National Aoronaullcs 
and Space Administration. Grant NGL-50-002-001. 

* National Science Foundation Graduate Fellow. and 


A 1 0o? <<fr 0 1 A 
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Table 1. 

Energy of H 2 {in cm~l). Here c = ((p 0 \ A ! ^o)** 1 and X is the value of the parameter in eq. (5) which optimizes the ex- 
pectation value of the Hamiltonian in refs. {!] and (9J. 



R 

X 

c 

Eq +■ Ei + R% - B {exact} 

Eq+Ei* XEz - S (exact) 

Eq + S 1 + cE 2 - & (exact) 

E (exact) 

Ground 

4 

1.36 

1.93 

427.0' 

177.8 

-214.2 

3592.6 

State 

V 

g 

6 

1.39 

2.00 

25.5 

9.4 

- 16.1 

378.9 

8 

1.92 

2.00 

3.8 

0.1 

- 0.2 

11.7 

Excited 

4 

1.42 

2.07 

169.6 

39.6 

-158.0 

- 1453.3 

State 

V 

u 

6 

1,89 

2.00 

23.5 

4.4 

2.2 

- 41.2 

8 

1.98 

2,00 

3.8 

0.3 

0.3 

4.3 


lame 

Energy of the One-Dimensional Model. Here c = = 2/(1 ± exp {-it 2 }} and X is the value of the parameter 

in eq. (5) which optimizes the expectation value of the Hamiltonian in ref. [2]. 



R 

X 

c 

Eq+Ei+S2~ E (exact) Eo+JSi+X-Eg-i? (exact) Eq 

+ + cE2 - E (exact) 

E (exact) 

Gerade 

0.8 

0.3302 

1.3095 

-0,1551 

0.0180 

- 0.2351 

0.3308 

Ground 

State 

1.0 

0.4528 

1.4621 

-0.0857 

0.0149 

- 0.1706 

0.2940 

1.2 

0.5964 

1.6169 

- 0.0469 

0.0161 

- 0.1432 

0.2806 


2.0 

1.3935 

1.9640 

0.0653 

0.0355 

- 0.0062 

0,3502 


2.6 

1.6236 

1.9977 

0,0473 

0,0272 

0.0151 

0.4307 


3.0 

1.6077 

1.9998 

0.0279 

0.0147 

0.0063 

0.4604 


5.0 

1.8314 

2.0000 

0.0072 

0.0014 

0,0002 

0.4895 

Unger ade 

0.8 

0.3457 

4.2310 

- 0.2404 

0.0426 

-1.6381 

1.1403 

“Excited 

State 

1.0 

0.4523 

3.1641 

- 0.0819 

0.0171 

- 0.374D 

0.9314 

1,2 

0.5536 

2.6209 

- 0.0343 

0.0088 

-0.1908 

0.7870 


2.0 

0.9494 

2.0373 

0.0025 

0.0037 

-0.0316 

0,5232 


2.6 

1.3268 

2.0024 

0,0139 

0.0056 

-0.0119 

0.4776 


3.0 

1.5031 

2.0002 

0.0164 

0.0062 

-0.0040 

0.4740 


5.0 

1.8315 

2.0000 

0.0072 

0.0014 

0.0002 

0,4895 
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»-( 1 


■*I^Q><»oU ... 


Substitution of eq, (4) into (3) and iteration yields 


E = Eq + Ei + .Eg * ^3 * * * * > 
where 

£ 1 = <0o|K4|0o>/<<#'oUl0o>. 

£2 - <<h I VPR a PV\ <fc,>/fo 0 U 1 <fo) . 

£3 - <<*>o I W 3 B 0 P(r- £i)££ 0 £^I *o>/foo U I 0 o> • 
The difficulty is that, for the systems consid- 
ered, 

E 2 

lim j — *=£ 1 , 

a— 2 


so that E a Eq + Ej + E2 is a bad approximation 
for R large. 


3. PERTURBATION-VARIATION APPROACH 

A better approximation was obtained [1,2] by 
computing the expectation value of H with the 
variational wave function 

f-^o + XPKoP^, 1 {5) 

where the variational parameter A was energy 
optimized, with the result . „ v 

E - -Eo + B\ + XE 2 , 

and 

lim \E 2 /e 2 = 1 • 

/t-**«* 

This result has been emphasized recently by 
Goscinski and Br&ndas [8], 


4. THE FEENBERG TRANSFORMATION [G] 

The computation of A requires a knowledge of 
E 2 and Eg [8], A procedure for removing the 
long-range defect with a knowledge of Eg only is 
as follows. Instead of iterating eq. (4), one iter- 
ates 

T a cPRqP + cPR 0 UT + <1 - c)T , (6) 

which is strictly equivalent to eq. (4) if c is an 
arbitrary constant, and obtains through second 
order 


E a Eq + Ej + cEg . 

This method is equivalent to replacing [6 ] (//q - Eq) 


by (l/c) (Hq-Eq) and Fby F-{l-c/c)(E 0 -Eo) ™ 
eq. (4). 

In order to obtain the correct long range be- 
havior, the constant c must satisfy 

lim cRgAg = 1 . CO 

For the systems previously considered, the 
choice 

c = (^oUlc&o)- 1 (8) 

satisfies eq. (7). In tables 1 and 2 we present a 
comparison between the variational approach and 
Feenberg transformation. 

For the systems considered, the choice of 
eq. (8) is consistent if the energy is computed 
through third order. Iteration of eq. (6) twice 
yields 

E =Eo + e(2~c)E 2 +c%3 . 

In these examples, c ~ 2 and E3 -* |E 2 as R — «*, 
so E -* Eq + 2E 2 which is equal to Eq + € 1 + e 2 for 
large R . 

In general, however, the computation of E3 
requires a large increase in the number and com- 
plexity of integrals over those needed to compute 
E 2 - For this reason the Feenberg transformation 
offers a practical advantage over the variational 
approach. The significance of the choice of c 
which we have made escapes us, however. 
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APPENDIX E. Direct Calculation of Second Urder Energies tor the 
Delta-Function Model. 

As was shown in Section 4.1, the second' order energies may be 
evaluated by a quadrature involving the function defined by 


(H m - € m ) + A±(H C0 - *E (E.i) 

| 4> t0> )> — 0. (E.2) 

The second order energies are given by 


. < t KlA±a\ l »-*E c ")\ 4 i6} > 

-E (EU-Hftv) ~ << 4 » £rt l A±\ 4 W > 7 




<*x \ hh'M4> (0) > 
< 4 > < «IA t |^“ , > 3 


± w <^\R t, H4 ce> >-ihE t0 -"£ w X :fc ^ttl4 ta > (e.5) 

~ E (HS) “ <4‘ e) l a ± 14 w >> 

Introducing the definitions of EE 

4**° , and A± , Eq. (*E. 1) may be written 


[£, * ncx+a) - (fX + [ S(x- a)4 ±E co ] e 


t [S(x+«.) + ± E t0 ]e 


(0 n -Ix-al 


where 


a = n h- } 
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or, omitting the (+) sub- and superscripts, 


-'lx + E“e- lx ‘“' t ui.a: 


with the matching conditions 


" ^ X L > 


" i X La- +e R * 


XL- a-0 - XC-a-) = o • 


X Cat) - X La-) — o . 


The general solution -of Fq\ (E.6) 


'W* -A 1 ** 1 " xe" x 


y- 


B°V + C { ‘*e - * 

|_+ xe* - c te> ) ft, xe x j 


[ £> to e K * 'D“ 0 Y <0 xe* 


X >«. 


; |xl<a 


(E.6) 


(E.7) 


*(E.v8) 


CE.9) 


* ( E.10) 


(E. 11) 
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yrtiere 


A m = i (e~ a ± e a ) 5 

B to> = ± -k e -a ; 

C fe) = ± e" a „ 

U«> = -k ( e a ± s“ a ) , 

C&5 




(I) ^ 


ie ^ e~ 2f ^ 


i +• ( i > 

and A 1 ' 1 , are chosen to satisfy Eq. (E.7) 

through (E.10). Thus, 


e ft e-°'\ / a ci) \ 


e’ a -e a e~ a o 


o 

-a 


e~ a e a 


~a «■<* 

*e ' -e o 




iU) 


s (») 


_ \ 

/•Y c,> [i+a* a,e wK ]+ e ^ 


(u 

±aY“> 


(E. 12) 


These equations are not linearly independent, however, as may he 
verified by comparing the sum of the first and third equations with 
the sum of the second and fourth equations, and making use of the 
definition of . Thus, one of the Eqs.(E.12) must be replaced 

by another which is linearly independent. This is just the freedom 
required to enable the normalization condition, Eq. (E.2), to be 
satisfied. Making use of Eq. (E.ll), Eq. (E.2) becomes 
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e R_a A 1 " 4 2 .Re" ft E° 5 4 e" a (e R -£-^) c“ } . 

(E.13) 

+e a-R - » i [_ I £ C 140 e -*] V‘° . 

Eq. (E.13) now replaces the first equation in Eqs. (E.12). Then, 

> 8 tl> , C C, \ can be obtained in a straightforward 

manner. The result is 

A U) - iv 6) e a [+ ( 2 --R.') + Ci-R.')e"" E> t (3+3R.+ R.' J -)e“ 2E: -)] 

4 4 e”** 0, [ I - (i+R.)eT 2R 3 

B <0 * :j:Y w e”* C+ (z.+sql + ie' e ' a j 

C co = ^^ c,) e' a It-Rt (?+SR-vK.' z )e~ B -l “ i C i+b.) eT 2 * - * 

D t0 = ^ * w e* jj+ R ±(i+R) 2 e -E l ~ t- R e ~ R 

This completes the determination of ^ 7C . 

In order to evaluate the second order energies, the following 
integrals are needed: 

<*X|H‘ 0 I4>.> = y'°[-| e _2R + $(3+?&+3R*)e~ ZR ] 

± ie~ Zlz ii ± (i + K.)e _R l 5 
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< i X|XH w |^‘ 0, > = -^ £ V R [li ( I+R.+ ft z )e -R- l +i Re* 

< ± X|H4 £a) > ^r !T u) [ : i :3 -^e'* a +(i+e)(3+S!i+R i )e“ 2s J 

•4-t e - ' 4 ' [l - 0+ R) 2 -e -ZR -3 , 

With these integrals and Eq. (E.3) through (E.5), the asymptotic 
dependencies given in Table 4.5 are obtained. 
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APPENDIX F. Other Publications of Phillip R. Certain 

Reprinted in this Appendix are the follpwing articles, which 
report work completed at" the Theoretical Chemistry Institute, 

Madison, Wisconsin, 

1. M A Single Variable Variational Approximation for a Perturbed 
Schrod-inger Equation" W. Byers Brown and Phillip R. Certain, 

Proc. Nat.' Acad. Sci. 57, 1206 (1967). 

2. 11 Comment on Van der Avoird‘$ Wave Operafqr Formalism", 

R. E. Johnson and 'P.. R. Certain, Chem. Phys. Letters, JL, 413 (1967). 

3. "Calculation of Matrix Elements for One-Dimensional Quantum- 
Mechanical Problems", A. S. Dickinson and P. R. Certain, J. Chem. 
Phys-., 49, 4209 (1968). 

4. "Localized Orbitals for Arbitrary Molecular Wave Functions", 
P. R. Certain and.J.,0. Hirschfelder , Chem. Phys. Letters , 2, 274 
(1968). 
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friiwi Uk« VnWKKMSti* <» vnr Nahonai, AvMtKm nr StithschU 
V«u 57, No, 5, !>{>. 1205-1212, May, 1007. 

A SINGLE VARIABLE VARIATIONAL APPROXIMATION 
FOR A PERTURBED SCHRQDINGBR EQUATION* 

By W. Byers Brown and Phillip R. Certain! 

THEORETICAL CHEMISTRY INSTITUTE, UNIVERSITY OP WISCONSIN, MADISON 
Communicated by Joseph 0 . Ilirschfeldcr, March 2/ h 1067 

1 he difficulty in applying the lUylcigh-Ritss variations] piefchod to molecular 
systems lies in choosing the trial wave function, £ In this paper we consider a 
variational approximation for the ground state of the perturbed Schrodingor 
equation 


(I7o + V — E)yf/ = o, (j) 

whore the perturbation V is a function of the configuration variables, and the 
ground state eigenfunction of Ih is known. We write the trial wave function 

0 * (2) 
but instead of guessing a specific form for F, we seek the best F which is a function of 
V only, F — F(V). Tins leads to a Schrodinger-type equation' 

(^o + V — E)F =s 0 Qfy 

in the single variable V, involving the approximate ground state eigenvalue E(^ E). 
Equation (3) is an ordinary linear differential equation which can be readily solved 
numerically. 

The analogous first-order variation-perturbation theory has been derived re- 
cently by Ivirfcman and Benston,* and called by them the extended average energy 
method. They present their derivation as an extension of the Unsold or average 
energy approximation, 2 * 3 which is equivalent to employing a first-order trial wave 
function of the form * - ft with / - — (F - E X )/AE, where Ag is an average 
excitation energy. The functional minimization method used in this paper leads 
to essentially the same results as those of Kirfcman and Benston. 1 

Perturbed Schrodinger Equation . — Consider a system described by the n orthog- 
onal curvilinear coordinates q = q h q 2 , . . ,q n . Let the unperturbed Hamiltonian be 


7/oa-V^i V 0i 

where U 0 (q) is a function of the coordinates, and 


V 2 « 2 - — 
*«i hbq. 


ww, 


where the fc,(q) are the metric scale factors, and h = h x k 2 , . .h n . Let ^o(q) and Eo 
be the lowest unperturbed eigenfunction (assumed real) and eigenvalue satisfying 


(Hn ” Ef»ft — 0. (,j) 

Wo arc interested in the solution ty.E) of the perturbed Schrfldingcr oquntion «) 
m the case m which the perturbation V is a function (not a differential operator) of 
the coordinates q; or more conveniently, in which ■ 


1206 
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V = V'{u) where u = w(q) with- range (a,b). (5) 

The variable*?* may or may not be one of the coordinates q. 

Varialionally Reduced . Schrddinger Equation . — Wo sock the best trial wave 
function of the form 


$ = ^(w)^o(q) . ((>) 

by means of the Baylcigh-Ritz variation principle SE = 0 with respect to arbitrary 
variations SF(u), where the functional E[F] is defined by 

(TIo + V - m> = 0. (7) 

By substituting for \p from (6), using the unperturbed Schrodinger equation (4) 
and rearranging, equation (7) may be rewritten 

{WF) 2 ± (V + E. - = 0. (S) 

Since?’ = F(u), (VF)' ! = (^j' (Yu)*, where (Yu)* = E -t Wo define 

two functions. of u, which involve only the unperturbed distribution: 


ifo(«0 = (f 0 , 5(u' — u)f 0 ), (9) 

QoW = tyo, S(u' - w)(Vu)Vo), (10) 

where 8(x) is the Dirac delta function. Pq(u) is the probability density for the 
variable u in the unperturbed ground state, and is- normalized to unity provided ^ 
has unit norm) 

** J Po(u)du - (Mo) = 1. 

In many cases (Vu) 2 will be a constant, say c; then Q a = cP 0 ; but in general Qo 
differs from Pq. 

Equation (8) can now be written as the single integral 

la dw Qo (3 +<y+ E °- H - o. 

By carrying out a functional variation 8F(u) and setting bB - 0 we obtain 

• f \Tu { q " 3 + < F + yi » - » + 1 [ft g tfj - o. 

. Since the variations 5F are arbitrary, we must have 


1 d J dF\ 

~2 du\ Q *dd + (y + B '-® P * F = <>> . (ID 

and Qo(a)F'(a) = Qi>(b)F'(b) = 0. (12) 

Equation (11) can be written in the form of the variationally reduced Schrodinger 
equation (3) where tho single-variable unperturbed Hamiltonian *is defined by 


fa - 


1 

o 


/V 
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The factor function F of interest is the ground state solution of (11) satisfying 
(12) which belongs to Lz(a,b) with weight factor P 0 (u); that is, FaFHu is finite. 
Note that P 0 ^ 0, and for the unperturbed ground state Po > 0 except possibly at 
the limits a } 6. Equation (11) is an ordinary linear differential equation which can 
bo readily solved by numerical methods. In a later section we discuss the explicit 
perturbation solution. 

By integrating equation (11) over the range (a,b) and using (12), we see that the 
eigenvalue is given "by 

£*& + f U P 0 VFdu/ f’ PoFdu. 

J a <J a 

This is a variationafiy reduced form of the so-called integral Hel I man n-Fey n man 
theorem. 4 

Generalizations. — There arc two obvious generalisations. 

(a) Perturbation V : The best approximation — F(n)ipa may still bo found 
even though the perturbation V = F(q) depends on variables other than u, and also 
if it contains differential operators in the Schrodinger representation. In the more 
general case equation (11) becomes 

- 2 1 £( <3 “£) +w+ «-^ f=o ' 

where 


W(u') - (\po, V8(u — u')$ Q ), 

and may involve a differential operator with respect to u\ 

(b) Factor function F: The procedure can clearly be generalized to allow F to 
depend on more than one variable, F = F{ u), u ss u lt Uz,. . However, this 
leads in general to a nonseparable partial differential equation of the form 

where, in an obvious notation, 

Q^(u) ® - u)(Vu a 'Vup)$o}. 

Perturbation Approach . — We now turn to the first-order perturbation solution of 
equation (11). It is instructive to derive the relevant equations from the Hyllornns 
variational principle for the second-order energy, P 2 , as a function of a trial first- 
order wave function 

St = <*,(*• ~ + 2 <&,(F - EM £ E 2i (13) 

where Pi = (^o,F^ 0 ) is the first-order energy. We expand F defined by equation 

{2} as a perturbation series - F = 1 -f / + where we have put F Zt .— 1, F t = 

/, so that 

“ f$ o. 

Then by substituting (14) into (13) and rearranging, wc obtain 


(14) 
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Eh ~ ~ (^o>(V/)Vo) + 2{iAo >(7 *"“ Wo>. (15) 

We again consider a perturbation 7(tt) as described by (5), and seek the best first- 
order trial wave function of the form 

q>. a®) 

Then (15) can be reduced to 

& =/>M£) 2+2(7 -^4 

where P${v) and Qo(u) are defined as before by equations (9) and (10). By making 
the functional J 2 j /] stationary with respect to arbitrary variations Sf(u) t wc derive 
the equations 

£(*£)- «*-*>* (17) 

and Qo(a)f'(a) ~ Qo(b)f(b) — 0. 

Equation (17) has the same form as the first-order perturbation equation for a 
separable variable u (or a simple one-dimensional case), and can be integrated 
directly 6, 3 to yield 



«-» + /.''§§“• 

(18) 

and hence, 

\r.m -• 

(19) 

where 

M (u f ) ~ 2(^o> (7 — Ei)A(tt' — u)\po) } 

(20) 


- 2 J U (7 - J3i)P 9 du, 



where A(x) is the Heaviside unit-stop function. 

Equation (18) is essentially the same as that derived by Kirtman and Bcnston 1 
by their extended average-energy method; our equation is simpler and slightly 
more general in that the variable u need not be a member of the q's. This result 
establishes that the extended average-energy method docs indeed give the best 
first-order trial wave function of the form (16) that is, the one giving the lowest 
value of Mi{>Ez exact) as claimed by Kirtman and Benston. 1 

Application to the Two-Eleciron Atom.— In this section we apply the foregoing 
theory to the ground state of the two-electron atom with atomic number Z. We 
take as our single variable u — rn, the interelectron distance. The simplest 
application is to tako the perturbation V to bo the electron repulsion term l/n*. 
In this case, we have (unit of length — Z Bohr) fa * 7r~ 5 exp {— ft — r») s ~ ~~% 2 f 

Ei * - Z. It follows from (9) that 
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P„ = + + 

and since (Vr K ) 2 — 2, Q 0 = 2P 0 . Unfortunately, it does not appear easy to solve 
(11) analytically (though a numerical solution could be easily obtained). We 
therefore turn to first-order perturbation theory. The integral (18) reduces to 

/ = |^ 5 u + In (u 2 + | u + + ^3 arctan (~^p)] + const *» 

and the integral (19) is equal to 16 - —0.135337 Hartreo. The accurate value 
is 6 * 7 Ei — —0.157666 Hartree. Those results agree essentially with those of 
Kirtman and Benston . 1 * 8 

A better starting point is to take the Hartree-Foclc wave function tyuv - <£(n)<£(r 2 ) 
as the unperturbed condition. We therefore seek the best correlation factor F in 
the trial wave function 

$ — F(r 32 )>/'fl-f(n/2). (21) 

We are even less able to solve the complete equation (11) analytically, but the first- 
order perturbation theory is still tractable, and not much harder. The first-order 
wave function can be written in the form & = fa HF + iAi corr where = {/hp^i) 
+ fiirto) } and & corr = Fx(r 12 ) to- The first-order Hartree factor is known an- 
alytically , 9 

f„,(r) = - i [3 f* (1 - e-y-'dt + 3(e -zr - l)/2r 

+ 2c- 2 ' - 5r + 23/4 - 3 In 2 j' 

The Hylleraas second-order energy, defined by equation (13), can be written in the 

13 9 / 3 \ 

form 10 E z = Bf r + & corr , where 9 B% ar = ~ ^ In ( 4 ) " -0.111003 

Hartree, and 

&T" = \ + 2<^o M'O'M, (22) 


where / corr = Fi(u) and the correlation perturbation (first-order) is 10 V = vr 1 + 

( e - 2r t — 1 )/n -{- (e _2r » — l)/r 2 + e~ 2r » + c _2r * + Since V depends on ri,r 2 as 

S 

well as u ~ Tn, we have to employ generalization (a). However, the resulting 
perturbation equations are again (18) and (19), with M defined by (20) . We obtain 


M ■ 


= [I + 1 “ - S + ? “ 3 - f 2 u< ] - [I + S “ + H c_< “ 

and carrying out tho integrations we find that 

= i h { _5,i + 18 ln (“ 3 + 1 * + S) - If V3 a '' ohi " C’vr ) + 
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136 e~ 2u ~ 1 
27 u + 


272 r* 
27 Jo 


25 

, + 17 


e~ 2a: d:r> 


+ const*, 


and 16 Ev — 0.032435 Hartree. This is 70 per cent of the accurate Hartree- 

Fock second-order correlation energy -0.04GG63 Hartree. Kirtman and Benston’s 
application to the total first-order wave function yielded only 51 per cent. 

It is interesting to compare the behavior of the two approximate first-order 
factors for small u and large u: 


1 5 

Zf ** const. + o « “ “ u 2 + 0(u 3 ) 
& 4S 


5 1 

16 U 8 In U const * + 


1 217 

^/corr « const. + “ U ~ — U 2 + 0(*4 3 ) 

5 15 

— J^w+g-lnw-h const, -f 0(u _1 ). 

We note that both satisfy the cusp condition for the exact first-order wave function 
factor / = rp x / ipO) Z The first terms in the asymptotic expansions for 

large u have opposite sign, however. 

Discussion. The trial wave function (21) is not the best form of its kind. It is 
better to start with a trial function 

$ = F(n2)<p(ri)(p(r 2 ), (23) 

and to vary both the correlation” factor F and the orbital function <p simultane- 
ously. A function of this form was first considered by Baber and Hasse, 11 whose 
work was extended by Green et al , 12 To the first order this approach is equivalent 
to employing a correlation wave function with the single excitations projected out; 15 
that is, 


& corr - [f(rn) - g(ri) - y(r 2 )}fo, (24) 


where g(xi) — f <£o 2 (r 2 t)/(rut)dr*. Roothaan and Weiss 13 have performed variational 
calculations for the total energy of two-electron atoms using trial wave functions of 
the form (23). The extrapolation of their results for various Z to Z = a> yields 14 
T ~ -0.0440 Hartree, or 94 per cent of the exact value. The interaction be- 
tween the factor F and the orbital <p is therefore very large and beneficial, as pointed 
out by Rosen fold and Konowalow. 14 


Summary —In a recent paper Kirtman and Houston 1 derived a first-order trial 
wave function & for a perturbed Hamiltonian //„ H- V of the form & = 

1 ho object of this paper is to extend thoir idou to the total SchrGdingcr equation, 
and to uso functional minimisation to obtain tho best trial wave function of the 
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form \p = F{V)t». The resulting equation for the first-order wave function agrees 
with theirs. The theory is applied to the first-order correlation problem for two 
electron atoms, and accounts for 70 per cent of the second-order correlation energy. 

* This reseaich was supported by National Aeronautics and Space Administration grant NsG-. 
275-62 and by a National Science Foundation graduate fellowship. 

* Kirlman, B. and M. L. Bcnslon, J. Chem. Phys., 46, 472 (1967). 

2 Unsold, A., Z. Physik, 43, 5G3 (1927). 

* Ilirschfelder, J. O., W. B. Biown, and S. T. Epstein, Advan. Quantum Chem., 1, 341 (1964). 

I Parr, R. G., J. Chem. Phys., 40, 3726 (1964). 

5 Young, W. H., and N. H. March, Phys. Rev., 109, 1S54 (1958); ibid., 106, 1151 (1957). 

6 Schcrr, C. W., and R. E. Knight, Rev. Mod. Phys., 3S, 436 (1963). 

7 Midtdal, J., Phys. Rev., 138, A 1010 (1965). 

8 Kirtman and Benston’s/ is too large by a factor of 2. 
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14 Roscnfeld, J. L. J., and D. D’. Konowalow, J. Chem. Phys., 41, 3556 (1964). 
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COMMENT ON VAN DER AVOIRD’S WAVE' OPERATOR FORMALISM 

R. E. JOHNSON and P. R. CERTAIN 

University of Wisconsin Theoretical Chemistry Institute, Madison, Wisconsin, USA 
Received 23 October 1967 


In regard to the papers on the perturbation 
theory of forces by Van der Avoird [1,2], it is 
worthwhile to point out that the matrix PRgP, 
defined * by eqs. 5 and 14, does not have an in- 
verse in the space P, then the proposed solution 
will not satisfy the original Schroedinger equa- 
tion. This follows since the set of equations 

PR 0 P{E-K)PT = PR 0 P 

will not have a unique solution, and the inverse 
of PR 0 P{E~H)P, expressed as an expansion in 
eq. 19, will not exist. 

For the case when E Q is the ground state of 
H Q) the R Q is negative definite so that PR Q P has 

* Here the equation numbers and notation are those of 

ref. [lj. 


an inverse and there is no difficulty **. For ex- 
cited states R 0 is no longer a definite operator 
so that the inverse of PR Q P might not exist. Al- 
though a slight change in the definition of tf 0 
would correct for this, convergence problems 
due to this difficulty may occur in the application 
of the theory to excited states. 


REFERENCES 

[1] A. Van der Avoird, J.Chem.Phys. (1 November 
1967). 

[2] A. Van der Avoird, Chem.Phys. Letters 1 (1967) 411 

** Preliminary calculations on the ground states of 
some simple potentials have shown that the theory 
gives reasonable results. 


November 1967 
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Calculation of Matrix Elements for One-Dimensional 
Quantum-Mechanical Problems* 

A. S. Dickinson! and P. R. Certain! 

University c/ Wisconsin Theoretical Chemistry Iiislitute, Madison , Wisconsin 53706 
(Received 10 May 1968) 

. A simple method proposed by Harris el al. using the techniques of transformation theory for the 
generation of the matrix elements of one-dimensional potential functions m a discrete, orthonormal basis 
is shown to be equivalent to Gaussian quadratures when Uie basis is constructed of orthogonal polynomials. 
The basis exp (*«0) on (— v, t) is also discussed. 


I. INTRODUCTION 


where the T n j and X/ are determined by 


A method has been proposed by Harris ei ah 1 and 
employed by several authors, 2 for the calculation of the 
matrix of a one-dimensional potential function V{x) 
in a discrete, orthonormal basis <&,(&), n—Q } 1, »•*, N, 
where the set <f> n is complete for «> 4 The variable 
x has the range (a, b ). The implementation, of the 
method requires the existence of a single-valued 
function u(%) on the range (a, &), in which case their 
prescription for the evaluation, of the matrix elements 
ofFis 

= 0<n,m<N, (1) 

i-o 


* This work received financial support from the National 
Aeronautics ami Space Administration Grant NsG-275-62. 

t "Present address: Department of Applied Mathematics, The 
Queen's Univcisity of Belfast, Belfast, N. Ireland. 

| National Science Foundation Graduate Fellow. 

1 D. O. Harris, G. G. Engerholm, and W. D. Gwinn, J. Chem. 
Phy9. 43, 1515 (1965). 

* (a) D. O. Harris, H. W. Harrington, A, C. Luntz, and W. D. 
Gwirra, J. Chem. Phys. 44, 3467 (1966); (b) D. F. Zetik and F. 
A. Matsen, J. Mol. Spectry. 24, 122 (1967) ; (c) P. F. Endres, J. 
Chem. Phys. 47, 798 (1967); (d) C. Schwartz, j. Comput. 
Phys.3,90 (1967). 


r N 

<Pn (x) <i>m (x) dX— 5nm~ 23 ( 2 ) 

j — 0 

r A r 

'An (#) # (a?) (frm (x) dx— lifitn ** 23 (3) 

j— 0 

i.e., the orthogonal matrix T~{T nm ) diagonalizes 

For the case where the u is tridiagonal, Harris cl al. 
have shown that V nm is exact for F, a polynomial in 
u of precise degree or less. For the 

general case where only il&t nonvanishing, it 
may be shown that V nm is exact for V of precise degree 
(2N—n~~ m) /V-f~ i or less. Hence, it is desirable to 
employ a function u whose matrix is tridiagonal. 

Since the formula (1) is very suggestive of an (iY-j- 1) - 
point mechanical quadrature, it may be compared with 
the Gaussian -type quadrature of the same order. In the 
present paper, the relationship of (i) to Gaussian 
quadratures is considered for; (i) the set </» rt constructed 
from orthogonal polynomials on {a } b) and; (ii) <f> H ~ 
exp (in 0) on (— x, tr), 

H. ORTHOGONAL POLYNOMIALS 

If u(x) —x, and the basis <f> n is obtained from the 
first Y-f-1 polynomials in x orthogonal with respect to 
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the positive weight function w(x) on (a, b ), i.e., if 
<f >» (*) *= [>(*) /h»l ll2 f » (*) > «=0, 1, *•«, N, (4) 

where 

f U ( x )fm (x)w(x) dx= hj nm 

"a 

and 

M*)=k n x’'+k n 'x*~'+---, (5) 

then the equivalence of (1) and a Gaussian quadrature 
may be explicitly demonstrated. 

Using the Gaussian quadrature defined by w(x) on 
(a, 6), the formula for V ntn is 3 


Expanding along the last row of (12), we sec that 
the D# satisfy the recursion relation 

ZV+i (X) = Dx(h)— (X). (13) 

It is easy to check that A(X) =/i(X) and Z> 2 (X) — / 2 (X) 
so that by (13), D N {\) =/y(X) for all N>0, and hence 
(11) is the same as (7). This demonstrates state- 
ment (i) . 

To show the validity of (ii) we make use of two 
identities which are easily derived from Eq. (9). The 
first is a special case of the Christo ffel-Darboux 
identity, 4 

E K-y.Mf.M = w, (i4) 

»-0 

and the second is 


where the are the zeros offx+i, 


(6) 

Efe r%MM<b)AA -' — swwr 1 


/w(a/">) =0, j= 0, 1, •••, N (7) 
and 

W A,, =-WW^ +1 W(% w )/w +2 (e^>), (8) 

where f{x)~df/dx. In the following the explicit 
dependence of the a/ iV) and upon N will be 

suppressed. 

To establish the identity of (1) and (6) it is necessary 
to prove the following relations: 

(i) X,[Eq. (3)]=a,£Eq. (7)],'* wher^ both are 
arranged monotonically; 

(ii) T^iWj/hy^iaj). 


, fn+lMIn (aj) +/„ («,•)/„+! (a,-) 


where a {) as are zeros of fx + i. 

Let S„y= (Ws/h n ) 1/2 /„ (a,-) . Then (14) implies that 

t SA=s,i, ( 16 ) 


and (15) implies that 
n tr 

T"/ SniSmk%nm~ XjSjfc. (I/) 


The statement (i) follows from the recursion rela- 
tion 4 among the /„ 

/«+i= (A n +B n x)f n ~C n f n - h (9) 

where 

Bn — ^n-fl/ knj 

Cn=BJt n /B n ^h n ^ (10) 

This greatly simplifies the construction of the secular 
equation to determine the X/, 

Aw(Xy) = 0, . 0, 1, • • • , N> (11) 

where 

ZMX)==(-)*7wdet(u-Xl), (12) 

and u and 1 arc of dimension NXN. 


* A. Ralston, A First Course in Numerical Analysis (McGraw- 
Hill Book Co., New York, 1965), pp. 85—111. 

4 G. Szeg6, Orthogonal Polynomials (American Mathematical 
Society Colloquium Publications, New York, 1959), Vol. 23, 
Sec. 3.2. 


Hence, the S n s are the elements of an orthogonal 
matrix which diagonalizes x. Since the eigenvalues of 
the latter matrix are nondegenerate, 3 the elements S n] 
are unique and thus equal to the T»s» 

Thus, since (1) is a Gaussian quadmture when tin: 
<j) n aic obtained from (4) and (5), V nm will be exact 
for V } a polynomial of precise degree p, where 
P< (ZN+l—n—m), as was shown by Harris ct al . 

For many of the basis sets satisfying (4) and (5) 
likely to be employed in problems of physical interest 
there exist corresponding quadratures with published 
values 6 of the aj, which substantially reduces the 
effort required to diagonalize x. 

III. THEJ3ASIS exp (htO) 

Where the’ basis functions <f> 0 do not satisfy (4) and 
(5), the transformation method may again be equiva- 
lent to a Gaussian quadrature, but not necessarily to 
the most efficient choice of weight function w(x). The 


1 A. H. Stroud and D. Sccrest, Gaussian Quadrature Formulas 
(Prentice-Hall, Inc., Englewood Cliffs, N.J., 1966). 
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basis expp(»— iV) 03 , —Tv<0<ir, «= 0, 1, •••, 2 N, 
discussed by Harris el al . may be taken as an example. 

It is assumed that V is a polynomial of precise 
degree p in u— cos 0, Since u is tridiagonal in this basis, 
the characteristic polynomials Djt, defined by (12), 
satisfy a three-term recursion relation. This is easily 
identified as that satisfied by the Chebyshev poly- 
nomials of the second kind, 8 U n (%). It should be noted 
that the zeros of U n can be written in closed form, and 
T can be obtained from statement (ii) above. There- 
fore, by (6) the matrix transformation technique 
evaluates the (2iV-fl) -point Gauss-Chebyshev quad- 
rature (of the second kind) of the integral 

- r V (co$8) U * (cosf?) U m (cosfl) sin toft, 0 <», m< 2 N, 
J Q 

which is exact for /><4iY+l— On the other 
hand, using the relations between the U n and the 
Chebyshev polynomials of the first kind T „ , V nm may 
be written 

7„ m = - f V(cosd) U n (cos0) Um(cosd) sin 7 8J9 
W o 

+*- 1 [* V(cosd)T n+mi . 2 (co$d)dd. 

The second integral on the right vanishes for p <n-\-m + 2 
by orthogonality. The two conditions on p for V nm to 
be exact may be combined to p<2N-fl— j 2 N—n—m |, 
for a basis consisting of 2iV-bl functions. 

The connection between the transformation method 
and Gaussian quadratures for this case may be clarified 
by applying the Wang transformation. 7 This is equiva- 
lent to choosing as a (2A 7 -f l)-term basis, the functions 

r n (cos0), 0<»<iV 


• Handbook of Mathematical Functions, M. Abramowitz and I. 
A. Stegun, Eds. (National Bureau of Standards, Washington, 
D.C., 1964), Chaps. 22 and 25 The Z>, distinguishes between the 
Chebyshev polynomials of the first kind. r n (costf) = cos«0, and of 
the second kind, f7 n _i(cos0) = (swHtfj/sinff, which otherwise 
satisfy the same recurrence relation. 

7 J. E, Wollrab, Rotational Spectra and Molecular Structure 
(Academic Press Inc., New York, 1967), p. 26. We are indebted 
to the referee for bringing this to our attention. 


and 

sin0Z7 n (cos0), 0<tt<iV— 1. 

These are, respectively, symmetric and antisymmetric 
with respect to 0— >— 0 and hence, the matrix V factors 
into an (N+ 1)X(AT+1) block and an NxN block. 
Since the basis now consists of orthogonal polynomials, 
the analysis of the preceding section may be applied 
separately to each block. This leads to the same condi- 
tions on p for V nm to be exact as given in the previous 
paragraph. 

IV. CONCLUSIONS 

The matrix transformation method provides a con- 
venient technique for generating integrals for a one- 
dimensional variational problem using standard matrix 
manipulations. As Harris el al. point out, the diago- 
nalization of u need only be performed once for a given 
N , if the X; and T are retained from problem to problem. 
In addition, for many cases of physical interest, a 
scaling parameter may be introduced into the basis*, 
enabling several different basis sets to employ the same 
Ay and T. For example, in the work of Zetik and Matsen 2b 
on the computation of vibrational-rotational energy 
levels, if the matrix of u= {R—Rd) where 

the parameters have their usual significance, is com- 
puted, then the values of the independent variable R 
may be obtained by a suitable choice of Ro and k . 

The method does, however, add an unnecessary 
inflexibility to the energy eigenvalue problem in that 
evaluation of the matrix elements and the convergence 
of the eigenvalues with increasing N need not be con- 
sidered simultaneously. The minimum number of 
values of V required for the accurate evaluation of V nm 
is not independent of n and m and need not equal the 
number of basis functions required for a satisfactory 
representation of the wavefunction. 

Note added in proof: The properties of the matrix u 
of Sec. II are also discussed by Wilf. 8 See also Golub and 
Welsh. 9 


8 H. S. Wilf, Mathematics for the Physical Sciences (John Wiley 
& Sons, Inc., New York, 1962), p. 55. 

9 G. II. Golub and J. H. Welscn, Calculation of Gauss Quadra- 
ture Rules, Computer Science Department Technical Reference 
CS81, Stanford University, Stanford, Calif., 1967 (unpublished). 
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Localized orbitals are defined for arbitrary wave functions for atoms- or molecules with, even or odd 
numbers of electrons by making a particular non-orthogonal transformation among the natural orbitals 
which preserves the diagonal nature of the electron density. Different definitions of localized orbitals 
are obtained, by relating them to different approximate forms for the wave function. The equilateral tri- 
angle configuration of Hg is treated as an example. 


1. INTRODUCTION 

Up to- now, the concept of localized orbitals 
[l-3],for atoms and molecules has been limited- 
to systems described by single- de ter minantal 
wave function. In the present paper we explore 
ways of extending this to molecules with odd or 
even numbers of electrons described by arbitrary 
wave- function. The extention is not unique- and is 
guided' by relating localized orbitals to particular 
approximate forms for the wave function. Mini- 
mizing the- energy obtained from a Hartree pro- 
duct of localized orbitals leads to the Edmiston- 
Ruedenberg criterion [2] for localization. Rela- 
ting' the localized orbitals to a valence bond func- 
tion leads to a different criterion. These two pos- 
sibilities are illustrated by considering the equi- 
lateral triangle configuration of Hg. 

The concept of localized orbitals was first de- 
veloped extensively by Lennard-Jones,. Pople, 
and Hall fl] for a system with high symmetry, 
whose wave function is- a single determinant of 
molecular orbitals. These workers exploited the 
freedom to make a unitary transformation among 
the 1 molecular orbitals' without affecting the value 
of the wave function, and defined "equivalent or- 
bitals" as linear combinations of the molecular or- 
bitals which have the property of being identical 
to each other except for orientation and position 
in. space. 

Edmiston and Ruedcnberg [2] and Foster and 
Boys [3] freed the concept of localized orbitals 

* This work received’ financial support from the Na- 
tional Aeronautics and Space Administration Grant 
NsG-275-62. 

** National Science Foundation Graduate Fellow. 


from the restriction that the system be symmet- 
ric, although they still required a. single deter- 
minant wave function. The technique- used was to 
define a function of the transformation from mo- 
lecular orbitals to localized orbitals, whose ex- 
treme points fix the transformation. This function 
provides the effective definition of "localization”. 


2. TRANSFORMATION TO LOCALIZED 
ORBITALS 

We wish now to free the concept from any re- 
striction on the form 'of the wave function. This 
we do by making use of the natural orbital expan- 
sion of the electron density p(r): 

p(r) =/ j//*p6(r- r^dr, 

= S\. x ?(r).. 

The Xi and Xf denote the ith natural orbital and 
its- occupation number,, respectively. 

When xf/ is a single' determinant of doubly- 
occupied orbital's, all X£ equal two and the form 
of p is invariant to an- orthogomtl transformation 
among the x*. * n more general case, the form 
of p is- invariant to the non-orthogonal transfor- 
mation 

n i u i-f T v4*f 

where T is unitary and' normalizes Then 

p(r) =Z/wu M |(r). 
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An essential feature of the scheme we are consid- 
ering is the preservation of the form of p as a 
"sum of squares" without any cross terms. The 
reason for this will become apparent in the ex- 
ample discussed below. 

The transformation 7 is thus far unspecified. 
We choose T by a straightforward generalization 
of the Edmiston-Ruedenberg criterion [2] for the 
single determinant case. Let Vij be the electro- 
static interaction energy of the charge distribu- 
tions » , and «.//?. 

3 3 r «f(r)u?(r') 


■ dr dr’ , 


The total self-energy of p is invariant to 7, but 
is equal to the sum of two terms, neither of which 
is invariant. 


i J dr dr' v if + £ Vji 

J | r-r'| . ti £<j y 


{invariant) 


i<j 

(maximize) (minimize) 


We choose 7 to minimize the interaction sum or, 
equivalently, to maximize the self-energy sum. 

The significance of this criterion may be seen 
by relating the localized orbitals u{ to an approx- 
imate wave function constructed from them. In 
particular, for a system of 2 N electrons, a Hart- 
ree (rather than a Hartree-Fock)~like wave func- 
tion may be constructed by doubly-occupying the 
first N natural orbitals. It may be shown that mi- 
nimizing the energy of this wave function is equi- 
valent to the Edmiston-Ruedenberg criterion for 
this case. 

This suggests that another way to obtain loca- 
lized orbitals is to relate them to a more accu- 
rate form for the wave function, such as a va- 
lence-bond function, or at least one which satis- 
fies the Pauli principle. We shall explore this 
possibility by means of the following example. 


3. EXAMPLE: Hg (equilateral triangle) 

We wish to consider the electron distribution 
of H$ in the equilateral triangle (R = 1.6575 Bohr) 
conngurnlion. The natural orbitals have recently 
been given by Christoff orsen and Shull [4j, based 
on Christoffersen's twelve configuration wave 
function [5], Truncating the natural orbital ex- 
pansion after three terms, we can write the wave 
function as 

«1,2) = x!xid)Xi(2) - Xf[x 2 (l)X 2 (2) + X 3 (1)X 3 (2)J , 
where the representation of the x* in terms of 
Slater orbitals is given in ref. [4] where = 


= 0,993594 andX 2 = 0.079906. The xi belongs to 
the totally symmetric representation of 03 ^, 
whereas X 2 ami X 3 belong to the doubly-degener- 
ate representation e\ 

3.1. Localized orbitals for H 3 

Now we can transform to three localized orbi- 
tals ?< a> « c by making use of the high sym- 
metry. Since the natural orbitals form a basis 
for an irreducible representation of £> 3 ^, the lo- 
calized orbitals form a basis for a reducible re- 
presentation. If C| denotes a counterclockwise 
rotation of 120 ° about the vertical axis, and cr a 
denotes reflection through the vertical plane con- 
taining nucleus a , then a reasonable choice for 
the symmetry properties of u a , u^, u c is 

u b = c Z u a > u c - c 3 u b * °a v a = u a • 

The 7 is the similarity transformation which ■* 
connects the representation spanned by xj, X 2 > 

X 3 with that spanned by u a> u c . Because of 
the degeneracy of e*, group theory does not com- 
pletely fix 7, but rather gives two possible sets 
of localized orbitals, which may be generated from 

u a = x l x l ± ( X l/ 2 )V3X 2 - X 3 ) (1) 

by application of C 3 . The plus sign corresjjonds 
to an atomic density polarized to a large extent 
toward the other atoms. The minus sign corres- 
ponds to a bond density polarized toward the non- 
bonded atom. The orbital which is the more loca- 
lized, i.e. which maximizes 

v aa m f k - r T 1 dr df 

is Its representation in terms of Slater orbi- 
tals is given in table 1 . 

3.2. Valence bond form 

The preceding analysis characterizes the to- 
tal density as a sum of three localized densities. 

A different sort of analysis is suggested by re- 
lating the natural orbital expansion of the wave 
function to a valence bond form, namely, 

, where 

<b ah = «(W 2) + M1M2) , 

(\f + 4\| fa = Xjxi + x 2 ^ x 2 " ♦ ( 2 ) 

We have simply rearranged the three term natu- 
ral orbital expansion of the wave function to have 
this form, and have required a, b and c to have the 
same symmetry properties as previously as- 
signed to a a , Uty u c . The orbital a is an atomic 
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Table 1 

Expansions of the orbitals *4, a, a a iy in terms of Slater 
orbitals for Hjj (equilateral triangle) 



Slater orbitals 

< 

(eq. 1) 

a 

(eq. 2) 

a ab 

(eq. 3) 


exp(-Cxr a ) 

+0.638 

+1.169 

+0.932 

Is 

expt-Cjrj,) 

+0.376 

+0.023 

-0.067 


oxp{-?ir c ) 

+0.376 

+0.023 

+0.433 


>a cxpt-Can*) 

-0.115 

-0.164 

-0.142 

2s 

>’b oxj >K 2 »'6) 

-0.088 

-0.046 

-0.040 


£|A)i >* c expf-^r^ 

-0.088 

-0.046 

-0.091 


^aexp(-C3 y fl) 

+0.067 

+0.140 

0.107 

2p 


+0.032 

-0.012 

-0.025 


(Cj>/J7p Z c exp(-? 3 r c ) 

+0.032 

-0.012 

0.041 


Cj = 1.20, £ 2 = 1-175, C 3 = 1.80. The p orbitals are di- 
rected from each nucleus to the center of the triangle. 
Thus Z a is dn'ected from nucleus a to the center of the 
triangle. 


proach, namely, although the individual terms 
and 4>' a b are quite different, the total density 
given by the two functional forms is the same, 
due to important contributions from "overlap 
densities". Preserving the form of p as we have 
done in defining localized orbitals avoids this 
difficulty. 

3.4. Comparison 

Finally, let us compare the three orbitals u (t , 
a, and a a f } by giving their expansion in terms of 
Slater orbitals centered on Lhc three nuclei. 11 is 
clear that a is the most localized having almost 
no contribution from the orbitals of the oilier nu- 
clei, although the meaning attached to this is dif- 
ferent from that attached to u+. 

Thus we see that .by relating the localized or- 
bitals to different functional forms for approxi- 
mate wave functions, we can arrive at different 
definitions for the localized orbitals. The inter- 
pretation given the localized orbitals then follows 
from the form of the wave function. 


orbital polarized to a slight extent toward the oth- 
er nuclei. Its representation in terms of Slater 
orbitals is given in table 1. 

3.3. Christoffersen and Shull 

Christoff ersen and Shull [4] have also written 
the wave function in a valence bond form. How- 
ever, they took the orbital a in to be differ- 
ent from a in <p ac . They wrote 

*(1.2) = JVW& + *ic+flic]. 

where 

(A 1 + 2A 2 ) \* " 4 X 1 + < A 2 /l ^)(X 2 -/3x 3 ). (3) 
The representation of a a b in terms of Slater or- 
bitals is given in table 1. Their result was a 
very delocalized description. The total wave func- 
tion is the same as ours, however, and we see 
no advantage to allowing this added flexibility. It 
does point up an unattractive feature of this ap- 
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